题目描述
因为SJY干的奇怪事情过多,SJY收到了休假的通知,于是他准备在都市间来回旅游。SJY有一辆车子,一开始行驶性能为0,每过1时间行驶性能就会提升1点。每个城市的道路都有性能要求。SJY一共有t时间休息,一开始他位于1号城市(保证1号城市道路要求为0),他希望在n号城市结束旅程。每次穿过一条城市间的路会花费1时间,当然他也可以停留在一个城市不动而花费1时间。当且仅当车子的行驶性能大于等于一个城市,我们才能到达那里。SJY希望知道,旅游的方案模10086后的答案。(只要在某一时刻通过的道路存在一条不相同,就算不同的方案)
题目解析
将每一步移动的转移矩阵都求出来,然后进行矩阵乘法。并用矩阵快速幂优化矩阵乘法。
对于刚开始的限制条件,不断更新转移矩阵就可以了
代码
#include<bits/stdc++.h>
#define N 105
#define M 10086
#define ll long long
using namespace std;
int n,m,t,maxh;
int h[N],a[N][N],f[N];
int mat[N][N];
void mul(int a[],int b[][N])
{
int c[N];memset(c,0,sizeof(c));
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
(c[i]+=a[j]*b[j][i])%=M;
memcpy(f,c,sizeof(c));
}
void mulself()
{
int c[N][N];memset(c,0,sizeof(c));
for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)for(int k=1;k<=n;k++)c[i][j]=(c[i][j]+mat[k][j]*mat[i][k]%M)%M;
memcpy(mat,c,sizeof(c));
}
int main()
{
scanf("%d%d%d",&n,&m,&t);
for(int i=1;i<=n;i++) scanf("%d",&h[i]),maxh=max(maxh,h[i]);
for(int i=1,x,y;i<=m;i++) scanf("%d%d",&x,&y),a[x][y]++;
for(int i=1;i<=n;i++) a[i][i]=1;
f[1]=1;
for(int i=0;i<=min(maxh,t);i++)
{
for(int j=1;j<=n;j++)
if(h[j]==i)
for(int k=1;k<=n;k++)
mat[k][j]=a[k][j];
mul(f,mat);
}
if(maxh>=t)
{
cout<<f[n];
return 0;
}
int x=t-maxh;
while(x)
{
if(x&1) mul(f,mat);
x>>=1;
mulself();
}
cout<<f[n];
}