zoukankan      html  css  js  c++  java
  • luoguP4707 重返现世

    收集邮票加强版,每个邮票不是等概率获得的了。

    而且是获得K个,如果把一个全集S集合找出其获得时间集合(显然获得时间两两不同)的话,那么就是第n-k+1大的期望!

    %%%Sooke

    min-max容斥扩展:

    推广到期望:

    只要求后面的东西

    对于集合T,设∑t∈T=SUM,那么,E(min(T))=m/SUM

    所以,只要知道SUM,就可以计算贡献

    所以,不妨把SUM放进状态里,记录贡献次数(就是-1和组合数那坨)

    k=n-k+1之后也很小

    f[p][i][j],当k=p时候,前i个,SUM=j的所有集合的贡献

    i不加入:<-f[p][i-1][j]

    i加入,[i-1][j]转移。p?之前的所有的这样的集合大小都+1了

    就是考虑用组合数来巧妙推出|T|->|T|+1

    然后就可以递推了!

    初值:f[0][0][0]=0,f[k][0][0]=-1,或者手动把i=1的情况做出来也可以

    滚动数组

    #include<bits/stdc++.h>
    #define reg register int
    #define il inline
    #define fi first
    #define se second
    #define mk(a,b) make_pair(a,b)
    #define numb (ch^'0')
    using namespace std;
    typedef long long ll;
    template<class T>il void rd(T &x){
        char ch;x=0;bool fl=false;
        while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
        for(x=numb;isdigit(ch=getchar());x=x*10+numb);
        (fl==true)&&(x=-x);
    }
    template<class T>il void output(T x){if(x/10)output(x/10);putchar(x%10+'0');}
    template<class T>il void ot(T x){if(x<0) putchar('-'),x=-x;output(x);putchar(' ');}
    template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('
    ');}
    
    namespace Miracle{
    const int N=1005;
    const int M=10000+5;
    const int K=14;
    const int mod=998244353;
    int n,k,m;
    int p[N];
    int f[2][K][M];
    int inv[M];
    int ad(int x,int y){
        return x+y>=mod?x+y-mod:x+y;
    }
    int main(){
        rd(n);rd(k);rd(m);
        for(reg i=1;i<=n;++i) rd(p[i]);
        k=n-k+1;
        inv[1]=1;
        for(reg i=2;i<=m;++i) {
            inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod;
        }
        int tmp=0;
        for(reg t=1;t<=k;++t) f[tmp][t][0]=-1;
        for(reg i=1;i<=n;++i){
            tmp^=1;
            memset(f[tmp],0,sizeof f[tmp]);
            for(reg t=1;t<=k;++t){
                for(reg j=0;j<=m;++j){
                    f[tmp][t][j]=f[tmp^1][t][j];
                    if(j>=p[i]){
                        f[tmp][t][j]=ad(f[tmp][t][j],ad(f[tmp^1][t-1][j-p[i]],mod-f[tmp^1][t][j-p[i]]));
                    }
                }
            }
        }
        ll ans=0;
        for(reg j=1;j<=m;++j){
            ans=ad(ans,(ll)f[tmp][k][j]*m%mod*inv[j]%mod);
        }
        cout<<ans;
        return 0;
    }
    
    }
    signed main(){
        Miracle::main();
        return 0;
    }
    
    /*
       Author: *Miracle*
    */
  • 相关阅读:
    C# MQTT M2MQTT
    C# MethodInvoker委托的使用
    linux打包/解压-tar
    linux挂载查看、添加与取消
    IE6、火狐不支持a:visited
    js和jquery中的触发事件
    MySQL5.7.9免安装版配置方法
    mysql数据库在Navicat Premium连接的时候出现1862错误
    mysql max_allowed_packet查询和修改
    服务器是windows时tomcat无法打印所有日志配置修改
  • 原文地址:https://www.cnblogs.com/Miracevin/p/10755638.html
Copyright © 2011-2022 走看看