【2012-12】【组合计数】棋盘游戏

最新推荐文章于 2024-07-26 20:03:18 发布

willinglive

最新推荐文章于 2024-07-26 20:03:18 发布

阅读量464

点赞数

分类专栏：数论|组合数学|数学

本文链接：https://blog.csdn.net/willinglive/article/details/42672899

版权

数论|组合数学|数学专栏收录该内容

50 篇文章 1 订阅

订阅专栏

简单题，坑点多（估计是蒟蒻写麻烦了）

另——zhonghaoxi的程序中小数据会有点问题

解个方程可以把y算出来，然后直接枚举，再求和。

但是不会算复杂度，想到枚举就是因为构造不出来这样的数据

//#define _TEST _TEST
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
/************************************************
Code By willinglive    Blog:http://willinglive.cf
************************************************/
#define rep(i,l,r) for(int i=(l),___t=(r);i<=___t;i++)
#define per(i,r,l) for(int i=(r),___t=(l);i>=___t;i--)
#define MS(arr,x) memset(arr,x,sizeof(arr))
#define LL long long
#define INE(i,u,e) for(int i=head[u];~i;i=e[i].next)
inline const int read()
{int r=0,k=1;char c=getchar();for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;
for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';return k*r;}
/
const LL mod=1000000007;
LL n,m,r,c,k;
LL ans;
/
LL pow_mod(LL a,LL b,LL p){LL t=1;for(;b;b>>=1,a=a*a%p)if(b&1)t=t*a%p;return t;}
LL C(LL n,LL m)
{
	m=max(m,n-m);
	if(n<m) return 0;
	LL a=1,b=1;
	rep(i,m+1,n) a=a*i%mod;
	rep(i,1,n-m) b=b*i%mod;
	return a*pow_mod(b,mod-2,mod)%mod;
}
/
void input()
{
    cin>>n>>m>>r>>c>>k;
}
void solve()
{
	for(LL x=r;x>=0;x-=2)
	{
		LL t1=c,t2=0;
		if(n-2*x) t1=t2=(k-m*x)/(n-2*x);
		for(LL y=t1;y>=t2;y-=2)
		{
			if(m*x+n*y-2*x*y!=k) continue;
			if((y+c)&1) continue;
			if(c-y<0) continue;
			ans+=C(n,x)*C(n+((r-x)>>1)-1,(r-x)>>1)%mod*C(m,y)%mod*C(m+((c-y)>>1)-1,(c-y)>>1)%mod;
			ans%=mod;
			//cout<<x<<" "<<y<<endl;
		}
	}
	cout<<ans<<endl;
}
/
int main()
{
    #ifndef _TEST
    freopen("std.in","r",stdin); freopen("std.out","w",stdout);
    #endif
    input(),solve();
    return 0;
}

但是不会算复杂度，想到枚举就是因为构造不出来这样的数据