Accept: 204    Submit: 627

Time Limit: 1000 mSec    Memory Limit : 65536 KB

 Problem Description

n个六边形排成一行,相邻两个六边形共用一条边,如下图所示:

记这个图形的生成树个数为t(n)(由于每条边都是不同的,不存在同构的问题)。那么t(1)=6,t(2)=35……

给出n,求mod 1000000007

 Input

第一行给出数据组数T。

之后每一行给出一个正整数n。

T大约为50000,n<=10^18。

 Output

每组数据输出一行答案。

 Sample Input

2212345678987654321

 Sample Output

41733521876

 Source

这题用矩阵快速幂做,先用dp[i][f]表示处理到第i个矩形,第i个矩形右边那条边会不会不去掉的方案数,那么dp[i][1]=d[i-1][0]+dp[i-1][1],dp[i][0]=4*dp[i-1][1]+5*dp[i-1][0].

我们可以构造矩阵【dp[i][1],dp[i][0],sum[i]  】*A=【dp[i+1][1],dp[i+1][0],sum[i+1]  】.
                   1 4 5
容易求出A=1 5 6
                   0 0 1
然后用矩阵快速幂就行了,这里要注意,这题卡常数,所以要先初始化64个矩阵的乘方,然后再算,而且注意取模不能取太多,不然会超时。

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<string>
#include<algorithm>
using namespace std;
typedef unsigned long long ll;
#define inf 0x7fffffff
#define pi acos(-1.0)
#define MOD 1000000007
struct matrix{
int n,m,i;
ll data[3][3];
}a,b,c,d[99],t; matrix multi(matrix &a,matrix &b){
matrix ans;
memset(ans.data,0,sizeof(ans.data));
ans.data[0][0]=(ans.data[0][0]+a.data[0][0]*b.data[0][0]);
ans.data[0][1]=(ans.data[0][1]+a.data[0][0]*b.data[0][1]);
ans.data[0][2]=(ans.data[0][2]+a.data[0][0]*b.data[0][2]);
ans.data[0][0]=(ans.data[0][0]+a.data[0][1]*b.data[1][0]);
ans.data[0][1]=(ans.data[0][1]+a.data[0][1]*b.data[1][1]);
ans.data[0][2]=(ans.data[0][2]+a.data[0][1]*b.data[1][2]);
ans.data[0][0]=(ans.data[0][0]+a.data[0][2]*b.data[2][0])%MOD;
ans.data[0][1]=(ans.data[0][1]+a.data[0][2]*b.data[2][1])%MOD;
ans.data[0][2]=(ans.data[0][2]+a.data[0][2]*b.data[2][2])%MOD; ans.data[1][0]=(ans.data[1][0]+a.data[1][0]*b.data[0][0]);
ans.data[1][1]=(ans.data[1][1]+a.data[1][0]*b.data[0][1]);
ans.data[1][2]=(ans.data[1][2]+a.data[1][0]*b.data[0][2]);
ans.data[1][0]=(ans.data[1][0]+a.data[1][1]*b.data[1][0]);
ans.data[1][1]=(ans.data[1][1]+a.data[1][1]*b.data[1][1]);
ans.data[1][2]=(ans.data[1][2]+a.data[1][1]*b.data[1][2]);
ans.data[1][0]=(ans.data[1][0]+a.data[1][2]*b.data[2][0])%MOD;
ans.data[1][1]=(ans.data[1][1]+a.data[1][2]*b.data[2][1])%MOD;
ans.data[1][2]=(ans.data[1][2]+a.data[1][2]*b.data[2][2])%MOD; ans.data[2][0]=(ans.data[2][0]+a.data[2][0]*b.data[0][0]);
ans.data[2][1]=(ans.data[2][1]+a.data[2][0]*b.data[0][1]);
ans.data[2][2]=(ans.data[2][2]+a.data[2][0]*b.data[0][2]);
ans.data[2][0]=(ans.data[2][0]+a.data[2][1]*b.data[1][0]);
ans.data[2][1]=(ans.data[2][1]+a.data[2][1]*b.data[1][1]);
ans.data[2][2]=(ans.data[2][2]+a.data[2][1]*b.data[1][2]);
ans.data[2][0]=(ans.data[2][0]+a.data[2][2]*b.data[2][0])%MOD;
ans.data[2][1]=(ans.data[2][1]+a.data[2][2]*b.data[2][1])%MOD;
ans.data[2][2]=(ans.data[2][2]+a.data[2][2]*b.data[2][2])%MOD;
return ans; /*
for(int i=0;i<3;i++){
for(int k=0;k<3;k++){
if(a.data[i][k]>0)
for(int j=0;j<3;j++){
ans.data[i][j]=(ans.data[i][j]+a.data[i][k]*b.data[k][j])%MOD;
}
}
}
return ans;
*/
} matrix fast_mod(ll n){
matrix ans;
ans.n=3;
ans.m=3;
memset(ans.data,0,sizeof(ans.data));
ans.data[0][0]=ans.data[1][1]=ans.data[2][2]=1;
int num=1;
while(n>0){
if(n&1){
ans=multi(ans,d[num]);
}
num++;
n>>=1;
}
printf("%I64d\n",(ans.data[0][2]+5*ans.data[1][2]+6*ans.data[2][2] )%MOD);
} void init()
{
int i,j;
d[1].n=d[1].m=3;
d[1].data[0][0]=1;
d[1].data[0][1]=4;
d[1].data[0][2]=5;
d[1].data[1][0]=1;
d[1].data[1][1]=5;
d[1].data[1][2]=6;
d[1].data[2][0]=0;
d[1].data[2][1]=0;
d[1].data[2][2]=1;
for(i=2;i<64;i++){
d[i]=multi(d[i-1],d[i-1]);
}
} int main()
{
ll n,m,i,j;
int T;
init();
scanf("%d",&T);
while(T--)
{
scanf("%I64d",&n);
if(n==0){
printf("0\n");
continue;
}
fast_mod(n-1);
}
return 0;
}

