1305 Freda的道路

 时间限制: 1 s
空间限制: 128000 KB
题目等级: 大师 Master

题目描写叙述 Description

Freda要到Rainbow的城堡去玩了。

我们能够觉得两座城堡位于同一条数轴上。Freda的城堡坐标是0,Rainbow的城堡坐标是N。正常情况下,Freda会朝着同一个方向(即Rainbow的城堡相对于Freda的城堡的方向)走若干步之后来到Rainbow的城堡。并且步长都为1或2。

但是。今天Freda在途中遇见了来自上海的小猫Resodo,惊奇之下,竟然有一步走反了方向!只是。Freda并没有神智不清。它仅仅有一步走反了方向, 并且这一步的步长也是1或2. 同一时候,Freda并不会路过Rainbow的城堡而不停下来。

当然,Freda是在途中遇到Resodo的。所以它不会在 自己家门口就走错方向。

举个例子。假设Rainbow的城堡坐标是3,那么以下两条路径是合法的:

0->1->2->1->3
0->1->-1->1->3

当然,还有其他的合法路径。以下这些路径则是不合法的:

0->-1->1->3 (Freda不可能第一步就走错方向)
0->1->3(Freda一定是有一步走错方向的)
0->2->1->0->2->3(Freda仅仅有一步是走错方向的)
0->-1->0->3(Freda每步的长度一定是1或2)
0->1->2->4->3(Freda不会越过Rainbow的城堡再回来)
0 -> 1 -> 2 -> 3 -> 2 -> 3(Freda一旦到达了Rainbow的城堡,就会停下来)

你如今须要帮助Freda求出。它一共同拥有多少种方法能够到达Rainbow的城堡呢?

输入描写叙述 Input Description

一行一个整数N,表示Rainbow城堡的坐标

输出描写叙述 Output Description

一行一个整数,表示Freda到Rainbow城堡的不同路径数。
由于这个数字可能非常大。你仅仅须要输出它mod 1000000007的结果。

例子输入 Sample Input

2

例子输出 Sample Output

5

数据范围及提示 Data Size & Hint

对于第一组例子,例如以下5条路径是合法的:

0->1->0->2
0->1->-1->0->1->2
0->1->-1->0->2
0->1->0->1->2
0->1->-1->1->2

数据范围与约定

对于10%的数据,N<=20.
对于70%的数据,N<=1000.
对于90%的数据,N<=1000000.
对于100%的数据。N<=10^15.

题解:

由于是讨论路径个数。所以首先想到的是dp。

设f[i]表示在不回头的情况下到达的路径数,g[i]表示在一次回头的情况下到达的路径数。

由题意可知:

f[i]=f[i-2]+f[i-1](当前仅仅能从两步前或是一步前转移来),即Fibonacci数列。

接下来考虑g[i]。发现g[i]由四种情况转移过来,即:

1>g[i-2](在转过的情况下。由两步前走来)

2>g[i-1](在转过的情况下。由一步前走来)

3>f[i+1](一步以后转回来)

4>f[i+2](两步以后转回来)

于是可得转移方程:

g[i]=g[i-2]+g[i-1]+f[i+1]+f[i+2]

裸的dp可得90 。

。。

最后一个点是极限数据,DP做不到。于是想到了矩阵乘。

对于以下两个矩阵:

A⎡⎣⎢⎢⎢0,1,0,01,1,1,10,0,0,10,0,1,1⎤⎦⎥⎥⎥
B⎡⎣⎢⎢⎢⎢g[n−2]g[n−1]f[n+1]f[n+2]⎤⎦⎥⎥⎥⎥

g[n]能够这样求得:

A∗B=C⎡⎣⎢⎢⎢⎢g[n−1]g[n]f[n+2]f[n+3]⎤⎦⎥⎥⎥⎥

于是我们能够这样求解g[n]:

An−2∗⎡⎣⎢⎢⎢⎢g[1]g[2]f[4]f[5]⎤⎦⎥⎥⎥⎥

Code:

DP:

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define p 1000000007LL
using namespace std;
long long n,f[1000010],g[1000010];
int main(){
scanf("%d",&n);
if (n==1){
printf("0\n");
return 0;
}
f[1]=f[2]=g[0]=g[1]=1;
for (int i=3; i<=n+1; i++)
f[i]=(f[i-1]+f[i-2])%p;
for (int i=2; i<=n+1; i++)
g[i]=(g[i-2]+g[i-1]+f[i+1]+f[i+2])%p;
printf("%d\n",g[n+1]);
return 0;
}

矩阵乘:

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define p 1000000007LL
using namespace std;
long long A[4][4]={{0,1,0,0},{1,1,1,1},{0,0,0,1},{0,0,1,1}};
long long B[4][4],C[4][4];
int main(){
long long n; bool f=true;
scanf("%lld",&n);
if (n<=2){
if (n==1) printf("0\n");
else printf("5\n");
return 0;
}
n-=2;
while (n){
if (n&1){
if (f){
for (int i=0; i<4; i++)
for (int j=0; j<4; j++)
B[i][j]=A[i][j];
f=false;
}
else {
memset(C,0,sizeof(C));
for (int i=0; i<4; i++)
for (int j=0; j<4; j++)
for (int k=0; k<4; k++)
C[i][j]=(C[i][j]+B[i][k]*A[k][j])%p;
for (int i=0; i<4; i++)
for (int j=0; j<4; j++)
B[i][j]=C[i][j];
}
}
n>>=1;
memset(C,0,sizeof(C));
for (int i=0; i<4; i++)
for (int j=0; j<4; j++)
for (int k=0; k<4; k++)
C[i][j]=(C[i][j]+A[i][k]*A[k][j])%p;
for (int i=0; i<4; i++)
for (int j=0; j<4; j++)
A[i][j]=C[i][j];
}
long long ans=(B[1][1]*5%p+B[1][2]*5%p+B[1][3]*8%p)%p;
printf("%lld\n",ans);
return 0;
}

