HDU6397 Character Encoding

考场上,卧槽这题怎么100+队做出来了,卧槽怎么200+队做出来了,卧槽怎么300+队做出来了。。。还有1010也是,这单调队列预处理DP+线段树不是挺难的吗,怎么最后也300+队做出来了。队友一开始NTT nlogn欢乐超时然后就是优化常数,另一个队友推公式。哎还是自己太菜,做完1010后挂了好久的机划了好久水,没有深入想构造题也没深入想1001,我在下午4点50的时候队友说为撒他的直接减不对,我一看C(k+m-1,m-1)-C(k-m-1-n,m-1),这不是显然不对吗,肯定要容斥一下,然而已经没时间了。。。最后构造和1001都没出。tsz:这不是经典水题吗?菜哭.jpg,最后想一下确实是个水题,就是所有情况减去1个大于等于n的情况,加上2个大于等于n的情况,减去3个大于等于n的情况,大于等于2n之类的不用考虑,反正他的约束同样是大于等于n,那么减去n之后之后随便选中间也会有大于等于n的情况出现,那他在原情况中就是大于等于2n。其实看复杂度也看得出,T=400,nlogn过不去,线性的组合数题,很有可能容斥了。

#include<cstdio>
#include<cstring>
#define maxl 200010
#define mod 998244353

int n,m,k,ans;
int fac[maxl],inv[maxl];

inline int qp(int a,int b)
{
	int ans=1,cnt=a;
	while(b)
	{
		if(b&1)
			ans=(1ll*ans*cnt)%mod;
		cnt=1ll*cnt*cnt%mod;
		b>>=1;
	}
	return ans;
}

inline void init()
{
	fac[0]=1;
	for(int i=1;i<maxl;i++)
		fac[i]=(1ll*fac[i-1]*i)%mod;
	inv[maxl-1]=qp(fac[maxl-1],mod-2);
	for(int i=maxl-2;i>=0;i--)
		inv[i]=(1ll*inv[i+1]*(i+1))%mod;
}

inline void prework()
{
	scanf("%d%d%d",&n,&m,&k);
}

inline int C(int n,int r)
{
	if(n<r)
		return 0;
	return (1ll*fac[n]*inv[n-r]%mod*inv[r])%mod;
}

inline void mainwork()
{
	int f=1,nn=k+m-1,rr=m-1;
	ans=C(k+m-1,m-1);
	for(int i=1;i<=m && nn>=rr;i++)
	{
		f=-f;nn-=n;
		ans=(ans+1ll*f*C(m,i)*C(nn,rr)%mod+mod)%mod;
	}
}

inline void print()
{
	printf("%d\n",ans);
}

int main()
{
	init();
	int t;
	scanf("%d",&t);
	for(int i=1;i<=t;i++)
	{
		prework();
		mainwork();
		print();
	}
	return 0;
}

猜你喜欢

转载自blog.csdn.net/liufengwei1/article/details/81737767