1. 程式人生 > >BZOJ5467 PKUWC2018Slay the Spire(動態規劃)

BZOJ5467 PKUWC2018Slay the Spire(動態規劃)

  即求所有情況的最大傷害之和。容易發現應該先打強化牌,至少打一張攻擊牌。同樣顯然的是強化牌和攻擊牌都應該按從大到小的順序打。進一步可以發現,只要還有強化牌,就應該使用(當然至少留一次攻擊的機會)。

  於是將強化牌和攻擊牌各自從大到小排序。顯然可以將其分開考慮。對強化牌,設f[i][j]為前i張牌抽到j張並打出的強化倍數之和,則顯然有f[i][j]=f[i-1][j]+f[i-1][j-1]·w[i]。這樣就搞定了強化牌可以打完的情況。同時設g[i]為抽i張打出k-1張的強化倍數之和,dp過程中通過f陣列計算,注意避免重複。對於攻擊牌也進行類似dp。然後列舉兩種牌各抽了幾張合併一下答案即可。注意細節。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 3010
#define P 998244353
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'
||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;} int gcd(int n,int m){return m==0?n:gcd(m,n%m);} int read() { int x=0,f=1;char c=getchar(); while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();} while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
return x*f; } int T,n,m,k,a[N],b[N],f[2][N][N],g[2][N],C[N][N],h[N]; int main() { #ifndef ONLINE_JUDGE freopen("bzoj5467.in","r",stdin); freopen("bzoj5467.out","w",stdout); const char LL[]="%I64d\n"; #else const char LL[]="%lld\n"; #endif T=read(); while (T--) { n=read(),m=read(),k=read();int lim=min(k-1,n); C[0][0]=1; for (int i=1;i<=n;i++) { C[i][0]=C[i][i]=1; for (int j=1;j<n;j++) C[i][j]=(C[i-1][j-1]+C[i-1][j])%P; } memset(f,0,sizeof(f));memset(g,0,sizeof(g));memset(h,0,sizeof(h)); for (int i=1;i<=n;i++) a[i]=read(); for (int i=1;i<=n;i++) b[i]=read(); sort(a+1,a+n+1),reverse(a+1,a+n+1); sort(b+1,b+n+1),reverse(b+1,b+n+1); f[0][0][0]=1; for (int i=1;i<=n;i++) { f[0][i][0]=1; for (int j=1;j<=lim;j++) f[0][i][j]=(f[0][i-1][j]+1ll*f[0][i-1][j-1]*a[i])%P; if (lim==0) for (int j=0;j<=n;j++) g[0][j]=C[n][j]; else for (int j=lim;j<=n;j++) g[0][j]=(g[0][j]+1ll*f[0][i-1][lim-1]*a[i]%P*C[n-i][j-lim])%P; } for (int i=1;i<=n;i++) { for (int j=1;j<=m;j++) f[1][i][j]=(f[1][i-1][j]+f[1][i-1][j-1]+1ll*b[i]*C[i-1][j-1])%P; for (int j=m-k+1;j<=n;j++) g[1][j]=(g[1][j]+(f[1][i-1][j-(m-k)-1]+1ll*b[i]*C[i-1][j-(m-k)-1])%P*C[n-i][m-k])%P; } for (int i=1;i<=min(n,m);i++) for (int j=1;j<=n;j++) h[i]=(h[i]+1ll*b[j]*C[n-j][i-1])%P; /*for (int i=0;i<=n;i++) cout<<f[0][n][i]<<' ';cout<<endl; for (int i=0;i<=n;i++) cout<<g[0][i]<<' ';cout<<endl; for (int i=0;i<=n;i++) cout<<f[1][n][i]<<' ';cout<<endl; for (int i=0;i<=n;i++) cout<<g[1][i]<<' ';cout<<endl; for (int i=0;i<=n;i++) cout<<h[i]<<' ';cout<<endl;*/ int ans=0; for (int i=max(0,m-n);i<=min(n,m);i++) if (i<=lim) ans=(ans+1ll*f[0][n][i]*g[1][m-i])%P;//抽了i張強化牌 全部出完 剩餘m-i張攻擊牌 選m-k張不出 else ans=(ans+1ll*g[0][i]*h[m-i])%P;//抽了i張強化牌 選lim張出 剩餘m-i張攻擊牌 出1張 printf("%d\n",ans); } return 0; }