BZOJ 1002 FJOI2007 轮状病毒 递推+高精度
题目大意:轮状病毒基定义如图。求有多少n轮状病毒
这个递推实在是不会……所以我选择了打表找规律
首先执行下面程序
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 110
using namespace std;
struct abcd{
int to,next;
bool ban;
}table[M<<2];
int head[M],tot=1;
int n,ans;
void Add(int x,int y)
{
table[++tot].to=y;
table[tot].next=head[x];
head[x]=tot;
}
int fa[M],v[M],q[M],r,h;
bool BFS()
{
int i;
r=h=0;
memset(v,0,sizeof v);
memset(fa,-1,sizeof fa);
q[++r]=0;
while(r!=h)
{
int x=q[++h];
for(i=head[x];i;i=table[i].next)
if(!table[i].ban)
{
if(table[i].to==fa[x])
continue;
if(v[table[i].to])
return 0;
fa[table[i].to]=x;
v[table[i].to]=1;
q[++r]=table[i].to;
}
}
if(r<=n)
return 0;
return 1;
}
void DFS(int x)
{
if(x+x>tot)
{
if( BFS() )
++ans;
return ;
}
table[x<<1].ban=table[x<<1|1].ban=0;
DFS(x+1);
table[x<<1].ban=table[x<<1|1].ban=1;
DFS(x+1);
}
int main()
{
int i;
while(1)
{
memset(head,0,sizeof head);
tot=1;ans=0;
cin>>n;
for(i=1;i<=n;i++)
Add(0,i),Add(i,0),Add(i,i%n+1),Add(i%n+1,i);
DFS(1);
cout<<ans<<endl;
}
}
够简单。够暴力吧
然后打表。1~14的答案例如以下
1 5 16 45 121 320 841 2205 5776 15125 39601 103680 271441 710645
奇数项
1 16 121 841 5776 39601 271441
开根号得
1 4 11 29 76 199 521
a[i]=a[i-1]*3-a[i-2]
偶数项
5 45 320 2205 15125 103680 710645
除以5得
1 9 64 441 3025 20736 142129
开根号得
1 3 8 21 55 144 377
a[i]=a[i-1]*3-a[i-2]
然后高精度递推即可了
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
struct abcd{
int x[100],cnt;
int& operator [] (int y)
{
return x[y];
}
void operator = (int y)
{
x[1]=y;
cnt=1;
}
}f[100];
abcd operator - (abcd x,abcd &y)
{
int i;
abcd z=f[0];
z.cnt=max(x.cnt,y.cnt);
for(i=1;i<=z.cnt;i++)
{
z[i]+=x[i]-y[i];
if(z[i]<0)
z[i+1]--,z[i]+=10;
}
while(z.cnt&&!z[z.cnt])
z.cnt--;
return z;
}
abcd operator * (abcd &x,abcd &y)
{
int i,j;
abcd z=f[0];
for(i=1;i<=x.cnt;i++)
for(j=1;j<=y.cnt;j++)
z[i+j-1]+=x[i]*y[j],z[i+j]+=z[i+j-1]/10,z[i+j-1]%=10;
z.cnt=x.cnt+y.cnt;
if(!z[z.cnt])
--z.cnt;
return z;
}
abcd operator * (abcd x,int y)
{
int i;
abcd z=f[0];
for(i=1;i<=x.cnt;i++)
z[i]+=x[i]*y,z[i+1]+=z[i]/10,z[i]%=10;
z.cnt=x.cnt;
if(z[z.cnt+1])
++z.cnt;
return z;
}
ostream& operator << (ostream &os,abcd x)
{
int i;
for(i=x.cnt;i;i--)
os<<x[i];
return os;
}
int n;
int main()
{
int i;
cin>>n;
f[1]=1;
f[2]=n&1?4:3;
for(i=3;i+i<=n+1;i++)
f[i]=f[i-1]*3-f[i-2];
cout<<f[n+1>>1]*f[n+1>>1]*(n&1?1:5)<<endl;
}
BZOJ 1002 FJOI2007 轮状病毒 递推+高精度的更多相关文章
- 【BZOJ】1002: [FJOI2007]轮状病毒 递推+高精度
1002: [FJOI2007]轮状病毒 Description 给定n(N<=100),编程计算有多少个不同的n轮状病毒. Input 第一行有1个正整数n. Output 将编程计算出的不同 ...
- 【BZOJ1002】[FJOI2007]轮状病毒 递推+高精度
Description 给定n(N<=100),编程计算有多少个不同的n轮状病毒. Input 第一行有1个正整数n. Output 将编程计算出的不同的n轮状病毒数输出 Sample Inpu ...
- bzoj 1002 [FJOI2007]轮状病毒 高精度&&找规律&&基尔霍夫矩阵
1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 2234 Solved: 1227[Submit][Statu ...
- BZOJ 1002: [FJOI2007]轮状病毒【生成树的计数与基尔霍夫矩阵简单讲解+高精度】
1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 5577 Solved: 3031[Submit][Statu ...
- BZOJ 1002 [FJOI2007]轮状病毒
1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 3106 Solved: 1724[Submit][Statu ...
- 生成树的计数(基尔霍夫矩阵):BZOJ 1002 [FJOI2007]轮状病毒
1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 3928 Solved: 2154[Submit][Statu ...
- 【刷题】BZOJ 1002 [FJOI2007]轮状病毒
Description 轮状病毒有很多变种,所有轮状病毒的变种都是从一个轮状基产生的.一个N轮状基由圆环上N个不同的基原子 和圆心处一个核原子构成的,2个原子之间的边表示这2个原子之间的信息通道.如下 ...
- bzoj 1002 [FJOI2007]轮状病毒——打表找规律
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1002 看 Zinn 的博客:https://www.cnblogs.com/Zinn/p/9 ...
- 1002: [FJOI2007]轮状病毒
1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 2609 Solved: 1450[Submit][Statu ...
随机推荐
- linux下kodi没有声音的解决
前几天,心血来潮,就安装了manjaro的pre3版本,由于是mini kde版本的,就随手安装了kodi,可以用来看视频,听音乐和看图片. 结果在所有插件都折腾好了之后发现,在屏幕的右上角有一个喇叭 ...
- MySql视图笔记(转载)
1. 视图的定义 视图就是从一个或多个表中,导出来的表,是一个虚拟存在的表.视图就像一个窗口(数据展示的窗口),通过这个窗口,可以看到系统专门提供的数据(也可以查看到数据表的全部数据),使 ...
- 使用jquery完成定时弹出广告图片
<script src="../js/jquery-1.8.3.js"></script> <script type="text/javas ...
- hdu 1395(欧拉函数)
2^x mod n = 1 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Tot ...
- SpringBoot优雅的全局异常处理
前言 本篇文章主要介绍的是SpringBoot项目进行全局异常的处理. SpringBoot全局异常准备 说明:如果想直接获取工程那么可以直接跳到底部,通过链接下载工程代码. 开发准备 环境要求 JD ...
- redis入门与集群部署
redis入门 redis入门级教程非常多,如http://www.runoob.com/redis/redis-backup.html,作为入门其实已经十分详细了,主要学习内容有如下几个方面吧 1. ...
- MyEclipse导入外部项目
1,File 2,Preferences 3,General----Existing----next 4,Browse选择要导入的项目---finash 5,导入后可能会出现很多error 检查项目的 ...
- django中日志使用学习记录
在setting中加入以下代码 LOGGING = { 'version': 1, 'disable_existing_loggers': True, 'formatters': { 'verbose ...
- antialiasing
Aliasing发生的原因在于,采样频率过低, 这样就无法重建更接近原来texture的图像, Temporal aliasing在时域上 采的texture 次数太少 不能重建真实 是帧与帧之间,s ...
- 【设计模式】工厂方法(FactoryMethod)模式
看不见PPT的请自行解决DNS污染问题. 相关类的代码: namespace FactoryPatternConsole.Model { public class Address { public s ...