标签：ch int 鳄鱼 base P2579 ans isdigit ZJOI2005 getchar

https://www.luogu.org/problemnew/show/P2579

1、注意正确把握周期时间点，从A[1]乘到A[11]再乘A[0]

2、矩阵乘法没有交换律，得先预处理Q=A[1]*A[2]*......*A[11]*A[0]，再用算出Q^[K/12]，对于剩下的，从A[1]开始乘完K%12个

#include<bits/stdc++.h>
#define LL long long
using namespace std;
int read(){
    char ch=getchar(); int w=1,c=0;
    for (;!isdigit(ch);ch=getchar()) if (ch=='-') w=-1;
    for (;isdigit(ch);ch=getchar()) c=(c<<3)+(c<<1)+(ch^48);
    return w*c;
}
const int mod=10000;
int n,m,S,E,K,P;
struct Matrix{
    int a[51][51];
    Matrix (){memset(a,0,sizeof a); }
    Matrix operator *(const Matrix &b) const {
        Matrix ret;
        for (int i=0;i<n;i++)
        for (int j=0;j<n;j++)
        for (int k=0;k<n;k++)
        ret.a[i][j]=(ret.a[i][j]+a[i][k]*b.a[k][j])%mod;
        return ret;
    }
}g[12],base,o,ans;
int main(){
    n=read(); m=read(); S=read(); E=read(); K=read();
    for (int i=1;i<=m;i++){
        int x=read(),y=read();
        o.a[x][y]=o.a[y][x]=1;
    }
    for (int i=0;i<12;i++) g[i]=o;
    P=read();
    while (P--){
        int T=read();
        for (int i=0;i<T;i++){
            int x=read();
            for (int j=i;j<12;j+=T){
                for (int k=0;k<n;k++)
                g[j].a[k][x]=0;
            }
        }
    }
    base=g[1];
    for (int i=2;i<12;i++) base=base*g[i];
    base=base*g[0];
    for (int i=0;i<n;i++) ans.a[i][i]=1; 
    int nn=K/12;
    for (;nn;nn>>=1,base=base*base)
    if (nn&1) ans=ans*base;
    for (int i=1;i<=K%12;i++) ans=ans*g[i];
    cout<<ans.a[S][E]<<"\n";
    return 0;
}

标签：ch,int,鳄鱼,base,P2579,ans,isdigit,ZJOI2005,getchar
来源： https://www.cnblogs.com/yuyue2005/p/10802145.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。