「PKUWC2018」Slay the Spire

Posted miracevin

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了「PKUWC2018」Slay the Spire相关的知识,希望对你有一定的参考价值。

对于这种题,感觉可以贪心

所以要先手玩,看看有没有最优决策

发现,如果k张,必定先打强化

关键的话是:“每个强化牌翻的倍数大于1”

意味着,多翻一倍少打一张攻击一定不劣。

如果m张,有i张是强化,选k个,

如果i<k,那么打i张强化,剩下打攻击

如果i>=k,打k-1张强化,1张攻击

发现,决策和分到的强化牌的数量有关系

从这方面dp

F(i,j)表示,分到i张强化,打出j张,所有的翻的倍率(卡牌乘积)的和

G(i,j)表示,分到i张攻击,打出j张,所有方案的(原始)伤害的和

那么答案就是:F(i,i)*G(m-i,k-i)和F(i,k-1)*G(m-i,1)(乘法分配律)

考虑怎么求F,G

有了一些牌,一定打最大的一些个

所以sort

然后,f[i][j]表示,前i个以i结尾,选择了j个强化牌所翻的倍率的和

dp,前缀和优化即可。

g[i][j]表示,前i个以i结尾,选择了j个攻击牌所有方案的伤害的和

同上处理

 

然后枚举强化牌i张,

F,G只需要计算n次,所以不用预处理。每次暴力求,枚举最后一个选择的是谁,剩下位置组合数来搞。

注意,f[i][0]是C(n,i),象征选择的方案数,因为乘法的单位元是1,不能没有。

没有强化牌认为翻1倍。

而g[i][0]就是0啦。

#include<bits/stdc++.h>
#define reg register int
#define il inline
#define numb (ch^‘0‘)
#define int long long
using namespace std;
typedef long long ll;
il void rd(int &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);
}
namespace Miracle{
const int N=3003;
const int mod=998244353;
int f[N][N],g[N][N],c[N][N];
int s[N],t[N];
int a[N],b[N];
int n,m,k;
void clear(){
    
}
bool cmp(int x,int y){
    return x>y;
}
int F(int x,int y){
    if(x<y) return 0;
    //if(x==0&&y==0) return 1;
    if(!y) return c[n][x];
    int ret=0;
    for(reg i=y;i<=n;++i){
        ret=(ret+(ll)f[i][y]*c[n-i][x-y])%mod;
    }
    return ret;
}
int G(int x,int y){
    if(x<y) return 0;
    //if(x==0||y==0) return 0;
    int ret=0;
    for(reg i=y;i<=n;++i){
        ret=(ret+(ll)g[i][y]*c[n-i][x-y])%mod;
    }
    return ret;
}
int main(){
    int T;
    //while(1);
    rd(T);
    c[0][0]=1;
    for(reg i=1;i<=3000;++i){
        c[i][0]=1;
        for(reg j=1;j<=i;++j){
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
        }
    }
    while(T--){
        rd(n);rd(m);rd(k);
        clear();
        for(reg i=1;i<=n;++i) rd(a[i]);
        for(reg j=1;j<=n;++j) rd(b[j]);
        sort(a+1,a+n+1,cmp);sort(b+1,b+n+1,cmp);
        
        memset(s,0,sizeof s);
        memset(t,0,sizeof t);
        t[0]=1;
        s[0]=1;
        for(reg i=1;i<=n;++i){
            for(reg j=i;j>=1;--j){
                f[i][j]=(ll)a[i]*s[j-1]%mod;    
            }
            for(reg j=i;j>=1;--j){
                s[j]=(s[j]+f[i][j])%mod;
            }
        }
        memset(s,0,sizeof s);
        memset(t,0,sizeof t);
        t[0]=1;
        for(reg i=1;i<=n;++i){
            for(reg j=i;j>=1;--j){
                g[i][j]=(s[j-1]+(ll)b[i]*c[i-1][j-1])%mod;
            }
            for(reg j=i;j>=1;--j){
                (t[j]+=t[j-1])%=mod;
                (s[j]+=g[i][j])%=mod;
            }
        }
//        for(reg i=1;i<=n;++i){
//            cout<<" ii "<<a[i]<<endl;
//            for(reg j=1;j<=i;++j){
//                cout<<" i j "<<i<<" "<<j<<" ff "<<f[i][j]<<" gg "<<g[i][j]<<endl;
//            }
//        }
        ll ans=0;
        for(reg i=0;i<m;++i){
            if(i<k){
                ans=(ans+(ll)F(i,i)*G(m-i,k-i)%mod)%mod;
            }else{
                ans=(ans+(ll)F(i,k-1)*G(m-i,1)%mod)%mod;
            }
        }
        printf("%lld
",ans);
    }
    return 0;
}

}
signed main(){
    Miracle::main();
    return 0;
}

/*
   Author: *Miracle*
   Date: 2019/1/9 7:21:30
*/

手玩出规律

列出状态

然后乘法分配律得到答案

然后求F,G

以上是关于「PKUWC2018」Slay the Spire的主要内容,如果未能解决你的问题,请参考以下文章

loj2538. 「PKUWC2018」Slay the Spire

spiredoc设置段落字体大小

spire.doc.free企业能用吗

解决使用spire操作文件出现红字:Evaluation Warning: The document was created with Spire.Doc for JAVA.

PKUWC 2018游记

PKUWC2018游记