HDU 4887 Endless Punishment (BSGS)

题意是给出S1，S2两个序列，对应下方01串的下标，每次操作为：查询S1内对应的01串的数字和，若是奇数则将01串第一位去掉前移，在最后添加一个1，反之添加一个0.

最后给出起始01串和最终01串，问你最少几次可以变化完成。

可以对操作写出一个转移矩阵，假设起始01串矩阵A，终点01串矩阵B，转移矩阵P,操作次数x；

写作，A*P^x=B, 这里我们令x=i*m-j ，   m=sqrt(所有可能状态数),   可以表示为A*P^(i*m)=B*P^j

那么我们只要将0~m-1的j存入到一个map，枚举i查找map中是否有对应的P^j使得其相等。

#include <cstdio>
#include <map>
#include <cmath>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define ll __int64
int s1[35], s2[35];
int seq1[35], seq2[35], n;
struct MAT1
{
    int a[35][35];
} mt1;
struct MAT2
{
    int b[35];
} mt2, x, y;
MAT2 mul(MAT1 r1,MAT2 r2)
{
    int i, j;
    MAT2 tp;
    memset(tp.b,0,sizeof(tp.b));
    for(i=0;i<=n;i++)
        for(j=0;j<=n;j++) tp.b[i]+=r1.a[i][j]*r2.b[j];
    for(i=0;i<=n;i++) tp.b[i]%=2;
    return tp;
}
MAT1 fast_pow(MAT1 r1,int cnt)
{
    MAT1 ret, tp;
    int i, j, k;
    memset(ret.a,0,sizeof(ret.a));
    for(i=0;i<=n;i++) ret.a[i][i]=1;
    while(cnt){
        if(cnt&1){
            memset(tp.a,0,sizeof(tp.a));
            for(i=0;i<=n;i++){
                for(j=0;j<=n;j++){
                    for(k=0;k<=n;k++){
                        tp.a[i][j]+=ret.a[i][k]*mt1.a[k][j];
                    }
                    tp.a[i][j]%=2;
                }
            }
            ret=tp;
        }
        memset(tp.a,0,sizeof(tp.a));
        for(i=0;i<=n;i++){
            for(j=0;j<=n;j++){
                for(k=0;k<=n;k++){
                    tp.a[i][j]+=mt1.a[i][k]*mt1.a[k][j];
                }
                tp.a[i][j]%=2;
            }
        }
        mt1=tp;
        cnt>>=1;
    }
    return ret;
}
map<int,int> mp;
int main()
{
    int i, j, S1, S2, m;
    ll tot, sta, ans;
    while(~scanf("%d%d%d",&n,&S1,&S2)){
        mp.clear();
        for(i=1;i<=S1;i++) scanf("%d",&s1[i]);
        for(i=1;i<=S2;i++) scanf("%d",&s2[i]);
        for(i=1;i<=n;i++) scanf("%d",&seq1[i]);
        for(i=1;i<=n;i++) scanf("%d",&seq2[i]);
        for(i=1;i<=n;i++) {
            if(seq1[i]!=seq2[i]) break;
        }
        if(i==n+1) {
            puts("0");
            continue;
        }
        memset(mt1.a,0,sizeof(mt1.a));
        for(i=0;i<n-1;i++) mt1.a[i][i+1]=1;
        for(i=1;i<=S1;i++) mt1.a[n-1][s1[i]-1]=1;
        for(i=1;i<=S2;i++) mt1.a[s2[i]-1][n]=1;
        mt1.a[n][n]=1;
        for(i=1;i<=n;i++) x.b[i-1]=seq1[i];
        for(i=1;i<=n;i++) y.b[i-1]=seq2[i];
        x.b[n]=1;
        y.b[n]=1;
        tot=1LL<<n;
        m=ceil(sqrt((double)tot))+1;
        sta=0;
        for(i=0;i<n;i++) {
            sta*=2;
            sta+=y.b[i];
        }
        mp[sta]=0;
        for(i=1;i<m;i++){
            y=mul(mt1,y);
            sta=0;
            for(j=0;j<n;j++) sta=sta*2+y.b[j];
            mp[sta]=i;
        }

        mt1=fast_pow(mt1,m);
        ans=-1;
        for(i=1;i<=m;i++){
            x=mul(mt1,x);
            sta=0;
            for(j=0;j<n;j++) sta=sta*2+x.b[j];
            if(mp.find(sta)!=mp.end()){
                ans=(ll)i*m-mp[sta];
                break;
            }
        }
        if(ans==-1){
            puts("poor sisyphus");
            continue;
        }
        printf("%I64d\n",ans);
    }
    return 0;
}