[Jzoj] 3847. 都市环游

题目描述

因为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];
}