数学专项matrix：UVa 11551

快速幂模板题。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=60;
const int mod=1000;
typedef int matrix[maxn][maxn];
int a[maxn],b[maxn];
int n,r;
void mat_mul(matrix A,matrix B,matrix res)
{
    matrix C;
    memset(C,0,sizeof(C));
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            for(int k=0;k<n;k++)
                C[i][j]=(C[i][j]+A[i][k]*B[k][j])%mod;
    memcpy(res,C,sizeof(C));
}
void mat_pow(matrix A,int m,matrix res)
{
    matrix a,r;
    memset(r,0,sizeof(r));
    memcpy(a,A,sizeof(a));
    for(int i=0;i<n;i++) r[i][i]=1;
    while(m)
    {
        if(m&1) mat_mul(r,a,r);
        m>>=1;
        mat_mul(a,a,a);
    }
    memcpy(res,r,sizeof(r));
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&r);
        for(int i=0;i<n;i++) scanf("%d",&a[i]);
        matrix A;
        memset(A,0,sizeof(A));
        for(int i=0;i<n;i++)
        {
            int t;
            scanf("%d",&t);
            for(int j=0;j<t;j++)
            {
                int p;
                scanf("%d",&p);
                A[i][p]=1;
            }
        }
        mat_pow(A,r,A);
        memset(b,0,sizeof(b));
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                b[i]=(b[i]+A[i][j]*a[j])%mod;
        for(int i=0;i<n;i++)
            if(i) printf(" %d",b[i]);
            else printf("%d",b[i]);
        puts("");
    }
    return 0;
}