Description
因为SJY干的奇怪事情过多,SJY收到了休假的通知,于是他准备在都市间来回旅游。SJY有一辆车子,一开始行驶性能为0,每过1时间行驶性能就会提升1点。每个城市的道路都有性能要求。SJY一共有t时间休息,一开始他位于1号城市(保证1号城市道路要求为0),他希望在n号城市结束旅程。每次穿过一条城市间的路会花费1时间,当然他也可以停留在一个城市不动而花费1时间。当且仅当车子的行驶性能大于等于一个城市,我们才能到达那里。SJY希望知道,旅游的方案模10086后的答案。(只要在某一时刻通过的道路存在一条不相同,就算不同的方案)
Input
第一行三个数n,m,t,表示有n个城市m条道路t时间。
第二行n个数,hi表示第i个城市的道路性能要求。
第三到m+2行,每行两个数u,v,表示城市u与城市v之间有一条单向道路连接(可能有重边)。
Output
包括一个数字,表示旅游的方案模10086。
Sample Input
5 17 7
0 2 4 5 3
1 2
2 1
1 3
3 1
1 4
4 1
4 5
5 4
5 3
4 1
2 1
5 3
2 1
2 1
1 2
2 1
1 3
Sample Output
245
Data Constraint
对于20%的数据,n<=10,t<=80;
对于50%的数据,n<=30,t<=80;
对于100%的数据,n<=70,m<=1000,t<=100000000,hi<=70。
题解
观察数据范围,n很小,限制的时间也很小,但总的时间非常大,由此想到矩阵乘法。
对于前面有限制的就直接暴力计算,
到了后面就已经没有限制的了,
每次的转移矩阵都是一样的,
即连边的矩阵。
矩阵乘法快速幂,就得到了答案。
code
#include<queue>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include <cstring>
#include <string.h>
#include <cmath>
#include <math.h>
#include <time.h>
#define ll long long
#define N 73
#define M 100003
#define db double
#define P putchar
#define G getchar
#define inf 998244353
#define mo 10086
using namespace std;
char ch;
void read(int &n)
{
n=0;
ch=G();
while((ch<'0' || ch>'9') && ch!='-')ch=G();
ll w=1;
if(ch=='-')w=-1,ch=G();
while('0'<=ch && ch<='9')n=(n<<3)+(n<<1)+ch-'0',ch=G();
n*=w;
}
int max(int a,int b){return a>b?a:b;}
int min(int a,int b){return a<b?a:b;}
ll abs(ll x){return x<0?-x:x;}
ll sqr(ll x){return x*x;}
void write(ll x){if(x>9) write(x/10);P(x%10+'0');}
struct node
{
int s[N][N];
}p;
int nxt[M],to[M],lst[N],tot;
int f[N][N],h[N],n,m,k,ans,mx,x,y;
node operator *(node a,node b)
{
node s;
memset(s.s,0,sizeof(s.s));
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++)
s.s[i][j]=(s.s[i][j]+a.s[i][k]*b.s[k][j])%mo;
return s;
}
node ksm(node x,int y)
{
node s;
memset(s.s,0,sizeof(s.s));
for(int i=1;i<=n;i++)
s.s[i][i]=1;
while(y)
{
if(y&1)s=s*x;
x=x*x;
y>>=1;
}
return s;
}
void ins(int x,int y)
{
nxt[++tot]=lst[x];
to[tot]=y;
lst[x]=tot;
}
int main()
{
read(n);read(m);read(k);
for(int i=1;i<=n;i++)
read(h[i]),mx=max(mx,h[i]),p.s[i][i]++;
for(int i=1;i<=m;i++)
read(x),read(y),p.s[x][y]++,ins(x,y);
f[0][1]=1;mx=min(mx,k);
for(int i=1;i<=mx;i++)
{
for(int j=1;j<=n;j++)
f[i][j]=f[i-1][j];
for(int j=1;j<=n;j++)
if(f[i-1][j]>0)
{
for(int t=lst[j];t;t=nxt[t])
if(h[to[t]]<=i)f[i][to[t]]=(f[i][to[t]]+f[i-1][j])%mo;
}
}
p=ksm(p,k-mx);
for(int i=1;i<=n;i++)
ans=(ans+f[mx][i]*p.s[i][n])%mo;
write(ans);
return 0;
}