【BZOJ 2323】 2323: [ZJOI2011]细胞 (DP+矩阵乘法+快速幂*)
2323: [ZJOI2011]细胞
Description
2222年,人类在银河系外的某颗星球上发现了生命,并且携带了一个细胞回到了地球。经过反复研究,人类已经完全掌握了这类细胞的发展规律:
这种细胞最初的形态是“长条形”,一端是头,一端是尾,中间是躯干。细胞内部含有一列密码(你可以认为它是这种细胞的DNA)。密码是一个长度为n的数字串,且仅含有1~9这9种数字,沿着细胞的躯干从头到尾排列着。
首先,细胞会经历一次分裂。细胞将沿躯干方向分裂成若干个球体,躯干将退化成丝状物,连接着相邻的球体。在分裂过程中,质量是均匀分布的。换句话说,若分裂成k个球体,每个球体的质量为原来的1/k。然而,密码的分布是不确定的。若分割成k个球体,密码会被切割成k段(每段长度至少为1),并按从头到尾的顺序分布在各个球体中。如图,为其中一种合法的一次分裂:
接下来,细胞会经历二次分裂。对于每个球体,其中会含有一小段密码(注意他是有序的),我们把它看作一个十进制的数T。这个球体会被分割成T个小球体,并排成一排,之间用躯干退化成的丝状物相连接,并且质量仍然是均匀分布的,每个小球体的质量都是原球体的1/T。至此,密码已经发挥了它的作用,便消失了。如图,为二次分裂:
最后,细胞会进行变异。相邻小球体之间的丝状物可能会退化掉,这两个小球体便会以相切的方式直接连接。显然,二次分裂后,除两端外的每个小球体都有两段丝状物与其连接(头尾两端的小球体只有一段丝状物与其相连)。对于每个小球体,必须至少退化一段与其相连的丝状物,否则这个结构不稳定,会继续变异。如图,为一种稳定的变异:
现在,我们想知道,对于一个给定密码的细胞,总共有多少种稳定的结构。两种结构被认为相同,当且仅当他们拥有相同个数的小球体,从头到尾每个小球体的质量相同,并且从头到尾每对相邻小球体之间的连接方式相同(都是通过丝状物相连或都是通过相切直接相连)。你只需要回答这个结果 mod 1000000007即可。
Input
第一行为一个正整数n,表示细胞密码的长度。
第二行共n个数字,为给定的细胞密码,中间没有空格。
Output
只包含一个整数,为细胞的种数 mod 1000000007的结果。
Sample Input
【样例输入一】1
1
【样例输入二】
1
5
【样例输入三】
2
11
Sample Output
【样例输出一】0
【样例输出二】
3
【样例输出三】
56
HINT
【数据规模】
对于5%的数据满足,n ≤ 6;
对于25%的数据满足,n ≤ 25;
对于60%的数据满足,n ≤ 100;
对于70%的数据满足,n ≤ 300;
对于100%的数据满足,n ≤ 1 000。
Source
【分析】
好题?
感谢数学试卷让我发现了这个可爱的斐波那契数列。。【hhh
然后【矩阵的指数不能mod phi[p]? 。。被奥爷爷d了。。。。
首先,第一步或第三步不同,则出来的结果一定不同【自己想为什么吧、、我就是这样觉得的
对于第三步,他说一个球两边一定有一边被缩了,其实就是n条边,让你选若干条不缩的边,他们不相邻(因为两端的点是有一条边,所以两端的边是不能选的)
这个模型就是数学试卷上的经典模型?f[n]=f[n-1]+f[n-2],就是枚举第n个选还是不选,就是斐波那契数列。
第三步的判断就搞定了。
对于第一第二步,就是把这个大数字分成很多部分,套上第三步的方案加入答案中,这个可以用DP处理。
f[x][0]表示x前面那条边没有选,
f[x][1]表示x前面那条边选了。
f[i][0]=f[i][0]+f[j][0]*F[nw-3]+f[j][1]*F[nw-2]
f[i][1]=f[i][1]+f[j][0]*F[nw-2]+f[j][1]*F[nw-1]
//F表示斐波那契数列, nw表示j+1~i表示的数。
然后问题来了,F[nw???]怎么求,nw是个极大的数额。。
【一开始我以为矩阵幂的指数可以mod phi 还乐呵呵地打了。。。
其实可以预处理,但是传统的快速幂的话,貌似,挺慢??而且你要高精??
传统的快速幂是以2为底的快速幂,他给你的高精数是10进制下的,不如把快速幂改成10为底的,那么你可以发现,有些过程可以合并,比如你算[l,r]的f,中间已经算了[l,r-1]的部分。
所以预处理就是O(n^2*一个不知道多大的常数)
好像别人打的都挺快的。。。可以看看其他解法:http://www.cnblogs.com/cghAndy/p/6594538.html
总时间复杂度:O(n^2*一个不知道多大的常数)
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define LL long long
#define Mod 1000000007
#define Maxn 1010 int n; struct node
{
LL w[][];
}t[],tt[]; LL a[Maxn],f[Maxn][];
char s[Maxn]; void mul(int z,int x,int y)
{
t[]=tt[];
for(int k=;k<=;k++)
for(int i=;i<=;i++)
for(int j=;j<=;j++)
t[].w[i][j]=(t[].w[i][j]+t[x].w[i][k]*t[y].w[k][j])%Mod;
t[z]=t[];
} void qpow(int x,int b)
{
t[]=tt[];
while(b)
{
if(b&) mul(,x,);
mul(x,x,x);
b>>=;
}
t[x]=t[];
} LL B[Maxn][Maxn][]; LL get_f(int x)
{
t[]=tt[];
t[]=tt[];
for(int i=x;i<=n;i++)
{
qpow(,);
t[]=t[];
qpow(,a[i]);
// mul(4,a[i]-10,4); t[]=t[];
B[x][i][]=t[].w[][];
mul(,,);
B[x][i][]=t[].w[][];
mul(,,);
B[x][i][]=t[].w[][];
}
} void init()
{
tt[].w[][]=;tt[].w[][]=;
tt[].w[][]=;tt[].w[][]=; tt[].w[][]=;tt[].w[][]=;
tt[].w[][]=;tt[].w[][]=; tt[].w[][]=;tt[].w[][]=;
tt[].w[][]=;tt[].w[][]=; t[]=tt[];t[]=tt[];
for(int i=;i<=;i++)
{
mul(i,i-,);
}
} int main()
{
init();
scanf("%d",&n);
scanf("%s",s+);
for(int i=;i<=n;i++) a[i]=s[i]-''; memset(f,,sizeof(f));
f[][]=; int i,j,k;
for(i=;i<=n;i++) get_f(i); for(i=;i<=n;i++)
for(j=;j<i;j++)
{
f[i][]=f[i][]+f[j][]*B[j+][i][];f[i][]%=Mod;
f[i][]=f[i][]+f[j][]*B[j+][i][];f[i][]%=Mod;
f[i][]=f[i][]+f[j][]*B[j+][i][];f[i][]%=Mod;
f[i][]=f[i][]+f[j][]*B[j+][i][];f[i][]%=Mod; }
printf("%d\n",f[n][]);
return ;
}
2017-03-21 21:16:25
【BZOJ 2323】 2323: [ZJOI2011]细胞 (DP+矩阵乘法+快速幂*)的更多相关文章
- ZOJ - 3216:Compositions (DP&矩阵乘法&快速幂)
We consider problems concerning the number of ways in which a number can be written as a sum. If the ...
- 【BZOJ-1009】GT考试 KMP+DP+矩阵乘法+快速幂
1009: [HNOI2008]GT考试 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 2745 Solved: 1694[Submit][Statu ...
- BZOJ-1875 HH去散步 DP+矩阵乘法快速幂
1875: [SDOI2009]HH去散步 Time Limit: 20 Sec Memory Limit: 64 MB Submit: 1196 Solved: 553 [Submit][Statu ...
- Qbxt 模拟赛 Day4 T2 gcd(矩阵乘法快速幂)
/* 矩阵乘法+快速幂. 一开始迷之题意.. 这个gcd有个规律. a b b c=a*x+b(x为常数). 然后要使b+c最小的话. 那x就等于1咯. 那么问题转化为求 a b b a+b 就是斐波 ...
- 洛谷 P4910 帕秋莉的手环 矩阵乘法+快速幂详解
矩阵快速幂解法: 这是一个类似斐波那契数列的矩乘快速幂,所以推荐大家先做一下下列题目:(会了,差不多就是多倍经验题了) 注:如果你不会矩阵乘法,可以了解一下P3390的题解 P1939 [模板]矩阵加 ...
- bzoj 3329: Xorequ【数位dp+矩阵乘法】
注意第一问不取模!!! 因为a+b=a|b+a&b,a^b=a|b-a&b,所以a+b=a^b+2(a&b) x^3x==2x可根据异或的性质以转成x^2x==3x,根据上面的 ...
- 矩阵乘法快速幂 codevs 1732 Fibonacci数列 2
1732 Fibonacci数列 2 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题解 查看运行结果 题目描述 Description 在“ ...
- 矩阵乘法快速幂 codevs 1250 Fibonacci数列
codevs 1250 Fibonacci数列 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题目描述 Description 定义:f0=f1=1 ...
- 矩阵乘法快速幂 cojs 1717. 数学序列
矩阵乘法模板: #define N 801 #include<iostream> using namespace std; #include<cstdio> int a[N][ ...
随机推荐
- Distributed Systems 分布式系统
先来扯淡,几天是14年12月31日了,茫茫然,2014就剩最后一天了.这两天国大都放假,我给自己安排了四篇博客欠账,这就是其中的第一篇,简单介绍一些分布式系统的一些概念和设计思想吧.后面三篇分别是Ne ...
- (64位)本体学习程序(ontoEnrich)系统配置说明文档
1系统环境 64位 Ubuntu 2 第三方依赖库配置 boost_1_44_0 #解压boost_1_44_0.tar.gz 到 /usr/local.如果出现权限问题,请用sudo执行该命令 ta ...
- 《JavaScript 实战》:JavaScript 实现图片切割效果
很久之前就在一个网站的截取相片的功能中看到这个效果,也叫图片裁剪.图片剪切(设置一下也可以做出放大镜等类似的效果).当时觉得很神奇,碍于水平有限,没做出来.前些日子突然想做一个透镜效果,就突然想到了这 ...
- ASP .NET登录界面用户验证码代码
//ASP .NET用户登录界面经常用到验证码代码如下 private void Page_Load(object sender, System.EventArgs e) { // 在此处放置用户代码 ...
- final关键字详解
java中,final关键字可以用来修饰类.方法和变量(包括成员变量和局部变量).下面就从这三个方面来了解一下final关键字的基本用法. 1.修饰类 当用final修饰一个类时,表明这个类不能被继承 ...
- Spring是什么+控制反转和依赖注入
Spring是一个开源框架,是一个轻量级的控制反转(IOC)和面向切面(AOP)的容器框架. 原因: (1)通过控制反转(IOC)达到松耦合,IOC也就是把控制权交出去,在使用中直接得到对象 (2)提 ...
- 短信API——短信验证码
简介 短信服务(Short Message Service.SMS)是指通过调用短信发送API,将指定短信内容发送给指定手机用户. 阿里云短信服务 阿里云短信服务产品介绍:https://www.al ...
- 【Learn】CSS定义
CSS基础语法 本文用于介绍CSS相关的知识,用于记录自己的学习笔记.由于我已经熟悉了部分的HTML,所以相关的概念也不在这里进行描述了,直接写自己的一些心得感悟. 1.CSS规则 CSS是由两个主要 ...
- Web服务器处理动态程序三种方式及Apache配置
模块.CGI.FastCGI三种方式介绍 以PHP脚本为例: 模块方式是指Web服务器通过libphp5.so模块调用PHP服务,模块将相关函数嵌入Web服务请求处理流程,不需要额外解释器进程.注意, ...
- sendEmail实现邮件报警
安装 wget http://caspian.dotconf.net/menu/Software/SendEmail/sendEmail-v1.56.tar.gz 或者点击下载 tar -xf sen ...