CF559C Gerald and Giant Ches

LIV.CF559C Gerald and Giant Chess

DP只要一与排列组合或是容斥等等东西结合在一起就会变得极其毒瘤……

我们设fi表示:走到第i个黑格子上,且之前没有走到任何一个黑格子时的方案数。

则我们如果将棋盘的右下角看作是第n+1个黑格子,fn+1就是答案。

我们将黑格子按照行优先,如果行相同列优先的顺序排序,这样就明确了DP顺序。

则我们有fi=CXi+YiXii=1i1fjC(XiXj)+(YiYj)XiXj

代码:

#include<bits/stdc++.h>
using namespace std;
const int N=2e5;
const int mod=1e9+7;
int n,fac[200100],inv[200100],f[2010];
pair<int,int>p[2010];
int ksm(int x,int y){
	int z=1;
	for(;y;y>>=1,x=(1ll*x*x)%mod)if(y&1)z=(1ll*z*x)%mod;
	return z;
}
int C(int x,int y){
	return 1ll*fac[x]*inv[y]%mod*inv[x-y]%mod;
}
int main(){
	fac[0]=1;
	for(int i=1;i<=N;i++)fac[i]=(1ll*fac[i-1]*i)%mod;
	inv[N]=ksm(fac[N],mod-2);
	for(int i=N-1;i>=0;i--)inv[i]=(1ll*inv[i+1]*(i+1))%mod;
	scanf("%d%d%d",&p[1].first,&p[1].second,&n),n++,p[1].first--,p[1].second--;
	for(int i=2;i<=n;i++)scanf("%d%d",&p[i].first,&p[i].second),p[i].first--,p[i].second--;
	sort(p+1,p+n+1);
	for(int i=1;i<=n;i++){
		f[i]=C(p[i].first+p[i].second,p[i].first);
		for(int j=1;j<i;j++)if(p[j].first<=p[i].first&&p[j].second<=p[i].second)f[i]=(f[i]-1ll*f[j]*C(p[i].first-p[j].first+p[i].second-p[j].second,p[i].first-p[j].first)%mod+mod)%mod;
	}
	printf("%d\n",f[n]);
	return 0;
}

posted @   Troverld  阅读(44)  评论(0编辑  收藏  举报
编辑推荐:
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
· 【杂谈】分布式事务——高大上的无用知识?
点击右上角即可分享
微信分享提示