clc;
clear
a=[001110];
syms?x;
fx=[];%多項式
fx（1）=1;
lc=[];%L
constantpoly=[1];
lc（1）=0;
sequenlength=length（a）;??%?序列長度
minpolycoefficient=[1];
minpolylen=length（minpolycoefficient）;
minpolydegree=minpolylen-1;
d=[];
count1=0;
count2=0;
for?n=0:sequenlength-1
????dtemp=0;
????for?j=0:minpolydegree
????????dtemp=mod（dtemp+minpolycoefficient（j+1）*a（n+1-minpolydegree+j）2）;
????end
????d=[d?dtemp];
????if?d（n+1）==0
????????fxtemp=fx（n+1）;
????????lctemp=lc（n+1）;
????else
????????for?i=0:n
????????????if?lc（i+1）==0
????????????????count1=count1+1;
????????????else
????????????????count1=count1;
????????????end
????????end
????????keyscript=[0];
????????if?count1==n+1
????????????fxtemp=1+x^（n+1）;
????????????lctemp=n+1;
????????????keyscript=[keyscript?n];
????????else
????????????for?k=0:（n-1）
????????????????if?lc（k+1）==lc（k+2）
????????????????????keyscript=keyscript;
????????????????else
????????????????????keyscript=[keyscript?k];
????????????????end
????????????end
????????????keyscript=[keyscript?n];
????????????keyscriptlength=length（keyscript）;
????????????ntemp=keyscript（keyscriptlength）;
????????????mtemp=keyscript（keyscriptlength-1）;
????????????if?mtemp==0
????????????????fxtemp=fx（ntemp+1）+x^ntemp;
????????????else
????????????????fxmcoefficient=sym2poly（fx（mtemp+1））;
????????????????fxmlen=length（fxmcoefficient）;
????????????????xnmcoefficient=sym2poly（x^（ntemp-mtemp））;
???