【JZOJ 5276】【清华集训2017模拟】神奇的玩具

Description

Solution

我说这题是暴力剪枝你信不信！
事实如此；
剪枝后快如闪电：

标准的暴力很简单，枚举每个商店选不选即可，
剪枝1：因为每个礼物只有两商店卖，如果一个商店不选，另一个就一定选了；
剪枝2：上面那个剪枝保证了每个礼物都会买到，那么每个礼物的代价至少是价格小的那个，这是一定要加上的的，如果代价小的商店不选，那么就得选另一个，代价就一定会增加，这个可以用最优性剪枝；
剪枝3：如果一个礼物只能在一个商店买到，这个商店一定选；

Code

#include <cstdio>
#include <cstdlib>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define min(q,w) ((q)>(w)?(w):(q))
#define max(q,w) ((q)<(w)?(w):(q))
using namespace std;
const int N=105;
int read(int &n)
{
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
int m,n,ans,m1,ansc;
int a[N][2][2];
int B[N][N];
int z1[N],z[N];
int bv[N];
void ss(int q,int e,int as)
{
    if(as>=ans)return;
    if(e>m1)return;
    if(q>m)
    {
        ans=min(ans,as);
        return;
    }
    if(z[q]||z1[q]<1){ss(q+1,e,as);return;}
    z[q]=1;
    ss(q+1,e+1,as);
    z[q]=0;
    fo(i,1,B[q][0])if(!z[B[q][i]])z[B[q][i]]=q,e++;
    ss(q+1,e,as+bv[q]);
    fo(i,1,B[q][0])if(z[B[q][i]]==q)z[B[q][i]]=0;
}
void doit()
{
    int q=0,w;
    fo(i,1,n)
    {
        a[i-q][0][0]=a[i][0][0];a[i-q][0][1]=a[i][0][1];
        a[i-q][1][0]=a[i][1][0];a[i-q][1][1]=a[i][1][1];
        if(1==z1[a[i][0][0]]&&z1[a[i][1][0]]==1)
        {
            q++;
            ansc+=min(a[i][0][1],a[i][1][1]);
            m1--;
            z[a[i][1][0]]=z[a[i][0][0]]=233333;
        }
    }
    n-=q;
    q=0;
    fo(i,1,n)
    {
        a[i-q][0][0]=a[i][0][0];a[i-q][0][1]=a[i][0][1];
        a[i-q][1][0]=a[i][1][0];a[i-q][1][1]=a[i][1][1];
        if((a[i][0][1]<a[i][1][1]&&z[a[i][0][0]])||(a[i][0][1]>a[i][1][1]&&z[a[i][1][0]]))
        {
            q++;
            ansc+=min(a[i][0][1],a[i][1][1]);
            z1[a[i][0][0]]--;
            z1[a[i][1][0]]--;
        }
    }
    n-=q;

}
int main()
{
    int q,w;
    read(n),read(m),read(m1);
    fo(i,1,n)
    {
        read(a[i][0][0]),read(a[i][0][1]);
        read(a[i][1][0]),read(a[i][1][1]);
    }
    q=0;
    ansc=0;
    fo(i,1,n)
    {
        a[i-q][0][0]=a[i][0][0];a[i-q][0][1]=a[i][0][1];
        a[i-q][1][0]=a[i][1][0];a[i-q][1][1]=a[i][1][1];
        if(a[i][0][0]==a[i][1][0])
        {
            ansc+=min(a[i][0][1],a[i][1][1]);
            q++;
            if(!z[a[i][0][0]])m1--;
            z[a[i][0][0]]=233333;
        }
    }
    n-=q;q=0;
    fo(i,1,n)z1[a[i][0][0]]++,z1[a[i][1][0]]++;
    doit();doit();doit();doit();
    if(m1<0){printf("-1\n");return 0;}
    fo(i,1,n)
    {
        if((q=a[i][0][0])<(w=a[i][1][0]))B[q][++B[q][0]]=w;
        else B[w][++B[w][0]]=q;
        if(a[i][0][1]<a[i][1][1])ansc+=a[i][0][1],bv[q]+=a[i][1][1]-a[i][0][1];
        else ansc+=a[i][1][1],bv[w]+=a[i][0][1]-a[i][1][1];
    }
    ans=1e9;
    ss(1,0,0);
    if(ans==1e9)printf("-1\n");
    else printf("%d\n",ans+ansc);
    return 0;
}