BZOJ 1002: [FJOI2007]轮状病毒 高精度

1002: [FJOI2007]轮状病毒

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 5410  Solved: 2944
[Submit][Status][Discuss]

Description

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

　　N轮状病毒的产生规律是在一个N轮状基中删去若干条边，使得各原子之间有唯一的信息通道，例如共有16个不
同的3轮状病毒，如下图所示

　　现给定n(N<=100)，编程计算有多少个不同的n轮状病毒

Input

　　第一行有1个正整数n

Output

　　计算出的不同的n轮状病毒数输出

Sample Input

3

Sample Output

16

看完想用基尔霍夫矩阵

但是。。。精度有点爆炸啊。。。

结果查题解。。。高精度。。。

这个递推式怎么得我也不会。。。练习练习我最厌恶的高精度吧。。。

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>

using namespace std;

typedef double db;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=110;

struct bignum
{int a[N];}f[N];

bignum operator - (const bignum &x,const bignum &y)
{
	bignum c=x;
	register int i;
	for(i=1;i<=y.a[0];i++)
	{
		c.a[i]-=y.a[i];
		if(c.a[i]<0)c.a[i]+=10,c.a[i+1]--;
	}
	while(0==c.a[c.a[0]])c.a[0]--;
	return c;
}
bignum operator + (const bignum &x,int y)
{
	bignum c=x;
	c.a[1]+=y;
	register int i;
	for(i=1;;++i)
	{
		if(c.a[i]>=10)
		{
	 		c.a[i+1]+=c.a[i]/10;c.a[i]%=10;
			continue;
		}
		c.a[0]=max(c.a[0],i);break;
	}
	return c;
}
bignum operator * (const bignum &x,int y)
{
	bignum c=x;
	register int i;
	for(i=1;i<=c.a[0];++i){c.a[i]*=y;}
	for(i=1;i<c.a[0];++i)
	{
		if(c.a[i]>=10)
		{
			c.a[i+1]+=c.a[i]/10;c.a[i]%=10;
		}
	}
	while(c.a[c.a[0]]>=10)
	{c.a[++c.a[0]]+=c.a[c.a[0]-1]/10;c.a[c.a[0]-1]%=10;}
	return c;
}

int main()
{
	int n=read();
	register int i;
	f[1].a[1]=f[1].a[0]=f[2].a[0]=1;f[2].a[1]=5;
	for(i=3;i<=n;++i)
	{f[i]=(f[i-1]*3+2)-f[i-2];}
	for(i=f[n].a[0];i>0;i--)putchar(f[n].a[i]+'0');puts("");
	return 0;
}