bzoj1002【fjoi2007】轮状病毒

最新推荐文章于 2019-02-18 23:00:00 发布

AaronPolaris

最新推荐文章于 2019-02-18 23:00:00 发布

阅读量1k

点赞数

分类专栏：高精度

本文链接：https://blog.csdn.net/AaronGZK/article/details/47091495

版权

高精度专栏收录该内容

3 篇文章 0 订阅

订阅专栏

1002: [FJOI2007]轮状病毒

Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 3153 Solved: 1749
[ Submit][ Status][ Discuss]

Description

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

Input

第一行有1个正整数n。

Output

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

Sample Input

Sample Output

根据一些推理可得到f[i]=3*f[i-1]-f[i-2]+2 (证明详见神犇vfk博客http://vfleaking.blog.163.com/blog/static/17480763420119685112649/)

然后用高精度就好了。

其实这道题我就是用来练习高精度了。。。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
using namespace std;
int a[105][105];
int n;
int main()
{
	scanf("%d",&n);
	if (n==1)
	{
		cout<<1<<endl;
		return 0;
	}
	if (n==2)
	{
		cout<<5<<endl;
		return 0;
	}
	memset(a,0,sizeof(a));
	a[1][0]=1;
	a[1][1]=1;
	a[2][0]=1;
	a[2][1]=5;
	F(i,3,n)
	{
		a[i][1]+=2;
		F(j,1,a[i-1][0])
		{
			a[i][j]=a[i][j]+3*a[i-1][j]-a[i-2][j];
			while (a[i][j]<0)
			{
				a[i][j]+=10;
				a[i][j+1]--;
			}
			a[i][j+1]+=a[i][j]/10;
			a[i][j]%=10;
		}
		int t=a[i-1][0]+1;
		if (a[i][t]==0) t--;
		while (a[i][t]>10)
		{
			t++;
			a[i][t]+=a[i][t-1]/10;
			a[i][t-1]%=10;
		}
		a[i][0]=t;
	}
	D(i,a[n][0],1) cout<<a[n][i];
	cout<<endl;
}