JZOJ 2679. 沼泽鳄鱼 【矩阵乘法】【7.6 A组 T2】【ZJOI2005】

题目大意：

潘塔纳尔沼泽地号称世界上最大的一块湿地，它地位于巴西中部马托格罗索州的南部地区。每当雨季来临，这里碧波荡漾、生机盎然，引来不少游客。
　　为了让游玩更有情趣，人们在池塘的中央建设了几座石墩和石桥，每座石桥连接着两座石墩，且每两座石墩之间至多只有一座石桥。这个景点造好之后一直没敢对外开放，原因是池塘里有不少危险的食人鱼。
　　豆豆先生酷爱冒险，他一听说这个消息，立马赶到了池塘，想做第一个在桥上旅游的人。虽说豆豆爱冒险，但也不敢拿自己的性命开玩笑，于是他开始了仔细的实地勘察，并得到了一些惊人的结论：食人鱼的行进路线有周期性，这个周期只可能是2，3或者4个单位时间。每个单位时间里，食人鱼可以从一个石墩游到另一个石墩。每到一个石墩，如果上面有人它就会实施攻击，否则继续它的周期运动。如果没有到石墩，它是不会攻击人的。
　　借助先进的仪器，豆豆很快就摸清了所有食人鱼的运动规律，他要开始设计自己的行动路线了。每个单位时间里，他只可以沿着石桥从一个石墩走到另一个石墩，而不可以停在某座石墩上不动，因为站着不动还会有其它危险。如果豆豆和某条食人鱼在同一时刻到达了某座石墩，就会遭到食人鱼的袭击，他当然不希望发生这样的事情。
　　现在豆豆已经选好了两座石墩Start和End，他想从Start出发，经过K个单位时间后恰好站在石墩End上。假设石墩可以重复经过（包括Start和End），他想请你帮忙算算，这样的路线共有多少种（当然不能遭到食人鱼的攻击）。

解题思路：

鱼的周期只为2，3，4，所以最多12步一个周期
我们可以把每一时刻的邻接矩阵相乘，便可知道路径的数目
当a[i,k]>0,b[k,j]>0时，即i—>k和k—>j都有路径时c[I,j]才有路径，且等于实际i—j的路径
然后矩乘随便搞搞

$Acceptedcode:$

#include<cstdio>
#include<cstring>

using namespace std;

typedef long long ll;

const int N = 60;
const ll P = 1e4;

bool use[N];
int n, fish, m, st, ed, K;
int fp[N][N];

struct Rect {
    ll a[N][N];
    ll *operator[](int i) { return a[i]; }
    void In() {
        for (int i = 0; i < n; ++i) a[i][i] = 1;
    }
    void Empty() { memset(a, 0, sizeof a); }
    Rect operator*(Rect b) {
        Rect c;
        c.Empty();
        for (int i = 0; i < n; ++i)
            for (int j = 0; j < n; ++j)
                for (int k = 0; k < n; k++)
                    if (!use[k])
                        (c[i][j] += a[i][k] * b[k][j] % P) %= P;
        return c;
    }
}f, ans, e;

Rect Ksm(Rect a, int y) {
    Rect ans;
    memset(use, 0, sizeof use);
    ans.In();
    while (y) {
        if (y & 1)
            ans = ans * a;
        a = a * a;
        y >>= 1;
    }
    return ans;
}

void Pre_Process(int x) {
    memset(use, 0, sizeof use);
    for (int i = 1; i <= fish; ++i) use[fp[i][(x-1)%fp[i][0]+1]] = 1;
}

int main() {
    scanf("%d %d %d %d %d", &n, &m, &st, &ed, &K);
    e.In(), ans.In();
    for (int i = 1, u, v; i <= m; ++i)
		scanf("%d%d", &u, &v), ++f[u][v], ++f[v][u];
    scanf("%d", &fish);
    for (int i = 1; i <= fish; ++i) {
        scanf("%d", &fp[i][0]);
        for (int j = 1; j <= fp[i][0]; ++j) scanf("%d", &fp[i][j]);
    }
    for (int i = 1; i <= 12; ++i) Pre_Process(i), e = e * f;
    ans = ans * Ksm(e, K / 12);
    for (int i = 1; i <= K % 12; ++i) Pre_Process(i), ans = ans * f;
    printf("%lld", ans.a[st][ed]);
}