51nod 1551 集合交易

题意

市场中有n个集合在卖。我们想买到满足以下要求的一些集合,所买到集合的个数要等于所有买到的集合合并后的元素的个数。
每个集合有相应的价格,要使买到的集合花费最小。
这里我们的集合有一个特点:对于任意整数k(k>0),k个集合的并集中,元素的个数不会小于k个。
现在让你去市场里买一些满足以上条件集合,可以一个都不买。

题解

这题最后一步没有做出来是很可惜的。。
首先,我们发现他的特点很神奇。。
于是我就思考这个有什么用。。
容易想到Hall定理
于是就猜想会不会先给每一个点匹配一下
然后怎么办呢?
我们就可以把每一个点看做是一个集合
你可以得到选一个点就可以选很多别的点
这就是很多个约束条件
但是我并不觉得这些条件是够的
在没有题解的情况下我主观地人为这个结论是对地,然后就没有往这个方向走了
然后想了一个看起来有可能对的,然后证了一下,是错了。。

然后看了题解
发现这些约束条件就是够的

然后尝试证明一下,发现还挺好证明这个东西的充分性

于是就没了

CODE:

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
const int MAX=(1<<30);
const int N=305;
int n;
bool a[N][N];
int g[N],G[N];
bool vis[N];
bool dfs (int x)
{
    if (vis[x]) return false;
    vis[x]=true;
    for (int u=1;u<=n;u++)
        if (a[x][u])
        {
            if (g[u]==-1||dfs(g[u]))
            {
                g[u]=x;
                return true;
            }
        }
    return false;
}
struct qq{int x,y,z,last;}e[N*N*10];int num,last[N];
int st,ed;
void init (int x,int y,int z)
{
//  printf("%d %d\n",x,y);
    num++;
    e[num].x=x;e[num].y=y;e[num].z=z;
    e[num].last=last[x];
    last[x]=num;

    num++;
    swap(x,y);z=0;
    e[num].x=x;e[num].y=y;e[num].z=z;
    e[num].last=last[x];
    last[x]=num;
}
int h[N];
bool Bt ()
{
    memset(h,-1,sizeof(h));h[st]=1;
    queue<int> q;
    q.push(st);
    while (!q.empty())
    {
        int x=q.front();q.pop();
        for (int u=last[x];u!=-1;u=e[u].last)
        {
            int y=e[u].y;
            if (e[u].z>0&&h[y]==-1)
            {
                h[y]=h[x]+1;
                q.push(y);
            }
        }
    }
    return h[ed]!=-1;
}
int dfs (int x,int f)
{
    if (x==ed) return f;
    int s1=0;
    for (int u=last[x];u!=-1;u=e[u].last)
    {
        int y=e[u].y;
        if (e[u].z>0&&h[y]==h[x]+1&&s1<f)
        {
            int lalal=dfs(y,min(e[u].z,f-s1));
            s1+=lalal;
            e[u].z-=lalal;
            e[u^1].z+=lalal;
        }
    }
    if (s1==0) h[x]=-1;
    return s1;
}
void dfs1 (int x)
{
    if (vis[x]) return ;
    vis[x]=true;
    x=g[x];
    for (int u=1;u<=n;u++)
        if (a[x][u])
            dfs1(u);
}
int s[N];
int main()
{
    num=1;memset(last,-1,sizeof(last));
    memset(a,false,sizeof(a));
    scanf("%d",&n);st=n+1;ed=st+1;
    for (int u=1;u<=n;u++)
    {
        int x;
        scanf("%d",&x);
        for (int i=1;i<=x;i++)  {int xx;scanf("%d",&xx);a[u][xx]=true;}
    }
    for (int u=1;u<=n;u++) 
    {
        scanf("%d",&s[u]);
        s[u]=-s[u];
    }
    memset(g,-1,sizeof(g));
    for (int u=1;u<=n;u++)  
    {
        memset(vis,false,sizeof(vis));
        dfs(u);
    }
    for (int u=1;u<=n;u++)
    {
        memset(vis,false,sizeof(vis));
        dfs1(u);
        for (int i=1;i<=n;i++)
            if (vis[i]&&u!=i)
                init(u,i,MAX);
    //  printf("%d\n",u);
    }
    //printf("YES\n");
    int ans=0;
    for (int u=1;u<=n;u++)
    {
        int x=g[u];
        if (s[x]<0) init(u,ed,-s[x]);
        if (s[x]>=0) 
        {
            init(st,u,s[x]);
            ans=ans+s[x];
        }
    }
    while (Bt()) ans=ans-dfs(st,MAX);
    printf("%d\n",-ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值