bzoj1297 [SCOI2009]迷路——拆点+矩阵快速幂

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1297

一看感觉是矩阵快速幂之类的,但边权不好处理啊;

普通的矩阵快速幂只能处理边权为1的,所以想办法把边权处理成1;

仔细一看还有一个条件是边权小于10;

所以拆点!把一个点拆成10个点表示到它不同的距离,那么和它相连的那些点就可以跟某个距离的点连边权为1的边;

虽然没有自己想出来,不过1A还是极好的!(因为太简单了)

代码如下:

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,T,w[15][15],tot,id[15][15],mod=2009;
struct Matrix{
    int a[105][105];
    Matrix(){memset(a,0,sizeof a);}
    Matrix operator * (const Matrix &y) const
    {
        Matrix x;
        for(int i=1;i<=tot;i++)
            for(int k=1;k<=tot;k++)
                for(int j=1;j<=tot;j++)
                    (x.a[i][j]+=a[i][k]*y.a[k][j])%=mod;
        return x;
    }
    void init(){for(int i=1;i<=tot;i++)a[i][i]=1;}
}f,ans;
Matrix pw(Matrix x,int k)
{
    Matrix ret; ret.init();
    for(;k;k>>=1,x=x*x)
        if(k&1)ret=ret*x;
    return ret;
}
int main()
{
    scanf("%d%d",&n,&T);
    char ch[15];
    for(int i=0;i<n;i++)
    {
        scanf("%s",&ch);
        for(int j=0;j<n;j++)
            w[i][j]=ch[j]-‘0‘;
    }
    for(int i=0;i<n;i++)
        for(int j=0;j<=9;j++)
        {
            id[i][j]=++tot;
            if(j)f.a[tot-1][tot]++;
        }
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            if(w[i][j])f.a[id[i][w[i][j]-1]][id[j][0]]++;
    ans.a[1][id[0][0]]=1;
    Matrix fn=pw(f,T);
    ans=ans*fn;
    printf("%d",ans.a[1][id[n-1][0]]);
    return 0;
}

原文地址:https://www.cnblogs.com/Zinn/p/9254628.html

时间: 2024-10-17 09:48:56

bzoj1297 [SCOI2009]迷路——拆点+矩阵快速幂的相关文章

BZOJ 1297: [SCOI2009]迷路 [矩阵快速幂]

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 包

BZOJ 1297 迷路(矩阵快速幂)

很容易想到记忆化搜索的算法. 令dp[n][T]为到达n点时时间为T的路径条数.则dp[n][T]=sigma(dp[i][T-G[i][n]]); 但是空间复杂度为O(n*T),时间复杂度O(n*n*T). 虽然本题的n<=10,但T最大可到1e9.行不通. 如果题目中的边的权值非0即1的话,显然1-n的长度为T的路径中数为 该图的邻接矩阵的T次幂. 实际上题目中的边权值<10. 可以用拆点的方法转化为边权值非0即1的情况. 即 将图中的每个点拆成至多9个点,首先将每个点的第i个点和第i+1

BZOJ1297 [SCOI2009]迷路

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 包

hdu 2842(矩阵快速幂+递推)

题意:一个中国环的游戏,规则是一个木棒上有n个环,第一个环是可以随意放上或拆下的,剩下的环x如果想放上或拆下必须前一个环x-1是放上的且前x-2个环全部是拆下的,问n个环最少多少次操作可以全部拆掉. 题解:需要进行递推,首先第一步肯定是要拆第n个环保证操作次数最少,因为后面的环是否存在对前面的环不造成影响,而先拆前面的如果要拆后面的环还是要把前面的放上,f(n)表示拆掉前n个环需要的最少操作次数,先拆第n个要拆前n-2个再拆第n个,花费f(n-2)+1,然后这时是00-0010,要拆第n-1个需

[BZOJ 2326][HNOI 2011]数学作业(矩阵快速幂)

蒟蒻线性代数太烂了...这个逼题居然卡了半天才做出来,弱的不行啊... 矩阵快速幂,把n这个len位数拆成len次分段快速幂就可以了. 注意取模的数字m<=1e9,所以矩阵乘法运算时要先对乘数取模,防止中间运算结果太大溢出,坑爹啊 代码: #include <iostream> #include <stdio.h> #include <stdlib.h> #include <string.h> #include <algorithm> #d

hdu 4878 ZCC loves words(AC自动机+dp+矩阵快速幂+中国剩余定理)

hdu 4878 ZCC loves words(AC自动机+dp+矩阵快速幂+中国剩余定理) 题意:给出若干个模式串,总长度不超过40,对于某一个字符串,它有一个价值,对于这个价值的计算方法是这样的,设初始价值为V=1,假如这个串能匹配第k个模式串,则V=V*prime[k]*(i+len[k]),其中prime[k]表示第k个素数,i表示匹配的结束位置,len[k]表示第k个模式串的长度(注意,一个字符串可以多次匹配同意个模式串).问字符集为'A'-'Z'的字符,组成的所有的长为L的字符串,

ACdream 1093 女神的正多面体 矩阵快速幂

题目大意:给你三种正多边形,给你起点s,终点e以及最多行走的步数k,问有多少种路径方案(路径中只要有一个顶点不同即视为不同). 题目分析: 可以通过矩阵快速幂求解. 为每个正多边形(最多三个)构建一个邻接矩阵A,然后第K步的方案数即为A^k. 结果即为A^1 + A^2 + A^3 + ...... + A^k. 对于这种形式的矩阵运算,我们可以把它拆分成: k为奇:ans = (A^1 + A^2 + ... + A^(k/2)) + (A^1 + A^2 + ... + A^(k/2)) *

二进制快速幂及矩阵快速幂

二进制快速幂 二进制快速幂虽然不难写,但是无奈总是会忘,所以还是在这里把板子写一下. 二进制快速幂很好理解: 假设我们要求a^b,那么其实b是可以拆成二进制的,该二进制数第i位的权为2^(i-1),例如当b==11时,11的二进制是1011,11 = 23×1 + 22×0 + 21×1 + 2o×1,因此,我们将a11转化为算 a2^0*a2^1*a2^3 int poww(int a, int b) { int ans = 1, base = a; while (b != 0) { if (

bzoj1875 [SDOI2009]HH去散步 矩阵快速幂

题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=1875 题解 如果没有这个"不能立刻沿着刚刚走来的路走回",那么这个题就是一个常规的矩阵乘法. 考虑一下这个限制怎么解决.因为限制的是边,我们不妨考虑和边有关的矩阵. 首先把一条无向边拆成两个有向边,如果边 \(A\) 的终点和边 \(B\) 的起点相同,那么我们就说从边 \(A\) 通向边 \(B\).但是,同源的有向边(也就是从同一条无向边拆成的两条有向边)之间不能建边. 但是为