[51Nod1850] 抽卡大赛

link

$solution:$

朴素 $dp$，暴力枚举选择 $i$ 号人的第 $j$ 张卡片，朴素 $dp$ 即可，时间复杂度 $O(n^4)$ 。

考虑对于朴素 $dp$ 的优化，发现其实是一个背包卷积的过程，考虑按 $A$ 值从大到小依次加入，每次维护新的 $P$ 值可以做到 $O(1)$ 。

设计生成函数 $F(x)$ 表示将 $1-n$ 的所有多项式卷在一起的答案，每次只要维护多项式除法与乘法即可，而对于每个多项式都是形如 $ax+b$ 的形式，所以直接暴力维护即可。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define mod 1000000007
#define int long long
using namespace std;
inline int read(){
    int f=,ans=;char c=getchar();
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}
    return f*ans;
}
const int MAXN=;
struct node{
    int val,id,res,p;
}g[MAXN*MAXN];
inline int ksm(int a,int b){
    int ans=;
    while(b){
        if(b&) ans*=a,ans%=mod;
        a*=a,a%=mod;
        b>>=;
    }return ans;
}
int n,tot,F[MAXN],G[MAXN],p[MAXN];
inline void mul(int x1,int x0){
    if(!x1) return;
    for(int i=;i<=n;i++) G[i]+=F[i]*x0,G[i]%=mod,G[i+]+=F[i]*x1,G[i+]%=mod;
    for(int i=;i<=n;i++) F[i]=G[i];memset(G,,sizeof(G));return;
}
int Mod(int x){return ((x%mod)+mod)%mod;}
inline void Div(int x1,int x0){
    if(!x1) return;int Inv=ksm(x1,mod-);
    for(int i=n-;i>=;i--) G[i]=(Inv*F[i+])%mod,F[i]=Mod(F[i]-x0*G[i]);
    for(int i=;i<=n;i++) F[i]=G[i];memset(G,,sizeof(G));return;
    return;
}
bool cmp(node x1,node x2){return x1.val>x2.val;}
int P[MAXN],Ans[MAXN],v[MAXN];
signed main(){
//    freopen("51nod_1850_11_in.txt","r",stdin);
    n=read();int inv100=ksm(,mod-);
    for(int i=;i<=n;i++){
        int num=read();
        int sum=;
        int L=++tot,R=;
        for(int j=;j<=num;j++){
            g[++tot].id=i;
            g[tot].val=read(),g[tot].res=((-read())*inv100)%mod,g[tot].p=read();
            sum+=g[tot].p;
            R=tot;
        }
        for(int j=L;j<=R;j++) g[j].p=(g[j].p*ksm(sum,mod-))%mod;
    }
    for(int i=;i<=n;i++) v[i]=read();
    F[]=;
    sort(g+,g+tot+,cmp);
    for(int i=;i<=tot;i++){
        int id=g[i].id;
        Div(P[id],Mod(-P[id]));
        P[id]+=g[i].p,P[id]=Mod(P[id]);P[id]=Mod(P[id]);
        for(int k=;k<n;k++) Ans[id]+=Mod(Mod(Mod(F[k]*v[k+])*g[i].res)*g[i].p),Ans[id]=Mod(Ans[id]);
        mul(P[id],Mod(-P[id]));
    }
    for(int i=;i<=n;i++) printf("%d\n",Ans[i]);return ;
}