fzu2198 快来快来数一数
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
Sample Output
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].
#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 快来快来数一数的更多相关文章
- fzu Problem 2198 快来快来数一数 (快速幂+优化)
题目链接: Problem 2198 快来快来数一数 题目描述: 给出n个六边形排成一排,a[i]代表i个六边形能组成的生成树个数,设定s[i]等于a[1]+a[2]+a[3]+....+a[i- ...
- OpenGL快问快答
OpenGL快问快答 本文内容主要来自对(http://www.opengl.org/wiki/FAQ)的翻译,随机加入了本人的观点.与原文相比,章节未必完整,含义未必雷同,顺序未必一致.仅供参考. ...
- 快读&快写模板【附O2优化】
快读&快写模板 快读快写,顾名思义,就是提升输入和输出的速度.在这里简单介绍一下几种输入输出的优劣. C++ cin/cout 输入输出:优点是读入的时候不用管数据类型,也就是说不用背scan ...
- P1908 逆序对——树状数组&离散化&快读快写の学习
题目简述: 对于给定的一段正整数序列,逆序对就是序列中 a_i>a_jai>aj 且 i<ji<j 的有序对. 输出序列中逆序对的数目. 知识补充: 树状数组: 这东西就是 ...
- HDU 1068 Girls and Boys(最大独立集合 = 顶点数 - 最大匹配数)
HDU 1068 :题目链接 题意:一些男孩和女孩,给出一些人物关系,然后问能找到最多有多少个人都互不认识. 转换一下:就是大家都不认识的人,即最大独立集合 #include <iostream ...
- 卡特兰数 Catalan数 ( ACM 数论 组合 )
卡特兰数 Catalan数 ( ACM 数论 组合 ) Posted on 2010-08-07 21:51 MiYu 阅读(13170) 评论(1) 编辑 收藏 引用 所属分类: ACM ( 数论 ...
- HDU 4160 Dolls (最小路径覆盖=顶点数-最大匹配数)
Dolls Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submiss ...
- catalan 数——卡特兰数(转)
Catalan数——卡特兰数 今天阿里淘宝笔试中碰到两道组合数学题,感觉非常亲切,但是笔试中失踪推导不出来后来查了下,原来是Catalan数.悲剧啊,现在整理一下 一.Catalan数的定义令h(1) ...
- html5 vedio 播放器,禁掉进度条快进快退事件
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- JS 实现计算一段文字中的字节数,字母数,数字数,行数,汉字数。
看到了匹配,第一个想到了用正则表达式,哈哈,果然很方便.不过正则表达式高深莫测!我还没有研究明白啊..目前学了点皮毛.代码如下: <!DOCTYPE html PUBLIC "-//W ...
随机推荐
- cookie和session会话技术
因为http协议是无状态的,也就是说每个客户端访问服务器端资源时,服务器并不知道该客户端是谁,所以需要会话技术识别客户端状态.会话技术是帮助服务器记住客户端状态的. 一次会话的开始是通过浏览器访问某个 ...
- 【Linux】Linux进程间通信的几种方式
一.进程间通信的目的 数据传输:一个进程需要将它的数据发送给另一个进程,发送的数据量在一个字节到几M字节之间 共享数据:多个进程要操作共享数据,一个进程对共享数据 信息传递:一个进程需要向另一个进程发 ...
- 记一次centos7重启后docker无法启动的问题
问题描述 在重新了centos7系统后,docker突然就启动不了了,查看报错信息 [root@localhost ~]# systemctl status docker.service ● dock ...
- ctfhub技能树—sql注入—过滤空格
手注 查询数据库 -1/**/union/**/select/**/database(),2 查询表名 -1/**/union/**/select/**/group_concat(table_name ...
- Linux删除文件后磁盘目录不释放
今天测试oracle数据库的时候,把表空间连带内容和数据文件一并删除了,但是删除之后,查看数据文件不存在了,但是目录的带下没有释放 SQL> drop tablespace users incl ...
- 入门OJ:扫雪
扫雪1 题目描述 大雪履盖了整个城市,市政府要求冬季服务部门尽快将一些街道(列在一份清单中)的积雪清除掉以恢复交通,整个城市由许多交叉路口和街道构成,当然任意两个交叉路口都是直接或间接连通的,清单给出 ...
- 机器学习7-模型保存&无监督学习
模型保存和加载 sklearn模型的保存和加载API from sklearn.externals import joblib 保存:joblib.dump(rf, 'test.pkl') 加载:es ...
- Javadoc命令与API
Javadoc命令与API Javadoc工具会抽取类,成员,方法的文档注释生成与这个类配套的API帮助文档 Javadoc命令和文档注释小公司基本不用,但我们应养成良好的编码习惯,所以还是了解一下 ...
- Mybatis Plus 3.4版本之后分页插件的变化
一.MybatisPlusInterceptor 从Mybatis Plus 3.4.0版本开始,不再使用旧版本的PaginationInterceptor ,而是使用MybatisPlusInter ...
- scrapy异步的爬虫框架简单的使用
scrapy异步的爬虫框架 异步的爬虫框架 高性能的数据解析,持久化存储,全栈数据的爬取,中间件,分布式 框架:就是一个集成好了各种功能且具有很强通用性的一个项目模板. 环境安装: Linux: pi ...