1002: [FJOI2007]轮状病毒

Description

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

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

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

Input

　　第一行有1个正整数n

Output

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

Sample Input

3

Sample Output

16

HINT

f[1]=1,f[2]=5,f[n]=3*f[n-1]-f[n-2]+2，然后用高精度求答案

#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;
const int maxm = 100005;
int a[maxm], b[maxm], f[maxm], la, lb, lf;
void cal();
int main()
{
	int n, i, j, k, sum;
	scanf("%d", &n);
	if (n == 1) { printf("1\n");return 0; }
	if (n == 2) { printf("5\n");return 0; }
	a[1] = 1;b[1] = 5;
	la = lb = 1;
	for (i = 3;i <= n;i++)
		cal();
	for (i = lf;i >= 1;i -- )
		printf("%d", f[i]);
	printf("\n");
	return 0;
}
void cal()
{
	int i, j, k, xx = 0;
	memset(f, 0, sizeof(f));
	lf = lb;f[1] = 2;
	for (i = 1;i <= lb;i++)
	{
		f[i] += b[i] * 3 + xx;
		xx = f[i] / 10;
		f[i] %= 10;
	}
	if (xx)	f[++lf] = xx;
	for (i = 1;i <= lf;i++)
	{
		f[i] -= a[i];
		if (f[i] < 0) f[i] += 10, f[i + 1]--;
	}
	while (f[lf] == 0) lf--;
	la = lb;
	for (i = 1;i <= lb;i++)	a[i] = b[i];
	lb = lf;
	for (i = 1;i <= lf;i++)	b[i] = f[i];
}