vijos1049 送给圣诞夜的礼品(矩阵倍增)

用01矩阵模拟位置的交换,m个操作为一组,倍增(k/m)次。remain的(k%m)次,单独做,最后乘起来。

#include 
#include 
#define N 105
int n,m,tot;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
struct Matrix{
    int mat[N][N];
    Matrix(bool t){
        memset(mat,0,sizeof(mat));
        if(t) for(int i=1;i<=n;++i) mat[i][i]=1;
    }
    Matrix operator*(Matrix b){
        Matrix res(0);
        for(int i=1;i<=n;++i)
            for(int j=1;j<=n;++j)
                for(int k=1;k<=n;++k)
                    res.mat[i][j]+=mat[i][k]*b.mat[k][j];
        return res;
    }
    Matrix operator^(int k){
        Matrix res(1),base(0);
        memcpy(base.mat,mat,sizeof(mat));
        for(;k;k>>=1,base=base*base)
            if(k&1) res=res*base;
        return res;
    }
}ans(0);
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();tot=read();
    Matrix a(1),tmp(1);
    int remain=tot%m;tot/=m;
    for(int i=1;i<=n;++i) ans.mat[1][i]=i;
    for(int i=1;i<=m;++i){
        Matrix t(0);
        for(int j=1;j<=n;++j){
            int x=read();
            t.mat[x][j]=1;
        }
        a=a*t;
        if(i==remain) tmp=a;
    }
    ans=ans*((a^tot)*tmp);
    for(int i=1;iprintf("%d ",ans.mat[1][i]);
    printf("%d\n",ans.mat[1][n]);
    return 0;
}

你可能感兴趣的:(其他oj,矩阵,模拟,快速幂)