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;
}