JZOJ3847. 【NOIP2014八校联考第1场第2试9.21】都市环游(travel)

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