Description
轮状病毒有很多变种,所有轮状病毒的变种都是从一个轮状基产生的。一个N轮状基由圆环上N个不同的基原子和圆心处一个核原子构成的,2个原子之间的边表示这2个原子之间的信息通道。如下图所示:
N轮状病毒的产生规律是在一个N轮状基中删去若干条边,使得各原子之间有唯一的信息通道,例如共有16个不同的3轮状病毒,如下图所示:
现给定n(N<=100),编程计算有多少个不同的n轮状病毒
Input
第一行有1个正整数n
Output
计算出的不同的n轮状病毒数输出
Sample Input
3
Sample Output
16
解题思路
求完全图的生成树个数。
需要基尔霍夫矩阵的相关知识。
递推式:F(n) = 3F(n - 1) - F(n - 2) + 2。
因为没有取模,此题需要高精度。
大佬写的详细教程
代码附上:
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<cstring>
#include<string>
#include<cstdlib>
using namespace std;
struct ldx
{
int a[101],len;
void operator = (ldx p)
{
len=p.len;
for(int i=1;i<=len;i++) a[i]=p.a[i];
}
void operator * (int p)
{
for(int i=1;i<=len;i++) a[i]*=p;
for(int i=1;i<len;i++)
if(a[i]>=10) {a[i+1]+=a[i]/10;a[i]%=10;}
while(a[len]>=10) {a[len+1]+=a[len]/10;a[len]%=10;len++;}
}
void operator - (ldx p)
{
for(int i=1;i<=len;i++) a[i]-=p.a[i];
for(int i=1;i<len;i++)
while(a[i]<0) {a[i]+=10;a[i+1]-=1;}
while(a[len]==0) len--;
}
void operator + (int p)
{
a[1]+=p;
for(int i=1;i<len;i++)
if(a[i]>=10) {a[i+1]+=a[i]/10;a[i]%=10;}
while(a[len]>=10) {a[len+1]+=a[len]/10;a[len]%=10;len++;}
}
}f[105];
int main()
{
int n;scanf("%d",&n);
f[2].len=1;f[2].a[1]=5;
f[3].len=2;f[3].a[1]=6;f[3].a[2]=1;
for(int i=4;i<=n;i++)
{
f[i]=f[i-1];f[i]*3;
f[i]-f[i-2];f[i]+2;
}
for(int i=f[n].len;i>=1;i--) printf("%d",f[n].a[i]);
return 0;
}