BZOJ1002: [FJOI2007]轮状病毒

bzoj 1002

Description

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

1

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

2

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

Input

第一行有1个正整数n

Output

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

Sample Input

3

Sample Output

16

solve

我们只要选择连在一起有几个就行了,相当于把n分成很多段连续的圆弧

可以由递推关系来解决
以下面n=4为例

3

上图为一种可行解,弧上的实、虚线可互相替换,十字上的实、虚线可互相替换
有4*4种方法

4

上图,弧上n=3的点集再连上4的,有1*f(3)种的情况

5

上图,相连的是同一点集,显然有2*f(2)种

同理,还有3*f(1)种


于是,有了f(4)=4*4+1*f(3)+2*f(2)+3*f(1)

同样的
f(i)=ii+1f(i-1)+2f(i-2)+3f(i-3)+...+(n-1)*f(1)


这时再将上式改为递推式
f(i)=2*i-1 + f(1) + f(2) + ... + f(i-1) + f(i-1)

怎么来的呢?
还以n=4为例

当前式 常数 f(1) f(2) f(3)
f(1) 1*2-1 0 0 0
f(2) 2*2-1 1 0 0
f(3) 3*2-1 1 1 0
f(4) 4*2-1 1 1 1
ANS(4)(求和) 4^2 3 2 1

然后就可以搞事情了,递推式推到底就AC了
注意高精

#include<cstdio>
#include<cstring>
#include<iostream>
#define MAXN 1000
using namespace std;

int n;
struct BIGNUM{
	int num[MAXN];
	int len;
	BIGNUM(){
		memset(num,0,sizeof(num));
		len=1;
	}
}SUM,f,ANS;

BIGNUM ADD(const BIGNUM &x,const BIGNUM &y){
	int len=max(x.len,y.len);
	int i;
	BIGNUM tmp;
	for(i=0;i<=len;i++){
		tmp.num[i]+=x.num[i]+y.num[i];
		tmp.num[i+1]=tmp.num[i]/10;
		tmp.num[i]%=10;
	}
	if(tmp.num[len])
		len++;
	tmp.len=len;
	return tmp;
}

void PRINT(const BIGNUM & x){
	int i;
	for(i=x.len-1;i>=0;i--)
		printf("%d",x.num[i]);
}

void solve(){
	int i;
	for(i=1;i<=n;i++){
		f.num[0]+=i*2-1;
		f=ADD(f,SUM);
		SUM=ADD(f,SUM);
		ANS=f;
	}
}

int main(){
	scanf("%d",&n);
	solve();
	PRINT(ANS);
	return 0;
}
posted @ 2018-01-21 13:16  zzzc18  阅读(106)  评论(0编辑  收藏  举报