fzu2198 快来快来数一数的更多相关文章

  1. fzu Problem 2198 快来快来数一数 (快速幂+优化)

    题目链接: Problem  2198  快来快来数一数 题目描述: 给出n个六边形排成一排,a[i]代表i个六边形能组成的生成树个数,设定s[i]等于a[1]+a[2]+a[3]+....+a[i- ...

  2. OpenGL快问快答

    OpenGL快问快答 本文内容主要来自对(http://www.opengl.org/wiki/FAQ)的翻译,随机加入了本人的观点.与原文相比,章节未必完整,含义未必雷同,顺序未必一致.仅供参考. ...

  3. 快读&快写模板【附O2优化】

    快读&快写模板 快读快写,顾名思义,就是提升输入和输出的速度.在这里简单介绍一下几种输入输出的优劣. C++ cin/cout 输入输出:优点是读入的时候不用管数据类型,也就是说不用背scan ...

  4. P1908 逆序对——树状数组&离散化&快读快写の学习

    题目简述: 对于给定的一段正整数序列,逆序对就是序列中 a_i>a_jai​>aj​ 且 i<ji<j 的有序对. 输出序列中逆序对的数目. 知识补充: 树状数组: 这东西就是 ...

  5. HDU 1068 Girls and Boys(最大独立集合 = 顶点数 - 最大匹配数)

    HDU 1068 :题目链接 题意:一些男孩和女孩,给出一些人物关系,然后问能找到最多有多少个人都互不认识. 转换一下:就是大家都不认识的人,即最大独立集合 #include <iostream ...

  6. 卡特兰数 Catalan数 ( ACM 数论 组合 )

    卡特兰数 Catalan数 ( ACM 数论 组合 ) Posted on 2010-08-07 21:51 MiYu 阅读(13170) 评论(1)  编辑 收藏 引用 所属分类: ACM ( 数论 ...

  7. HDU 4160 Dolls (最小路径覆盖=顶点数-最大匹配数)

    Dolls Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submiss ...

  8. catalan 数——卡特兰数(转)

    Catalan数——卡特兰数 今天阿里淘宝笔试中碰到两道组合数学题,感觉非常亲切,但是笔试中失踪推导不出来后来查了下,原来是Catalan数.悲剧啊,现在整理一下 一.Catalan数的定义令h(1) ...

  9. html5 vedio 播放器,禁掉进度条快进快退事件

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  10. JS 实现计算一段文字中的字节数,字母数,数字数,行数,汉字数。

    看到了匹配,第一个想到了用正则表达式,哈哈,果然很方便.不过正则表达式高深莫测!我还没有研究明白啊..目前学了点皮毛.代码如下: <!DOCTYPE html PUBLIC "-//W ...

随机推荐

  1. cookie和session会话技术

    因为http协议是无状态的,也就是说每个客户端访问服务器端资源时,服务器并不知道该客户端是谁,所以需要会话技术识别客户端状态.会话技术是帮助服务器记住客户端状态的. 一次会话的开始是通过浏览器访问某个 ...

  2. 【Linux】Linux进程间通信的几种方式

    一.进程间通信的目的 数据传输:一个进程需要将它的数据发送给另一个进程,发送的数据量在一个字节到几M字节之间 共享数据:多个进程要操作共享数据,一个进程对共享数据 信息传递:一个进程需要向另一个进程发 ...

  3. 记一次centos7重启后docker无法启动的问题

    问题描述 在重新了centos7系统后,docker突然就启动不了了,查看报错信息 [root@localhost ~]# systemctl status docker.service ● dock ...

  4. ctfhub技能树—sql注入—过滤空格

    手注 查询数据库 -1/**/union/**/select/**/database(),2 查询表名 -1/**/union/**/select/**/group_concat(table_name ...

  5. Linux删除文件后磁盘目录不释放

    今天测试oracle数据库的时候,把表空间连带内容和数据文件一并删除了,但是删除之后,查看数据文件不存在了,但是目录的带下没有释放 SQL> drop tablespace users incl ...

  6. 入门OJ:扫雪

    扫雪1 题目描述 大雪履盖了整个城市,市政府要求冬季服务部门尽快将一些街道(列在一份清单中)的积雪清除掉以恢复交通,整个城市由许多交叉路口和街道构成,当然任意两个交叉路口都是直接或间接连通的,清单给出 ...

  7. 机器学习7-模型保存&无监督学习

    模型保存和加载 sklearn模型的保存和加载API from sklearn.externals import joblib 保存:joblib.dump(rf, 'test.pkl') 加载:es ...

  8. Javadoc命令与API

    Javadoc命令与API Javadoc工具会抽取类,成员,方法的文档注释生成与这个类配套的API帮助文档 Javadoc命令和文档注释小公司基本不用,但我们应养成良好的编码习惯,所以还是了解一下 ...

  9. Mybatis Plus 3.4版本之后分页插件的变化

    一.MybatisPlusInterceptor 从Mybatis Plus 3.4.0版本开始,不再使用旧版本的PaginationInterceptor ,而是使用MybatisPlusInter ...

  10. scrapy异步的爬虫框架简单的使用

    scrapy异步的爬虫框架 异步的爬虫框架 高性能的数据解析,持久化存储,全栈数据的爬取,中间件,分布式 框架:就是一个集成好了各种功能且具有很强通用性的一个项目模板. 环境安装: Linux: pi ...