c白色15.fc 2009年1月15日按顺序读取,并通过重复追加扩展c卷曲数。c与Gijswijt的序列A090822有关,更直接地与A094004有关。c用于生成A116909,方法是给它一个包含两行的输入文件c 22 500c 232232322323223223223转2323223c参见笔记本Latt 94 p 97c隐式整数(a-z)整数a(100000)c读取起始序列的长度和所需的字数读取(*,*)M0,M1c读取初始序列读取(*,200)(a(i),i=1,M0)200格式(60i1)写入(*,*)“开始”写入(*,*)“M0,M1=”,M0,M2写入(*,200)(a(i),i=1,M0)100格式(30i2)c获得下学期do 1 n=M0+1,M1c写入(*,*)“下一个n=”,nkrec=1yrec=1c try next value of ly=y的长度do 2 ly=1,n-1如果((krec+1)*ly.gt.n-1)转到7c写(*,*)“trying ly=”,lyc计算这个ly的kdo 3 k=2,M1do 4 i=1,y什么=n-(k-1)*ly-i如果(a(what).ne.a(n-i))转到54继续c到达了匹配行的末尾3继续c永远不会到这里写入(*,*)“错误#3”停止5继续k=k-1c写入(*,*)“k=”,k如果(k.le.krec)转到6c k的新记录krec=kyrec=lyc写(*,*)“将k增加到”,krecc写(*,*)“with ly=”,ly6继续2继续7继续c如果krec=1,则退出c如果(krec.eq.1)转到299a(n)=krec1继续写入(*,100)(a(i),i=1,M1)do 30 n=150030写入(*,31)n,a(n)31格式(i6、i4)停止结束