AGC028D Chords

AGC028D

在一个有 2 n 2n 2n个点的圆上,有 k k k对点被钦定右边,剩下 2 ( n − k ) 2(n-k) 2(nk)个点任意配对。

问所有的方案中,所有的连通块的个数之和。

n ≤ 300 n\le 300 n300


思维僵化。。。

考虑分别对每个连通块计算它的贡献。每个连通块以它的左右端点来统计。

f l , r f_{l,r} fl,r表示: l l l r r r在一个连通块中,并且 [ l , r ] [l,r] [l,r]内的点连向的点都在 [ l , r ] [l,r] [l,r]内,这样的方案数。

g i g_i gi表示 i i i个点任意配对的方案数。

那么答案为 ∑ i , j f i , j g [ i , j ] 外 未 固 定 的 点 \sum_{i,j} f_{i,j}g_{[i,j]外未固定的点} i,jfi,jg[i,j]

考虑 f i , j f_{i,j} fi,j的转移,可以用任意配对的方案数,减去 i i i j j j不在同一个连通块内的方案数。具体来说枚举 k k k,将区间分成 [ i , k ] [i,k] [i,k] [ k + 1 , j ] [k+1,j] [k+1,j]


using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 610
#define ll long long
#define mo 1000000007
int n,K;
int to[N],c[N];
bool bz[N][N];
int f[N][N],g[N];
int main(){
	scanf("%d%d",&n,&K);
	for (int i=1;i<=K;++i){
		int a,b;
		scanf("%d%d",&a,&b);
		if (a>b) swap(a,b);
		to[a]=b,to[b]=a;
		for (int i=a+1;i<=b;++i)
			for (int j=b;j<=2*n;++j)
				bz[i][j]=1;
		for (int j=a;j<=b-1;++j)
			for (int i=1;i<=a;++i)
				bz[i][j]=1;
	}
	for (int i=1;i<=2*n;++i)
		c[i]=c[i-1]+(to[i]==0);
	g[0]=1;
	for (int i=2;i<=2*n;i+=2)
		g[i]=(ll)g[i-2]*(i-1)%mo;
	for (int l=2*n;l>=1;--l){
		for (int r=l+1;r<=2*n;r+=2){
			if (bz[l][r]) continue;
			ll s=0;
			for (int i=l+1;i<r;i+=2)
				(s+=(ll)f[l][i]*g[c[r]-c[i]])%=mo;
			f[l][r]=(g[c[r]-c[l-1]]-s+mo)%mo;
		}
	}
	ll ans=0;
	for (int i=1;i<=2*n;++i)
		for (int j=i+1;j<=2*n;j+=2)
			if (f[i][j])
			(ans+=(ll)f[i][j]*g[c[2*n]-(c[j]-c[i-1])])%=mo;
	printf("%lld\n",ans);
	return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值