2023牛客OI赛前集训营-提高组(第二场) 集合

题目大意

给定正整数 n n n,求有 n n n个元素的集合 { 1 , 2 , … , n } \{1,2,\dots,n\} {1,2,,n}的所有非空子集和的乘积模 998244353 998244353 998244353后的结果。

1 ≤ n ≤ 200 1\leq n\leq 200 1n200


题解

依题意,非空子集和的取值范围为 [ 1 , n × ( n + 1 ) 2 ] [1,\dfrac{n\times (n+1)}{2}] [1,2n×(n+1)]

考虑 D P DP DP。设 f i , j f_{i,j} fi,j表示有 i i i个元素的集合 { 1 , 2 , … , i } \{1,2,\dots,i\} {1,2,,i}中非空子集和为 j j j的子集数量。那么,可以用 O ( n 3 ) O(n^3) O(n3)的时间复杂度求出 f i , j f_{i,j} fi,j

由扩展欧拉定理,可得在求 f i , j f_{i,j} fi,j时模数为 ϕ ( 998244353 ) \phi(998244353) ϕ(998244353),也就是 998244352 998244352 998244352

那么,答案为

a n s = ∑ i = 1 n × ( n + 1 ) / 2 i f n , i ans=\sum\limits_{i=1}^{n\times (n+1)/2}i^{f_{n,i}} ans=i=1n×(n+1)/2ifn,i

这个式子的模数为 998244353 998244353 998244353

于是,问题就解决了。

时间复杂度为 O ( n 3 ) O(n^3) O(n3)

code

#include
using namespace std;
const long long mod=998244353;
int n,l;
long long ans=1,f[205][40005];
long long mi(long long t,long long v){
	if(!v) return 1;
	long long re=mi(t,v/2);
	re=re*re%mod;
	if(v&1) re=re*t%mod;
	return re;
}
int main()
{
	scanf("%d",&n);
	l=n*(n+1)/2;
	f[0][0]=1;
	for(int i=1;i<=n;i++){
		for(int j=0;j<=l;j++){
			f[i][j]=f[i-1][j];
			if(j>=i) f[i][j]=(f[i][j]+f[i-1][j-i])%(mod-1);
		}
	}
	for(int i=1;i<=l;i++){
		ans=ans*mi(i,f[n][i])%mod;
	}
	printf("%lld",ans);
	return 0;
}

你可能感兴趣的:(题解,题解,c++)