标签：Slay int WZY Spire 攻击 强化 牌出 CDY PKUWC2018

题目链接

题意分析

这个题其实不是期望

就是一共有\(C_{2n}^m\)种情况 每一种情况选择\(k\)张牌 然后求最大攻击值的总和

我们考虑

当前抽出了选出了\(i\)张强化牌 \(m-i\)张攻击牌

首先 可以肯定的是 能出强化牌就尽量出强化牌

我们去枚举\(i\)

如果\(i<k\) 那么就出\(i\)张强化牌 \(m-i\)张攻击牌

如果\(i≥k\) 那么就出\(k-1\)张强化牌 \(1\)张攻击牌

\(CDY(i,j)\)表示i张强化牌出\(j\)张 所有方案强化的倍率之和

\(WZY(i,j)\)表示i张攻击牌出\(j\)张 所有方案强化的攻击力之和

二者分别对应\(CDY(i,i)* WZY(m-i,k-i)\)以及\(CDY(i,k-1)* WZY(m-i,1)\)

根据乘法分配律

现在考虑如何计算 \(CDY\)以及\(WZY\)

首先 对于相同数量的牌 由于跟顺序没有关系 所以我们\(sort\)之后贪心选择最大即可

\(f(i,j)\)表示选了\(i\)张强化牌并且最靠前的是第\(j\)张牌

同理 \(g(i,j)\)同理

详细见代码

CODE：

/*-------------OI使我快乐-------------*/
int T,n,m,k,ans;
int cdy[M],wzy[M];
int C[N][N];
int dp[M][M][2],sum[M];
IL int CDY(int x,int y)
{
    if(x<y) return 0;
    if(!y) return C[n][x];
    int res=0;
    for(R int j=x-y+1;j<=n-y+1;++j)
    res=(res+1ll*dp[y][j][0]*C[j-1][x-y]%mod)%mod;//唯一合法的情况 就是把不选的(x-y)个放在(j-1)中
    return res;
}
IL int WZY(int x,int y)
{
    if(x<y) return 0;
    int res=0;
    for(R int j=x-y+1;j<=n-y+1;++j)
    res=(res+1ll*dp[y][j][1]*C[j-1][x-y]%mod)%mod;
    return res;
}
int main()
{
//  freopen(".in","r",stdin);
//  freopen(".out","w",stdout);
    read(T);
    for(R int i=0;i<=3000;++i) C[i][0]=1;
    for(R int i=1;i<=3000;++i)
     for(R int j=1;j<=i;++j)
      C[i][j]=(1ll*C[i-1][j-1]+1ll*C[i-1][j])%mod;//预处理出杨辉三角
//  for(R int i=1;i<=10;++i)
//   for(R int j=1;j<=i;++j)
//    printf("%d%c",C[i][j],(j==i ? '\n':' '));
    while(T--)
    {
        read(n);read(m);read(k);ans=0;
        memset(dp,0,sizeof dp);
        for(R int i=1;i<=n;++i) read(cdy[i]);
        for(R int i=1;i<=n;++i) read(wzy[i]);
        sort(cdy+1,cdy+n+1);
        sort(wzy+1,wzy+n+1);
        for(R int i=1;i<=n;++i)
        {//dp[i][j][0]=f[i][j]
            dp[1][i][0]=cdy[i];
            sum[i]=(sum[i-1]+cdy[i])%mod;
        }//我们利用前缀和优化
        for(R int i=2;i<=n;++i)
        {
            for(R int j=1;j<=n-i+1;++j)
            dp[i][j][0]=1ll*cdy[j]*(sum[n]-sum[j]+mod)%mod;//首先是乘法原理
            for(R int j=1;j<=n;++j)
            sum[j]=(sum[j-1]+dp[i][j][0])%mod;
        }
        for(R int i=1;i<=n;++i)
        {
            dp[1][i][1]=wzy[i];
            sum[i]=(sum[i-1]+wzy[i])%mod;
        }
        for(R int i=2;i<=n;++i)
        {
            for(R int j=1;j<=n-i+1;++j)
            dp[i][j][1]=(1ll*wzy[j]*C[n-j][i-1]%mod+(sum[n]-sum[j]+mod)%mod)%mod;
        //由于这里是加法 所以当前贡献就是当前攻击值*过继的状态数 C(n-j,i-1)    
            for(R int j=1;j<=n;++j)
            sum[j]=(sum[j-1]+dp[i][j][1])%mod;
        }
        for(R int i=0;i<m;++i)
        {
            if(i<k) ans=(ans+1ll*CDY(i,i)*WZY(m-i,k-i)%mod)%mod;
            else ans=(ans+1ll*CDY(i,k-1)*WZY(m-i,1)%mod)%mod;
        }
        printf("%d\n",ans);
    }
//  fclose(stdin);
//  fclose(stdout);
    return 0;
}

HEOI 2019 RP++

标签：Slay,int,WZY,Spire,攻击,强化,牌出,CDY,PKUWC2018
来源： https://www.cnblogs.com/LovToLZX/p/10590688.html

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