hdu 2157 How many ways??(矩阵经典应用)

题目链接：hdu 2157 How many ways??

题意：

给你一个邻接矩阵，问你从a到b点，可走重复路径，恰好经过k个节点的路径数。

题解：

其实就是一个裸的矩阵快速幂，每走一次，相当于一次乘法。

#include<bits/stdc++.h>
#define mst(a,b) memset(a,b,sizeof(a))
#define F(i,a,b) for(int i=a;i<=b;++i)
using namespace std;
typedef long long ll;

const int mat_N=21;
int mo=1000,n,m;
struct mat{
    ll c[mat_N][mat_N];
    void init(){mst(c,0);}
    mat operator*(mat b){
        mat M;int N=n-1;M.init();
        F(i,0,N)F(j,0,N)F(k,0,N)M.c[i][k]=(M.c[i][k]+c[i][j]*b.c[j][k])%mo;
        return M;
    }
    mat operator^(ll k){
        mat ans,M=(*this);int N=n-1;ans.init();
        F(i,0,N)ans.c[i][i]=1;
        while(k){if(k&1)ans=ans*M;k>>=1,M=M*M;}
        return ans;
    }
}A,B;

int main()
{
    while(scanf("%d%d",&n,&m),n+m)
    {
        A.init();
        int x,y,t;
        F(i,1,m)
        {
            scanf("%d%d",&x,&y);
            A.c[x][y]=1;
        }
        scanf("%d",&t);
        while(t--)
        {
            int a,b,k;
            scanf("%d%d%d",&a,&b,&k);
            B=A^k;
            printf("%lld
",B.c[a][b]);
        }
    }
    return 0;
}