CF954F Runner's Problem(动态规划,矩阵快速幂)

题面

CodeForces

翻译:

有一个\(3\times M\)的田野

一开始你在\((1,2)\)位置

如果你在\((i,j)\)位置

在不出界的前提下,可以走到\((i+1,j),(i+1,j±1)\)

有\(n\)段障碍,障碍不能走

询问从\((1,2)\)到达\((M,2)\)的方案数

\(n<=10^4,M<=10^{18}\)

题解

发现\(M\)的范围非常大

很容易往矩阵快速幂的方向考虑

如果知道上一行的方案,以及这一行的状态

很容易可以列出转移矩阵

所以,将所有的障碍段离散

将\(3\)行分离,检查当前段的障碍组成

构建出转移矩阵

分段做快速幂就行了

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 11111
#define MOD 1000000007
inline ll read()
{
RG ll x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
struct Maxtrix
{
ll s[5][5];
void clear(){memset(s,0,sizeof(s));}
void init(){clear();s[1][1]=s[2][2]=s[3][3]=1;}
void M000(){clear();s[1][1]=s[1][2]=s[2][1]=s[2][2]=s[2][3]=s[3][2]=s[3][3]=1;}
};
Maxtrix operator*(Maxtrix a,Maxtrix b)
{
Maxtrix ret;ret.clear();
for(int i=1;i<=3;++i)
for(int j=1;j<=3;++j)
for(int k=1;k<=3;++k)
ret.s[i][j]=(ret.s[i][j]+1ll*a.s[i][k]*b.s[k][j]%MOD)%MOD;
return ret;
}
Maxtrix fpow(Maxtrix a,ll b)
{
Maxtrix s;s.init();
while(b){if(b&1)s=s*a;a=a*a;b>>=1ll;}
return s;
}
struct Block{ll l,r;int a;}blk[MAX];
ll tot,top;
ll S[MAX<<2];
ll c[4][MAX<<1];
ll n,Q;
ll M;
int main()
{
n=read();M=read();
for(int i=1;i<=n;++i)
{
int a=read();ll l=read(),r=read();
S[++top]=l-1,S[++top]=r;
blk[++tot]=(Block){l,r,a};
}
S[++top]=1;S[++top]=M;
sort(&S[1],&S[top+1]);
top=unique(&S[1],&S[top+1])-S-1;
for(int i=1;i<=n;++i)
{
ll L=lower_bound(&S[1],&S[top+1],blk[i].l)-S;
ll R=lower_bound(&S[1],&S[top+1],blk[i].r)-S;
c[blk[i].a][L]++;c[blk[i].a][R+1]--;
}
Maxtrix ans;ans.clear();ans.s[2][1]=1;
ll ss[4];ss[1]=ss[2]=ss[3]=0;
for(int i=2;i<=top;++i)
{
ll len=S[i]-S[i-1];
Maxtrix now;now.M000();
for(int j=1;j<=3;++j)
{
ss[j]+=c[j][i];
if(ss[j])now.s[j][1]=now.s[j][2]=now.s[j][3]=0;
}
now=fpow(now,len);
ans=now*ans;
}
cout<<ans.s[2][1]<<endl;
return 0;
}

CF954F Runner's Problem(动态规划,矩阵快速幂)的更多相关文章

  1. HDU1757-A Simple Math Problem,矩阵快速幂,构造矩阵水过

    A Simple Math Problem 一个矩阵快速幂水题,关键在于如何构造矩阵.做过一些很裸的矩阵快速幂,比如斐波那契的变形,这个题就类似那种构造.比赛的时候手残把矩阵相乘的一个j写成了i,调试 ...

  2. poj 3744 Scout (Another) YYF I - 概率与期望 - 动态规划 - 矩阵快速幂

      (Another) YYF is a couragous scout. Now he is on a dangerous mission which is to penetrate into th ...

  3. A Simple Math Problem(矩阵快速幂)(寒假闭关第一题,有点曲折啊)

    A Simple Math Problem Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Other ...

  4. Educational Codeforces Round 60 (Rated for Div. 2) - D. Magic Gems(动态规划+矩阵快速幂)

    Problem   Educational Codeforces Round 60 (Rated for Div. 2) - D. Magic Gems Time Limit: 3000 mSec P ...

  5. dutacm.club Water Problem(矩阵快速幂)

    Water Problem Time Limit:3000/1000 MS (Java/Others)   Memory Limit:163840/131072 KB (Java/Others)Tot ...

  6. LightOJ 1070 Algebraic Problem:矩阵快速幂 + 数学推导

    题目链接:http://lightoj.com/volume_showproblem.php?problem=1070 题意: 给你a+b和ab的值,给定一个n,让你求a^n + b^n的值(MOD ...

  7. A Simple Math Problem (矩阵快速幂)

    Lele now is thinking about a simple function f(x).  If x < 10 f(x) = x.  If x >= 10 f(x) = a0 ...

  8. HDU——4291A Short problem(矩阵快速幂+循环节)

    A Short problem Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)T ...

  9. hdu 1757 A Simple Math Problem (矩阵快速幂)

    Description Lele now is thinking about a simple function f(x). If x < 10 f(x) = x. If x >= 10 ...

  10. hdu 2604 Queuing(动态规划—>矩阵快速幂,更通用的模版)

    题目 最早不会写,看了网上的分析,然后终于想明白了矩阵是怎么出来的了,哈哈哈哈. 因为边上的项目排列顺序不一样,所以写出来的矩阵形式也可能不一样,但是都是可以的 //愚钝的我不会写这题,然后百度了,照 ...

