POJ 2289 Jamis's Contact Group

一开始题目意思理解错了，其实是要求每个group最多能有几个人，我去找最大group的编号了。。。囧,size居然被我理解成那个意思。。。

正确理解题目意思后，二分答案，枚举最大值，源点与每个人连容量为1的边，每个人与他可以呆的组连一条边，然后将每个组与汇点连线，容量为枚举的最大值，如果满流就减小最大值，否则扩大最大值

代码：

#include<iostream>
using namespace std;
#define min(x,y) (x)<(y)?(x):(y)
#define inf 99999999
#define MAX 1600
struct node
{
int x,y,c,op,next;
}g[1000000];
int map[MAX][605],dis[MAX],pre[MAX];
int num[MAX],cur[MAX],adj[MAX],in[MAX];
int m,n,S,T,vn,max_flow,e,l,r;
void add(int x,int y,int c)
{
e++;
g[e].x=x; g[e].y=y; g[e].c=c; g[e].op=e+1; g[e].next=adj[x]; adj[x]=e;
e++;
g[e].x=y; g[e].y=x; g[e].c=0; g[e].op=e-1; g[e].next=adj[y]; adj[y]=e;
}
void relabel(int k)
{
int i,mm=vn;
cur[k]=adj[k];
for(i=adj[k];i;i=g[i].next)
{
if(g[i].c>0)
mm=min(mm,dis[g[i].y]);
}
dis[k]=mm+1;
}
void change()
{
int nf=inf,i;
for(i=T;i!=S;i=g[pre[i]].x)
{
nf=min(nf,g[pre[i]].c);
}
for(i=T;i!=S;i=g[pre[i]].x)
{
g[pre[i]].c-=nf;
g[g[pre[i]].op].c+=nf;
}
max_flow+=nf;
}
void sap()
{
int i;
for(i=1;i<=vn;i++)
{
dis[i]=0;
num[i]=0;
cur[i]=adj[i];
}
num[0]=vn;
i=S;
while(dis[S]<vn)
{
for(;cur[i];cur[i]=g[cur[i]].next)
{
if(g[cur[i]].c>0&&dis[i]==dis[g[cur[i]].y]+1)
break;
}
if(!cur[i])
{
num[dis[i]]--;
if(num[dis[i]]==0)
break;
relabel(i);
num[dis[i]]++;
if(i!=S)
i=g[pre[i]].x;
}
else
{
pre[g[cur[i]].y]=cur[i];
i=g[cur[i]].y;
if(i==T)
{
change();
i=S;
}
}
}
}
bool check(int limit)
{
int i,j;
e=0;
memset(adj,0,sizeof(adj));
S=n+m+1;
T=n+m+2;
vn=T;
for(i=1;i<=n;i++)
for(j=1;j<=m;j++)
{
if(map[i][j])
{
add(i,n+j,1);
}
}
for(i=1;i<=n;i++)
add(S,i,1);
for(i=1;i<=m;i++)
{
add(i+n,T,limit);
}
max_flow=0;
sap();
if(max_flow==n)
return true;
else
return false;
}
void partition()
{
int i,mid,ans;
l=0;
r=-1;
for(i=1;i<=m;i++)
if(in[i]>r)
r=in[i];
while(l<=r)
{
mid=(l+r)/2;
if(check(mid))
{
ans=mid;
r=mid-1;
}
else
l=mid+1;
}
printf("%d/n",ans);
}
int main()
{
int i,t;
char name[20],c;
while(scanf("%d%d",&n,&m)&&m&&n)
{
memset(map,0,sizeof(map));
memset(in,0,sizeof(in));
for(i=1;i<=n;i++)
{
scanf("%s",&name);
t=0;
getchar();
while(scanf("%c",&c))
{
if(c==10)
{
map[i][++t]=1;
in[t]++;
t=0;
break;
}
if(c==' ')
{
map[i][++t]=1;
in[t]++;
t=0;
}
else
t=t*10+c-'0';
}
}
partition();
}
return 0;
}