ICPC Greater New York Region 2019 J. Unicycles 递推分析计数

很好的一道递推计数题。

比赛时我以为是MST计数。。想多了

首先分析只有一个环的情况有哪些？

1：中间点不在换上，其余所有点构成一个环，有n个方案，中间点连n个外圈点。

2：中间点在换上，大小固定时，环的方案也固定（最多n个环）。

发现方案数一定是n*S（n），（把其中一种环旋转得到n种）。

考虑2种环固定时的情况。

我们发现每个点只有3种连法，连中间，和顺时针连边，逆时针连边。

由于环大小固定环固定，只有除环外的点对方案数有贡献。

设f[i]，表示环外i个点的方案数。

显然f[0]=1,f[1]=3.//这与n无关，环外点固定时方案是固定的。

对于i+1，即环外点多一个的情况。

我们考虑在顺时针末尾加这一个点。

这个点连中点时，其余点连边任意即f[i]个

这个点连顺时针点时，其余点也是f[i]

这个点连逆时针的点时。考虑其逆时针的一个点，它顺时针时会有重复

所以这个的方案数时：f[i]-f[i-1]

综上：f[i+1]=3*f[i]-f[i-1]

 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
const double PI= acos(-1.0);
const int M = 4400+7;
/*
int head[M],cnt;
void init(){cnt=0,memset(head,0,sizeof(head));}
struct EDGE{int to,nxt,val;}ee[M*2];
void add(int x,int y){ee[++cnt].nxt=head[x],ee[cnt].to=y,head[x]=cnt;}
*/
const int mod=100007;
ll f[M],sm[M]; 
int main()
{
	f[0]=1,f[1]=3,sm[1]=5;
	for(int i=2;i<=4100;i++)
	f[i]=(f[i-1]*3-f[i-2]+mod)%mod,sm[i]=(sm[i-1]+f[i])%mod;
	int n;
	cin>>n;
	cout<<sm[n-2]*n%mod<<endl;
	return 0;
}