随机推荐

  1. python-全栈开发-前方高能-内置函数

    python_day_14 13. 前方高能-内置函数 ⼀. 本节主要内容: 1. 内置函数 什么是内置函数? 就是python给你提供的. 拿来直接⽤的函数, 比如print., input等等. ...

  2. Mac下布置appium环境

    1.下载或者更新Homebrew:homebrew官网 macOS 不可或缺的套件管理器 $ /usr/bin/ruby -e "$(curl -fsSL https://raw.githu ...

  3. 英特尔近日发布最新版实感™ SDK R5 (v7)

    英特尔开发人员专区 原文地址 英特尔® 实感™ SDK 的 7.0.23.8048 版本(也称为 R5)现已推出.您将看到的主要变化包括: 支持英特尔® 实感™ SR300 摄像头:应于 2016 年 ...

  4. LJ语录

    "保持安静,不要打扰我睡觉." ( 半分钟后) "哎呦喂~可以睡觉喽~" "考场上遇到这种题目,你们应该高兴." "因为考场上那群 ...

  5. 【Python入门学习】闭包&装饰器&开放封闭原则

    1. 介绍闭包 闭包:如果在一个内部函数里,对在外部作用域的变量(不是全局作用域)进行引用,那边内部函数被称为闭包(closure) 例如:如果在一个内部函数里:func2()就是内部函数, 对在外部 ...

  6. 5.airflow问题

    1. Traceback (most recent call last): File "/usr/bin/airflow", line 28, in <module> ...

  7. 面对30页左右的运放数据手册datasheet,你需要知道如何看懂

    1.输入失调电压(Input Offset Voltage) VOS    若将运放的两个输入端接地,理想运放输出为零,但实际运放输出不为零.此时,用输出电压除以增益得到的等效输入电压称为输入失调电压 ...

  8. Divide by three, multiply by two(DFS+思维)

    Polycarp likes to play with numbers. He takes some integer number x, writes it down on the board, an ...

  9. Sorting a Three-Valued Sequence(三值排序)

    Description 排序是一种很频繁的计算任务.现在考虑最多只有三值的排序问题.一个实际的例子是,当我们给某项竞赛的优胜者按金银铜牌序的时候. 在这个任务中可能的值只有三种1,2和3.我们用交换的 ...

  10. Alpha阶段中间产物——GUI Prototype、WBS及PSP

    作业地址:https://edu.cnblogs.com/campus/nenu/SWE2017FALL/homework/1224 内容: GUI Prototype 我的书架 我的书架→添加图书 ...