51NOD 1486 大大走格子

 

    把第1个和第n个点加进来之后,设f[i][j]为到了第i个点,且经过的不合法的点的奇偶性是j的方案数。最后答案就是 f[n+2][0] - f[n+2][1] 。

    考虑我们算的其实是至少经过枚举的点的方案数,所以 给奇数加上权 -1,偶数加上权 1之后,只有一个不合法点都没走的方案会贡献1,其他都被消去了(考虑除了第一行的杨辉三角的某一行 奇数列的和 = 偶数列的和)

    

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=200000,maxn=2005,ha=1000000007;
int jc[N+5],ni[N+5],n,ans,f[maxn][2];
struct node{
	int x,y;
	bool operator <(const node &u)const{
		return x==u.x?y<u.y:x<u.x;
	}
}a[maxn];
inline int add(int x,int y){ x+=y; return x>=ha?x-ha:x;}
inline int ksm(int x,int y){ int an=1; for(;y;y>>=1,x=x*(ll)x%ha) if(y&1) an=an*(ll)x%ha; return an;}
inline int C(int x,int y){ return x<y?0:jc[x]*(ll)ni[y]%ha*(ll)ni[x-y]%ha;}

inline void init(){
	jc[0]=1;
	for(int i=1;i<=N;i++) jc[i]=jc[i-1]*(ll)i%ha;
	ni[N]=ksm(jc[N],ha-2);
	for(int i=N;i;i--) ni[i-1]=ni[i]*(ll)i%ha;
}

inline void solve(){
	f[1][1]=1;
	for(int i=2,now;i<=n;i++)
	    for(int j=1;j<i;j++) if(a[j].y<=a[i].y){
	    	now=C(a[i].x+a[i].y-a[j].x-a[j].y,a[i].x-a[j].x);
	    	f[i][0]=add(f[i][0],f[j][1]*(ll)now%ha);
	    	f[i][1]=add(f[i][1],f[j][0]*(ll)now%ha);
		}
	ans=add(f[n][0],ha-f[n][1]);
}

int main(){
	int uu,vv; init();
	scanf("%d%d%d",&uu,&vv,&n);
	for(int i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y);
	a[++n]=(node){1,1},a[++n]=(node){uu,vv};
	sort(a+1,a+n+1),solve();
	printf("%d\n",ans);
	return 0;
}

  

转载于:https://www.cnblogs.com/JYYHH/p/8877612.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值