线性反馈移位寄存器-LFSR

LFSR

badmonkey的博客
最开始了解到LFSR的时候是在学习MT19937伪随机数生成器的时候，当时也是初步了解。也没有用代码实现过，最近做了几道相关的题，在这里记录一下。

原理

证明过程不太会，简单说一下大致的原理。LFSR是用于生成随机数序列的，每一个LFSR都是有级数的，级数决定了潜在的循环周期（即随机数序列的周期）。对于每个LFSR还需要有一个线性函数，用于生成随机数。比如

四级LFSR	参数
级数	4
最长循环周期	2^4-1
线性函数	x4+x2+x+1

这里的线性函数必须是二元域上的n阶不可约多项式，线性函数的作用有点类似掩码的功能，对上面的线性函数相当于mask 1011（x的系数非零则为1，不考虑常数项，从左到右次数依次递减）

---

对上面的例子，选取一个4位的seed用于生成随机数，假设seed为1010
step1
线性函数充当掩码的功能选取seed的对应位，即 seed & mask => 1010 & 1011 选取了第1，3，4位分别为 1，1，0 ，则有 1 xor 1 xor 0 => 0
step2
将上一步得到的结果和seed后三位进行拼接，所以新的随机数就是0100

---

同理下一个随机数为1000，若经过k次迭代后生成了和seed一样的的数，则称k为为循环周期，显然周期越长为好，MT19937的周期就很长可以达到2^19937-1。

题目

CISCN 2018 oldstreamgame

# coding=utf-8
flag = "flag{xxxxxxxxxxxxxxxx}"
assert flag.startswith("flag{")
assert flag.endswith("}")
assert len(flag)==14

def lfsr(R,mask):
    output = (R << 1) & 0xffffffff
    i=(R&mask)&0xffffffff
    lastbit=0
    while i!=0:
        lastbit^=(i&1)
        i=i>>1
    output^=lastbit
    return (output,lastbit)
R=int(flag[5:-1],16)
mask = 0b10100100000010000000100010010100

f=open("key","w")
for i in range(100):
    tmp=0
    for j in range(8):
        (R,out)=lfsr(R,mask)
        tmp=(tmp << 1)^out
    f.write(chr(tmp))
f.close()

相信看过原理之后应该可以理解代码，关键代码为

for i in range(100):
    tmp=0
    for j in range(8):
        (R,out)=lfsr(R,mask)
        tmp=(tmp << 1)^out
    f.write(chr(tmp))</