2023牛客OI赛前集训营-提高组（第二场）

2023牛客OI赛前集训营-提高组（第二场）

T1

题意，给定正整数$n$,计算$n$个元素的集合$1,2,\cdots ,n$,求所有非空子集和的乘积取模$998244353$后的结果。

其中$1\le n\le 200$

首先，该集合的非空子集有$2^n-1$个，但子集和的最大值只有$\frac{n\ast (n-1)}{2}$

，所以肯定有重复的子集和。我们设$f_{i,j}$表示从前面$i$个数中取到值为$j$的方案数。显而易见有两种情况，选和不选，那么$f_{i,j}=(f_{i-1,j}+(j>=i)\ast f_{i-1,j-i})\%(mod-1)$。

此处记得模$(mod-1)$而不是$mod$，下面给出详细证明。

费马小定理:$a^{p-1}\equiv 1(modp)$

利用同乘性可以推出$a^{n\ast (p-1)}\equiv 1(modp)$

所以在$mod$意义下，想使$x^a\equiv x^b(modp)$，就得让$a\equiv b(mod(p-1))$ ，至此证毕。

时间复杂度为$\mathcal{O}(n^3)$。

#include<bits/stdc++.h>
using namespace std;
long long n,f[500+10][125250+10],maxn=0,ans=1;
const long long mod=998244353;
long long read()
{
    long long s=0,w=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
            w=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
        s=s*10+(ch-'0'),ch=getchar();
    return s*w;
}
long long ksm(long long a,long long b)
{
    long long sum=1;
    while(b)
    {
        if(b&1)
            sum=sum*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return sum;
}
int main()
{
    n=read();
    f[0][0]=1;
    maxn=(n+1)*n/2;
    for(int i=1;i<=n;++i)
        for(int j=0;j<=maxn;++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<=maxn;++i)
        ans=ans*ksm(i,f[n][i])%mod;
    printf("%lld",ans);
    return 0;
}

T2

手推了几组数据就会发现，存在一个区间$[l,r]$，$\sum _{i=l}^r{v}_i>(r-l+1+d)\ast k$，就输出$NO$，将这个式子转换成

$\sum _{i=l}^r(v_i-k)>k\ast d$，用$w_i$表示$v_i-k$，也就是$\sum _{i=l}^r{w}_i>k\ast d$，前面就可以用线段树维护最大子段和。

时间复杂度为$\mathcal{O}(q\log n)$。

#include<bits/stdc++.h>
using namespace std;
int n,q,k,d;
struct node
{
    int l,r;
    long long sum,maxl,maxr,maxn;
}t[5000000+10];
void build(int root,int l,int r)
{
    t[root].l=l;
    t[root].r=r;
    if(l==r)
    {
        t[root].sum=-k;
        return ;
    }
    int mid=(l+r)>>1;
    build(root*2,l,mid);
    build(root*2+1,mid+1,r);
    t[root].sum=t[root*2].sum+t[root*2+1].sum;
}
void add(int root,int l,int r,int pos,int val)
{
    if(l==r)
    {
        t[root].sum+=val;
        t[root].maxn=t[root].maxl=t[root].maxr=max(0ll,t[root].sum);
        return ;
    }
    int mid=(l+r)>>1;
    if(pos<=mid)
        add(root*2,l,mid,pos,val);
    else    
        add(root*2+1,mid+1,r,pos,val);
    t[root].maxl=max(t[root*2].maxl,t[root*2].sum+t[root*2+1].maxl);
    t[root].maxr=max(t[root*2+1].maxr,t[root*2+1].sum+t[root*2].maxr);
    t[root].maxn=max(t[root*2].maxr+t[root*2+1].maxl,max(t[root*2].maxn,t[root*2+1].maxn));
    t[root].sum=t[root*2].sum+t[root*2+1].sum;
}
int read()
{
    int s=0,w=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
            w=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
        s=s*10+(ch-'0'),ch=getchar();
    return s*w;
}
int main()
{   
    n=read(),q=read(),k=read(),d=read();
    build(1,1,n);
    while(q--)
    {
        int x=read(),y=read();  
        add(1,1,n,x,y);
        if(t[1].maxn<=1ll*k*d)
            puts("YES");
        else
            puts("NO");
    }
    return 0;
}

T3、T4

未补。