题目大意:轮状病毒基定义如图。求有多少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 轮状病毒 递推+高精度的更多相关文章

  1. 【BZOJ】1002: [FJOI2007]轮状病毒 递推+高精度

    1002: [FJOI2007]轮状病毒 Description 给定n(N<=100),编程计算有多少个不同的n轮状病毒. Input 第一行有1个正整数n. Output 将编程计算出的不同 ...

  2. 【BZOJ1002】[FJOI2007]轮状病毒 递推+高精度

    Description 给定n(N<=100),编程计算有多少个不同的n轮状病毒. Input 第一行有1个正整数n. Output 将编程计算出的不同的n轮状病毒数输出 Sample Inpu ...

  3. bzoj 1002 [FJOI2007]轮状病毒 高精度&&找规律&&基尔霍夫矩阵

    1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 2234  Solved: 1227[Submit][Statu ...

  4. BZOJ 1002: [FJOI2007]轮状病毒【生成树的计数与基尔霍夫矩阵简单讲解+高精度】

    1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 5577  Solved: 3031[Submit][Statu ...

  5. BZOJ 1002 [FJOI2007]轮状病毒

    1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 3106  Solved: 1724[Submit][Statu ...

  6. 生成树的计数(基尔霍夫矩阵):BZOJ 1002 [FJOI2007]轮状病毒

    1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 3928  Solved: 2154[Submit][Statu ...

  7. 【刷题】BZOJ 1002 [FJOI2007]轮状病毒

    Description 轮状病毒有很多变种,所有轮状病毒的变种都是从一个轮状基产生的.一个N轮状基由圆环上N个不同的基原子 和圆心处一个核原子构成的,2个原子之间的边表示这2个原子之间的信息通道.如下 ...

  8. bzoj 1002 [FJOI2007]轮状病毒——打表找规律

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1002 看 Zinn 的博客:https://www.cnblogs.com/Zinn/p/9 ...

  9. 1002: [FJOI2007]轮状病毒

    1002: [FJOI2007]轮状病毒 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 2609  Solved: 1450[Submit][Statu ...

随机推荐

  1. linux下kodi没有声音的解决

    前几天,心血来潮,就安装了manjaro的pre3版本,由于是mini kde版本的,就随手安装了kodi,可以用来看视频,听音乐和看图片. 结果在所有插件都折腾好了之后发现,在屏幕的右上角有一个喇叭 ...

  2. MySql视图笔记(转载)

    1.       视图的定义 视图就是从一个或多个表中,导出来的表,是一个虚拟存在的表.视图就像一个窗口(数据展示的窗口),通过这个窗口,可以看到系统专门提供的数据(也可以查看到数据表的全部数据),使 ...

  3. 使用jquery完成定时弹出广告图片

    <script src="../js/jquery-1.8.3.js"></script> <script type="text/javas ...

  4. hdu 1395(欧拉函数)

    2^x mod n = 1 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Tot ...

  5. SpringBoot优雅的全局异常处理

    前言 本篇文章主要介绍的是SpringBoot项目进行全局异常的处理. SpringBoot全局异常准备 说明:如果想直接获取工程那么可以直接跳到底部,通过链接下载工程代码. 开发准备 环境要求 JD ...

  6. redis入门与集群部署

    redis入门 redis入门级教程非常多,如http://www.runoob.com/redis/redis-backup.html,作为入门其实已经十分详细了,主要学习内容有如下几个方面吧 1. ...

  7. MyEclipse导入外部项目

    1,File 2,Preferences 3,General----Existing----next 4,Browse选择要导入的项目---finash 5,导入后可能会出现很多error 检查项目的 ...

  8. django中日志使用学习记录

    在setting中加入以下代码 LOGGING = { 'version': 1, 'disable_existing_loggers': True, 'formatters': { 'verbose ...

  9. antialiasing

    Aliasing发生的原因在于,采样频率过低, 这样就无法重建更接近原来texture的图像, Temporal aliasing在时域上 采的texture 次数太少 不能重建真实 是帧与帧之间,s ...

  10. 【设计模式】工厂方法(FactoryMethod)模式

    看不见PPT的请自行解决DNS污染问题. 相关类的代码: namespace FactoryPatternConsole.Model { public class Address { public s ...