【CSA49F】【XSY3317】card 博弈论 DP

题目大意

　　不会博弈论的 yww 在和博弈论大师 yxq 玩一个游戏。

　　有 \(n\) 种卡牌，第 \(i\) 种卡牌有 \(b_i\) 张。

　　yww 会先把全部 \(B=\sum_{i=1}^nb_i\) 张卡分红两堆，每堆 \(\frac{B}{2}\) 张。保证 \(B\) 是偶数。

　　他们会轮流从第一堆中取卡牌，每次取一张，yww 先取，直到取完为止。

　　而后他们会轮流从第二堆中取卡牌，每次取一张，yxq 先取，直到取完为止。

　　取完卡牌后，他们会计算本身的得分。假设某人在某一堆中取了 \(x\) 张第 \(i\) 种卡牌，那么就能得到 \(\lfloor\frac{x}{a_i}\rfloor c_i\) 分。

　　每一个人的最终得分是这我的在两堆中的得分之和。

　　yxq 想最小化 yww 的得分。

　　做为一名博弈论大师，yxq 每步都会执行最优策略。

　　yww 不会博弈论，因此请你帮 yww 求出他最多能得到多少分。

　　记 \(A=\sum_{i=1}^na_i,B=\sum_{i=1}^nb_i\)；

　　\(1\leq a_i\leq A\leq 2000,1\leq b_i\leq B\leq 500000,2\mid B,1\leq n\leq 2000,1\leq c_i\leq 3000\)；

题解

　　考虑对于一堆牌，yww 先手，他能得到多少分。

　　对于一种牌 \(i\)，若是 \(b_i\equiv -1 \pmod {2a_i}\)，那么先开始拿这种牌的人能够拿到 \(\lfloor \frac{b_i}{2a_i}\rfloor+1\) 张，其余状况都只能拿到 \(\lfloor\frac{b_i}{2a_i}\rfloor\) 张。

　　记 \(b_i\equiv -1 \pmod {2a_i}\) 的牌为特殊的牌，按照 \(c_i\) 从大到小排序，记为 \(d_1,d_2,\ldots,d_k\)，那么最终先手的得分是 \(\sum\limits_i \lfloor\frac{b_i}{2a_i}\rfloor c_i+\sum\limits_{2\nmid i} c_{d_i}\)，后手的得分是 \(\sum\limits_i \lfloor \frac{b_i}{2a_i}\rfloor c_i+\sum\limits_{2\mid i} c_{d_i}\)。

　　这样就能够设计DP状态了：

　　\(f_{i,j,p1,p2}\) 为前 \(i\) 种牌，第一堆分了 \(j\) 张，第一堆有 \(p1\) 种特殊的牌，第二堆有 \(p2\) 种特殊的牌，yww 的最大收益。

　　转移时枚举第 \(i\) 种牌分多少到第一堆。

　　复杂度为 \(O(B^2)\)。

　　注意到收益只与每种牌分到第一堆的牌数 \(\bmod {2a_i}\) 有关，那么DP的时候就能够只枚举 模 \(2a_i\) 的值就行了。

　　还有一个问题，第一堆牌能不能凑出 \(\frac{B}{2}\) 张？

　　对于一种牌，假设咱们把 \(k\) 张牌放到了第一堆，那么有 \(\lfloor\frac{b_i-k}{2a_i}\rfloor\) 组 \(2a_i\) 张牌能够随意分配。这个东西等于 \(\lfloor\frac{b_i}{2a_i}\rfloor\) 或 \(\lfloor\frac{b_i}{2a_i}\rfloor -1\)。咱们伪装它等于 \(\lfloor\frac{b_i}{2a_i}\rfloor -1\)。

　　只可能有 \(O(\sqrt{A})\) 种不一样的 \(a_i\)，随便DP一下就行了。

　　记 \(g_{i,j}\) 为用了前 \(i\) 种 \(a_i\) 组出 \(j\) 张卡牌，第 \(i\) 种 \(a_i\) 最少要多少份（每份 \(2a_i\) 张）（\(-1\) 表示组不出）。

