FWT板子

来源:互联网 发布:2016华语电影知乎 编辑:程序博客网 时间:2024/05/21 06:44

感觉这东西用的不算多。所以先把板子记下来吧

void FWT(LL *a,int n){    for (int i=1;i<n;i<<=1)     for (int p1=i<<1,j=0;j<n;j+=p1)      for (int k=0;k<i;k++) {        LL x=a[j+k]; LL y=a[j+k+i];        a[j+k]=(x+y)%p;         a[j+k+i]=(x-y+p)%p;        //xor: a[j+k]=x+y,a[j+k+i]=x-y        //and: a[j+k]=x+y        //or : a[j+k+i]=x+y      }}void UFWT(LL *a,int n){    for (int i=1;i<n;i<<=1)     for (int p1=i<<1,j=0;j<n;j+=p1)      for (int k=0;k<i;k++){        LL x=a[j+k]; LL y=a[j+k+i];        a[j+k]=(x+y)%p*ret%p;        a[j+k+i]=((x-y)*ret%p+p)%p;        //xor: a[j+k]=(x+y)/2,a[j+k+i]=(x-y)/2        //and: a[j+k]=x-y        //or : a[j+k+i]=y-x      }}void solve(LL *a,LL *b,int n){    FWT(a,n); FWT(b,n);    for (int i=0;i<n;i++) a[i]=a[i]*b[i]%p;    UFWT(a,n);}