刚开始把建立一个矩阵,用以表示任意两个顶点之间是否有边,如果有矩阵上就为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