【矩阵乘】【DP】【codevs 1305】Freda的道路的更多相关文章

  1. Codevs 1305 Freda的道路(矩阵乘法 DP优化)

    1305 Freda的道路 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 大师 Master 题目描述 Description Freda要到Rainbow的城堡去玩了.我们可以认 ...

  2. 矩阵优化dp

    链接:https://www.luogu.org/problemnew/show/P1939 题解: 矩阵优化dp模板题 搞清楚矩阵是怎么乘的构造一下矩阵就很简单了 代码: #include < ...

  3. bzoj 3120 矩阵优化DP

    我的第一道需要程序建矩阵的矩阵优化DP. 题目可以将不同的p分开处理. 对于p==0 || p==1 直接是0或1 对于p>1,就要DP了.这里以p==3为例: 设dp[i][s1][s2][r ...

  4. HDU - 2294: Pendant(矩阵优化DP&前缀和)

    On Saint Valentine's Day, Alex imagined to present a special pendant to his girl friend made by K ki ...

  5. [六省联考2017]组合数问题 (矩阵优化$dp$)

    题目链接 Solution 矩阵优化 \(dp\). 题中给出的式子的意思就是: 求 nk 个物品中选出 mod k 为 r 的个数的物品的方案数. 考虑朴素 \(dp\) ,定义状态 \(f[i][ ...

  6. 洛谷P3502 [POI2010]CHO-Hamsters感想及题解(图论+字符串+矩阵加速$dp\&Floyd$)

    洛谷P3502 [POI2010]CHO-Hamsters感想及题解(图论+字符串+矩阵加速\(dp\&Floyd\)) 标签:题解 阅读体验:https://zybuluo.com/Junl ...

  7. [BZOJ 4818/LuoguP3702][SDOI2017] 序列计数 (矩阵加速DP)

    题面: 传送门:https://www.lydsy.com/JudgeOnline/problem.php?id=4818 Solution 看到这道题,我们不妨先考虑一下20分怎么搞 想到暴力,本蒟 ...

  8. 树形DP codevs 1814 最长链

    codevs 1814 最长链  时间限制: 1 s  空间限制: 256000 KB  题目等级 : 钻石 Diamond 题目描述 Description 现给出一棵N个结点二叉树,问这棵二叉树中 ...

  9. Codevs 1183 泥泞的道路

    1183 泥泞的道路 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 传送门 题目描述 Description CS有n个小区,并且任意小区之间都有两条单向道路 ...

随机推荐

  1. spring cloud 学习(10) - 利用springfox集成swagger

    对绝大多数程序员而言,写接口文档是一件痛苦的事情,相对文档,他们更愿意写代码.最理想的情况就是:代码即文档!服务开发完成后,部署上去文档就自动生成,没错,这就是springfox + swagger要 ...

  2. 利用webBrowser获取页面iframe中的内容

    1.获取frame的document HtmlDocument htmlDoc = webBrowser1.Document;  htmlDoc = webBrowser1.Document.Wind ...

  3. C#消息队列(RabbitMQ)零基础从入门到实战演练

    一.课程介绍 如果您从工作中之听过但未有接触过消息对队列(MQ),如果你接触过一点关于MQ的知识,如果没有这么的多如果的话......,那么阿笨将通过本次<C#消息队列零基础从入门到实战演练&g ...

  4. AngularJS自定义Directive不一定返回对象

    AngularJS中,当需要自定义Directive时,通常返回一个对象,就像如下的写法: angular.module('modulename') .directive('myDirective', ...

  5. 使用RemObjects Pascal Script

    摘自RemObjects Wiki 本文提供RemObjects Pascal Script的整体概要并演示如何创建一些简单的脚本. Pascal Script包括两个不同部分: 编译器 (uPSCo ...

  6. IIS7下配置ASP+ACCESS环境

    先要设置应用程序池(Application Pool)为Classic .NET AppPool,而不是默认的Default AppPool,可以在网站目录里对每个站点设置,也可以在站点进行单独设置. ...

  7. Linux进程ID号--Linux进程的管理与调度(三)

    转自:http://blog.csdn.net/gatieme/article/category/6225543 日期 内核版本 架构 作者 GitHub CSDN 2016-05-12 Linux- ...

  8. C#泛型委托Predicate、Action、Func

    Predicate Predicate泛型委托:表示定义一组条件并确定指定对象是否符合这些条件的方法.此委托由 Array 和 List 类的几种方法使用,用于在集合中搜索元素.通过查看源码发现 Pr ...

  9. java根据jar包反编译后修改再打包回jar的做法

    1. 得到一个待要修改的jar包 2. 我的环境是windows,然后解压这个jar包,得到一堆class文件,这时候就找到你需要的那个class文件 3. 我首先是使用jd-gui工具看一下这个cl ...

  10. hex string 换转

    hex2string std::stringstream R; R << std::hex << 0x1254; DWORD Str2Hex(string str){ retu ...