( VIJOS )VOJ 1049 送给圣诞夜的礼品 矩阵快速幂

https://vijos.org/p/1049

 

非常普通的矩阵快速幂...

但是我

第一次写忘了矩阵不能交换律...

第一二次提交RE直到看到题解才发现这道题不能用递归快速幂...

第三次提交成了c编译错误...

第四次提交WA发现写循环快速幂的时候少清零了一个f...

所以提交了五次才终于对了,什么垃圾的代码能力...通过率杀手...

希望下次能捡起我的脑子....

代码很简单 如下

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #include<cmath>
 #include<queue>
 using namespace std;
 const int maxn=;
 const double eps=1e-;
 int n,m;
 struct mat{
     int e[maxn][maxn];
 };
 mat pro(mat x,mat y){
     mat z;
     for(int i=;i<=n;i++){
         for(int j=;j<=n;j++){
             z.e[i][j]=;
             for(int k=;k<=n;k++){
                 z.e[i][j]+=x.e[i][k]*y.e[k][j];
             }
         }
     }return z;
 }
 mat pow(mat x,int k){
     mat z;
     bool f=;
     while(k>){
         if(k%!=){
             if(f){
                 z=x;
                 f=;
             }else{
                 z=pro(x,z);
             }
         }k/=;
         x=pro(x,x);
     }
     return z;
 }
 int main(){
     int k;
     scanf("%d%d%d",&n,&m,&k);
     int x,z,w;
     z=k/m; w=k%m; mat c;
     for(int i=;i<=n;i++){
         for(int j=;j<=n;j++){
             c.e[i][j]=;
         }
     }
     mat a=c; mat d=c;
     for(int i=;i<=m;i++){
         mat b=c;
         if(i==){
             for(int j=;j<=n;j++){
                 scanf("%d",&x);
                 a.e[j][x]=;
             }
             if(w>=){
                 d=a;
             }
             continue;
         }
         for(int j=;j<=n;j++){
             scanf("%d",&x);
             b.e[j][x]=;
         }
         a=pro(b,a);
         if(i<=w){
             d=a;
         }
     }
     if(z==){
         a=d;
     }else if(w==){
         c=a;
         a=pow(c,z);
     }else{
         c=a;
         a=pro(d,pow(c,z));
     }
     for(int i=;i<=n;i++){
         for(int j=;j<=n;j++){
             if(a.e[i][j]!=){
                 printf("%d ",j);
                 break;
             }
         }
     }
     cout<<endl;
     return ;
 }