博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
hdu 2842 Chinese Rings 矩阵快速幂
阅读量:5767 次
发布时间:2019-06-18

本文共 1857 字,大约阅读时间需要 6 分钟。

分析:

后面的环能不能取下来与前面的环有关,前面的环不被后面的环所影响。所以先取最后面的环

设状态F(n)表示n个环全部取下来的最少步数

先取第n个环,就得使1~n-2个环属于被取下来的状态,第n-1个环属于未被取下来的状态。那么F(n) = F(n-2) + 1 + ... (这里的1表示取下第n个需要一步)

即F(n)可以为F(n-2) + 1与某些数的和。取下n后,1~n-2为取下的状态,n-1为未被取下的状态。如果我们想取下n-1,那么n-2要为未被取下来的状态且1~n-3为被取下的状态。

但是这里不能直接将n-2变成为未取下的状态,你想把n-2的状态改变,就得使n-3为未被取下来的状态且1~n-4为被取下的状态;你想把n-3的状态改变,就得使n-4为未被取下来的状态且1~n-5为被取下的状态.

这个一直递推下去你要使1~n-2全部属于未被取下的状态才能继续取n-1。

将1~n全部取下来有下列四步

第一步:取1~n-2,需要F(n-2)步

第二步:取第n个,需要1步

第三步:恢复1~n-2,需要F(n-2)步 (取下来的最少步数和恢复的最少步数是一样的,它们是对称的)

第四步:将1~n-1个全部取下来

所以:F(n) = F(n-1) + 2*F(n-2) + 1

矩阵快速幂套模板啦!

#include 
#include
#include
#include
using namespace std;typedef long long ll;const int Mod = 200907;const int N = 5;int msize;struct Mat{ int mat[N][N];};Mat operator *(Mat a, Mat b){ Mat c; memset(c.mat, 0, sizeof(c.mat)); for(int k = 0; k < msize; ++k) for(int i = 0; i < msize; ++i) if(a.mat[i][k]) for(int j = 0; j < msize; ++j) if(b.mat[k][j]) c.mat[i][j] = ((ll)a.mat[i][k] * b.mat[k][j] + c.mat[i][j])%Mod; return c;}Mat operator ^(Mat a, int k){ Mat c; memset(c.mat,0,sizeof(c.mat)); for(int i = 0; i < msize; ++i) c.mat[i][i]=1; for(; k; k >>= 1) { if(k&1) c = c*a; a = a*a; } return c;}int main(){ int n; msize = 3; Mat A; A.mat[0][0] = 1, A.mat[0][1] = 2, A.mat[0][2] = 1; A.mat[1][0] = 1, A.mat[1][1] = 0, A.mat[1][2] = 0; A.mat[2][0] = 0, A.mat[2][1] = 0, A.mat[2][2] = 1; while(~scanf("%d",&n) && n) { if(n==1) { puts("1"); continue; } Mat ans = A^(n-2); printf("%d\n", (ans.mat[0][0]*2 + ans.mat[0][1] + ans.mat[0][2])%Mod); } return 0;}

 

转载于:https://www.cnblogs.com/pach/p/7256925.html

你可能感兴趣的文章
springboot系列十 Spring-Data-Redis
查看>>
excel进行矩阵计算
查看>>
基于Android平台的动态生成控件和动态改变控件位置的方法
查看>>
BOM
查看>>
iOS: Block的循环引用
查看>>
css详解1
查看>>
MySQL类型转换
查看>>
HashSet HashMap 源码阅读笔记
查看>>
变量声明提升1
查看>>
随笔2013/2/19
查看>>
Windows Phone的Silverlight Toolkit 安装及其使用
查看>>
DBS:同学录
查看>>
Mysql备份系列(1)--备份方案总结性梳理
查看>>
[CareerCup] 1.6 Rotate Image 翻转图像
查看>>
Python中的画图初体验
查看>>
Java程序员的日常 —— 响应式导航Demo
查看>>
objective-c内存管理基础
查看>>
sap关于价值串的说法(转载)
查看>>
Migration to S/4HANA
查看>>
sed 对目录进行操作
查看>>