　　时间复杂度为 \(O(A^2+B\sqrt A)\)

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<functional>
#include<cmath>
#include<vector>
#include<assert.h>
//using namespace std;
using std::min;
using std::max;
using std::swap;
using std::sort;
using std::reverse;
using std::random_shuffle;
using std::lower_bound;
using std::upper_bound;
using std::unique;
using std::vector;
typedef long long ll;
typedef unsigned long long ull;
typedef double db;
typedef std::pair<int,int> pii;
typedef std::pair<ll,ll> pll;
void open(const char *s){
#ifndef ONLINE_JUDGE
    char str[100];sprintf(str,"%s.in",s);freopen(str,"r",stdin);sprintf(str,"%s.out",s);freopen(str,"w",stdout);
#endif
}
void open2(const char *s){
#ifdef DEBUG
    char str[100];sprintf(str,"%s.in",s);freopen(str,"r",stdin);sprintf(str,"%s.out",s);freopen(str,"w",stdout);
#endif
}
int rd(){int s=0,c,b=0;while(((c=getchar())<'0'||c>'9')&&c!='-');if(c=='-'){c=getchar();b=1;}do{s=s*10+c-'0';}while((c=getchar())>='0'&&c<='9');return b?-s:s;}
void put(int x){if(!x){putchar('0');return;}static int c[20];int t=0;while(x){c[++t]=x%10;x/=10;}while(t)putchar(c[t--]+'0');}
int upmin(int &a,int b){if(b<a){a=b;return 1;}return 0;}
int upmax(int &a,int b){if(b>a){a=b;return 1;}return 0;}
const int N=2010;
const int M=500010;
int f[N][2][2][4*N];
int s[N];
int qs[N];
int c[M];
int g[M];
struct info
{
    int a,q,v;
};
info a[N];
int cmp(info a,info b)
{
    return a.v>b.v;
}
int n;
int main()
{
    open("49F");
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d%d%d",&a[i].a,&a[i].q,&a[i].v);
    sort(a+1,a+n+1,cmp);
    for(int i=1;i<=n;i++)
    {
        s[i]=s[i-1]+4*a[i].a;
        qs[i]=qs[i-1]+a[i].q;
    }
    memset(f,0x80,sizeof f);
    f[0][0][0][0]=0;
    for(int i=1;i<=n;i++)
        for(int l1=0;l1<=1;l1++)
            for(int l2=0;l2<=1;l2++)
            {
                int x=a[i].q/(2*a[i].a)-1;
                for(int k=0;k<=a[i].q;k++)
                {
                    int v1=k/(2*a[i].a);
                    int v2=(a[i].q-k)/(2*a[i].a);
                    if(v2<x)
                        break;
                    int temp=(v1+v2)*a[i].v;
                    int p1=l1,p2=l2;
                    if(v1*2*a[i].a+2*a[i].a-1==k)
                    {
                        p1?0:temp+=a[i].v;
                        p1^=1;
                    }
                    if(v2*2*a[i].a+2*a[i].a-1==a[i].q-k)
                    {
                        p2?temp+=a[i].v:0;
                        p2^=1;
                    }
                    for(int j=0;j<=s[i-1];j++)
                        f[i][p1][p2][j+k]=max(f[i][p1][p2][j+k],f[i-1][l1][l2][j]+temp);
                }
            }
    for(int i=1;i<=n;i++)
        if(a[i].q/(2*a[i].a)-1>0)
            c[a[i].a]+=a[i].q/(2*a[i].a)-1;
    memset(g,-1,sizeof g);
    g[0]=0;
    for(int i=1;i<=s[n];i++)
        if(c[i])
        {
            for(int j=0;j<=qs[n];j++)
                if(~g[j])
                    g[j]=0;
                else if(j>=2*i&&(g[j-2*i]>=0&&g[j-2*i]<c[i]))
                    g[j]=g[j-2*i]+1;
                else
                    g[j]=-1;
        }
    int ans=0;
    for(int i=0;i<=s[n]&&i<=qs[n]/2;i++)
        if(~g[qs[n]/2-i])
            for(int i1=0;i1<=1;i1++)
                for(int i2=0;i2<=1;i2++)
                    ans=max(ans,f[n][i1][i2][i]);
    printf("%d\n",ans);
    return 0;
}