hdu 5411
分析:
首先我們知道了鄰接矩陣A,那么A^k代表的就是長度為k的路徑有多少個。那么結果就是A^0+A^1+A^2+...+A^m。然后我們可以構造一個矩陣來幫助我們求解。
X = [A,E;0,E]
?==> 然后X^m 的矩陣的右上角的矩陣代表的就是A^0+A^1+A^2+...+A^m。
代碼:
#include<iostream> #include<cstdio> #include<cstring> using namespace std;const int maxn = 105; const int mod = 2015; struct Matrix {int m[maxn][maxn];Matrix(){memset(m,0,sizeof(m));} }; int n,m; Matrix I;Matrix mul(Matrix a,Matrix b) {Matrix c;for(int i = 1; i <= 2*n; i++)for(int j = 1; j <= 2*n; j++){for(int k = 1; k <= 2*n; k++){c.m[i][j] = (c.m[i][j] + a.m[i][k]*b.m[k][j]) % mod;}}return c; }Matrix power(int k,Matrix A) {Matrix ans = I;while(k){if(k & 1) ans = mul(ans,A);A = mul(A,A);k = k >> 1;}return ans; }void init() {for(int i = 1; i < maxn; i++)I.m[i][i] = 1; }int main() { int t;init();cin>>t;while(t--){cin>>n>>m;Matrix A;for(int i = 1; i <= n; i++){int k,a;cin>>k;while(k--){cin>>a;A.m[i][a] = 1;}}for(int i=1; i<=n; i++) A.m[i][i+n] = 1; for(int i=n+1; i<=2*n; i++) A.m[i][i] = 1; int sum = 0;A = power(m,A);for(int i = 1; i <= n; i++)for(int j = n+1; j <= 2*n; j++)sum = (sum + A.m[i][j]) % mod;cout<<sum+1<<endl;}return 0; }
總結
- 上一篇: MiniDao_1.6.4 版本发布,轻
- 下一篇: 费马素性测试和米勒—拉宾素性测试