ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

CTS2019D1T1 随机立方体

2022-04-12 16:33:32  阅读:139  来源: 互联网

标签:const CTS2019D1T1 int getchar 1ll 随机 立方体 5e6 mod


做题时:难题啊,我不会啊

看题解:woc,简单题啊,我咋不会啊

讲一下卡住的步骤吧,就是在钦定了极大值的位置和选完数之后,如何算一个类似于十字架的方案。

我们考虑最大值,好的,我们递归到了一个子问题,做完了。

CTS2019 都是些什么神仙题啊!

#include<bits/stdc++.h>
using namespace std;
#define inf 1e9
const int N=5e6;
const int maxn=5e6+10;
const int mod=998244353;
inline int read(){
	int x=0,f=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){x=(x<<1)+(x<<3)+c-'0';c=getchar();}
	return x*f;
}
int n,m,l,k,T,f[maxn],g[maxn],c[maxn];
int fac[maxn],ifc[maxn],inv[maxn];
inline int ksm(int x,int y){
	int res=1;
	while(y){
		if(y&1)res=1ll*res*x%mod;
		x=1ll*x*x%mod;y>>=1;
	}return res;
}
int main(){
	inv[1]=fac[0]=ifc[0]=1;
	for(int i=2;i<=N;i++)inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
	for(int i=1;i<=N;i++)fac[i]=1ll*fac[i-1]*i%mod;
	for(int i=1;i<=N;i++)ifc[i]=1ll*ifc[i-1]*inv[i]%mod;
	T=read();
	while(T--){
		n=read(),m=read(),l=read(),k=read();
		int lim=min(min(n,m),l);g[0]=1;
		for(int i=1;i<=lim;i++)
			g[i]=1ll*g[i-1]*(n-i+1)%mod*(m-i+1)%mod*(l-i+1)%mod;
		int all=1ll*n*m%mod*l%mod,now=1,ans=0;
		for(int i=1;i<=lim;i++)
			c[i]=all+mod-1ll*(n-i)*(m-i)%mod*(l-i)%mod,now=1ll*now*c[i]%mod;
		now=ksm(now,mod-2);
		for(int i=lim;i>=1;i--)
			f[i]=1ll*g[i]*now%mod,now=1ll*now*c[i]%mod;
		for(int i=k;i<=lim;i++)
			if((i-k)&1)ans=(ans-1ll*fac[i]*ifc[k]%mod*ifc[i-k]%mod*f[i])%mod;
			else ans=(ans+1ll*fac[i]*ifc[k]%mod*ifc[i-k]%mod*f[i])%mod;
		printf("%d\n",(ans+mod)%mod);
	}
	return 0;
}

标签:const,CTS2019D1T1,int,getchar,1ll,随机,立方体,5e6,mod
来源: https://www.cnblogs.com/syzf2222/p/16135712.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有