矩阵快速幂<求有向图路径为k的方案（每条边可以取多次）>

刚开始把建立一个矩阵，用以表示任意两个顶点之间是否有边，如果有矩阵上就为1，反之为0。

那么此时如果 这个矩阵乘这个矩阵，意思就成了这个矩阵u到w长度为1的个数乘上w到v长度为1的个数，也就成了长度为2的个数的多少（边取得任意多次）。

此时得到的k=2的矩阵，这个矩阵乘长度为1的矩阵还是这个矩阵u到w长度为2的个数乘上w到v长度为1的个数，也就是长度为3的矩阵个数

这么乘可以用快速幂求出。

#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
const int MOD=1e4+7;
struct ttt{
    int map1[20][20];
};
int n;
ttt mul(ttt &a,ttt &b){ //矩阵A*B
    ttt c;
    int i,j,k;
    for(i=1;i<=n;i++){
        for(j=1;j<=n;j++){
            c.map1[i][j]=0;
            for(k=1;k<=n;k++){
                c.map1[i][j]+=a.map1[i][k]*b.map1[k][j];
                c.map1[i][j]%=MOD;
            }
        }
    }
    return c;
}
ttt mul_self(ttt &a){
    ttt b;
    int i,j,k;
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++){
           b.map1[i][j]=0;
           for(k=1;k<=n;k++){
                b.map1[i][j]+=a.map1[i][k]*a.map1[k][j];
                b.map1[i][j]%=MOD;
           }
        }
        return b;
}
ttt pow1(ttt &a,int k){
    ttt b;
    int i,j;
    memset(b.map1,0,sizeof(b.map1));
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            if(i==j)
                b.map1[i][j]=1;
            while(k){
                //cout << "k=" <<k <<endl;
                if(k&1)
                    b=mul(b,a);
                    k=k>>1;
                    a=mul_self(a);
            }
            return b;
}

int main(){
    freopen("in.txt","r",stdin);
    int i,j,k,l,f1,f2,f3,t1,t2,t3;
    int r,c,m;
    int a1,b1,c1;
    int T;
    ttt a;
    memset(a.map1,0,sizeof(a.map1));
    cin >> n >> m;//n个结点，m条边
    for(i=1;i<=m;i++){
        cin >>t1>> t2;
        a.map1[t1][t2]=1;
    }
    a=pow1(a,2);
    int sum1=0;
    for(i=1;i<=n;i++)
    for(j=1;j<=n;j++){
        sum1+=a.map1[i][j];
    }
    c