当前位置: 代码迷 >> 综合 >> 51nod 1551 集合交易
  详细解决方案

51nod 1551 集合交易

热度:38   发布时间:2023-10-29 06:42:42.0

题意

市场中有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;
}