[Bzoj1297][Scoi2009 ]迷路 （矩阵乘法 + 拆点）

1297: [SCOI2009]迷路

---

 

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1385  Solved: 993
[Submit][Status][Discuss]

Description

---

 

windy在有向图中迷路了。 该有向图有 N 个节点，windy从节点 0 出发，他必须恰好在 T 时刻到达节点 N-1。 现在给出该有向图，你能告诉windy总共有多少种不同的路径吗？ 注意：windy不能在某个节点逗留，且通过某有向边的时间严格为给定的时间。

Input

---

 

第一行包含两个整数，N T。 接下来有 N 行，每行一个长度为 N 的字符串。 第i行第j列为'0'表示从节点i到节点j没有边。 为'1'到'9'表示从节点i到节点j需要耗费的时间。

Output

---

 

包含一个整数，可能的路径数，这个数可能很大，只需输出这个数除以2009的余数。

Sample Input

---

 

【输入样例一】

2 2
11
00

 

---

【输入样例二】

5 30
12045
07105
47805
12024
12345

 

Sample Output

---

 

【输出样例一】

852

 

---

【样例解释一】

0->0->1

 

---

【输出样例二】

852

 

---

HINT

 

---

 

30%的数据，满足 2 <= N <= 5 ； 1 <= T <= 30 。 100%的数据，满足 2 <= N <= 10 ； 1 <= T <= 1000000000 。

 

分析：

 

---

 

  先用dp的思想分析可以分析出转移方程dp[i][t] = ∑dp[j][t  -  time（i，j）]；但是t太大，并且时间也不允许；

  想了一下发现转移方程是固定的，可以很自然联想矩阵乘法 + 快速幂。

  但是写的时候出现困难了，发现每个点对间有些彼此到达时间各不相同不好处理。于是就困住了

 

正解：

---

 

 

可以把一个点拆出9个虚拟结点出来，node[i][j] 向 node[i][j + 1]点向虚拟结点连边为1，然后所有到达这个点的都连边向node[i][0],所有以这个点出发能到距离为j的点q从node[i][j - 1]连边1向node[q][0]；

于是原图换成了一个10 * 10的矩阵，矩乘 + 快速幂求出即行。

复杂度还是很低的。

View Code

 

 

 

 

 

 

矩阵乘法：

---

 

其实关于dp的东西如果发现是固定转移并且是刷表什么的，都可以联想到转换成矩乘来做。

AC代码：

---

 

 

# include <iostream>
# include <cstdio>
# include <cstring>
using namespace std;
const int mod = 2009;
int n;
long long t;
struct fi{
    int data[102][102];
    fi(){memset(data,0,sizeof(data));
    }
}A,T;
inline fi operator * (fi & a,fi & b){
    fi t;
    for(int i = 0;i < n * 9;i++){
        for(int j = 0;j < n * 9;j++){
            t.data[i][j] = 0;
            for(int k = 0;k < n * 9;k++){
                t.data[i][j] += a.data[i][k] * b.data[k][j]; 
            }
            if(t.data[i][j] >= mod)t.data[i][j] %= mod;
        }
    }
    return t;
} 
inline fi operator +(fi & a,fi & b){
    fi t;
    for(int i = 0;i < n * 9;i++){
        for(int j = 0;j < n * 9;j++){
            t.data[i][j] = a.data[i][j] + b.data[i][j];
            if(t.data[i][j] >= mod)t.data[i][j] %= mod;
        }
    }
    return t;
}
inline void cmd(fi & A,fi & T,long long k){
    while(k){
        if(k & 1LL)A = A * T;
        T = T * T;
        k >>= 1LL;
    }
    return;
}
char str[102];
void build(){
    scanf("%d %lld",&n,&t);
    int x;
    for(int i = 0;i < n;i++){
        scanf("%s",str);
      for(int j = 0;j < n;j++){
         x = str[j] - '0';
         if(!x)continue;
         x -= 1;
         T.data[i * 9 + x][j * 9] = 1;
      }
      for(int j = 0;j < 8;j++){
        T.data[i * 9 + j][i * 9 + j + 1] = 1;
      }
    }
    A.data[0][0] = 1;   
     cmd(A,T,t);
    printf("%d\n",A.data[0][(n - 1) * 9]);
}
int main(){
    build();
}

 

转载于:https://www.cnblogs.com/lzdhydzzh/p/7652584.html