P2579 [ZJOI2005]沼泽鳄鱼(邻接矩阵,快速幂)

Posted ajmddzp

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了P2579 [ZJOI2005]沼泽鳄鱼(邻接矩阵,快速幂)相关的知识,希望对你有一定的参考价值。

题目简洁明了(一点都不好伐)

照例,化简题目

给一张图,每一个时间点有一些点不能走,(有周期性),求从起点第k秒恰好在终点的方案数,可重复,不可停留。

额dp实锤

于是就被打脸了....

有一种东西叫做邻接矩阵,还有一种东西叫做矩阵乘法,bk201大仙曾经讲过,能用邻接矩阵的k次方求这个东西。

那,难度下降了很多了。

但是,对于那周期是3,4,6的鳄鱼怎么办呢?

答案就是:

观察3,4,6,lcm(3,4,6)=12,所以,当k%12==0,鳄鱼们都回到了原状态。

所以,我们就暴力地对于每一种情况建一个邻接矩阵,把不能走的点刨掉,然后把他们依次乘起来(因为不满足交换律)

然后把它给k次方就行了。。。

值得注意的地方:

1、因为点是从0开始,所以根据我的习惯,都加上了1;

2、不满足交换律不满足交换律不满足交换律

3、k不一定整除12,所以在剩余的那一部分要手动乘一下初始矩阵

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=55;
const int mod=10000;
int m,n,st,en,K,fn;
int map[maxn][maxn];
int unable[maxn][maxn];
struct node

    int a[maxn][maxn];
    node()
    
        memset(a,0,sizeof(a));
    
q,w[55];
inline node operator * (const node& x,const node& y)

    node z;
    for(int i=1;i<=n;i++)
    
        for(int k=1;k<=n;k++)
        
            for(int j=1;j<=n;j++)
            
                z.a[i][j]=(z.a[i][j]+x.a[i][k]*y.a[k][j])%mod;
            
        
    
    return z;

node ksm(node x,int y)

    node z;
    for(int i=1;i<=n;i++)
    z.a[i][i]=1;
    while(y)
    
        if(y&1)
        z=z*x;
        x=x*x;
        y>>=1;
    
    return z;


int main()

    scanf("%d%d%d%d%d",&n,&m,&st,&en,&K);
    st++;en++;
    for(int i=1;i<=m;i++)
    
        int x,y;
        scanf("%d%d",&x,&y);
        x++;y++;
        map[x][y]=map[y][x]=1;
    
    scanf("%d",&fn);
    while(fn--)
    
        int x;
        scanf("%d",&x);
        for(int j=0;j<x;j++)
        
            int y;
            scanf("%d",&y);
            y++;
            for(int k=j;k<12;k+=x)
            
                unable[k][y]=true;
            
        
    

    for(int i=1;i<=n;i++)
    
        q.a[i][i]=1;
    
    for(int i=0;i<12;i++)
    
        for(int j=1;j<=n;j++)
        
            for(int k=1;k<=n;k++)
            
                w[i].a[j][k]=(map[j][k]&&!unable[i][k]);
            
        
    
    for(int i=1;i<12;i++)
    
        q=q*w[i];
    
    q=q*w[0];
    int k=K/12,t=K-k*12;
    node ans;
    for(int i=1;i<=n;i++)
    
        ans.a[i][i]=1;
    
    ans=ans*ksm(q,k);
    for(int i=1;i<=t;i++)
    ans=ans*w[i];
    printf("%d ",ans.a[st][en]);
    return 0;

(完)

以上是关于P2579 [ZJOI2005]沼泽鳄鱼(邻接矩阵,快速幂)的主要内容,如果未能解决你的问题,请参考以下文章

ZJOI2005沼泽鳄鱼

bzoj 1898 [Zjoi2005]Swamp 沼泽鳄鱼——矩阵快速幂

bzoj1898 [Zjoi2005]Swamp 沼泽鳄鱼——矩阵快速幂

BZOJ 1898: [Zjoi2005]Swamp 沼泽鳄鱼 [矩阵乘法]

BZOJ1898[Zjoi2005]Swamp 沼泽鳄鱼 矩阵乘法

[luogu2579 ZJOI2005] 沼泽鳄鱼(矩阵快速幂)