BZOJ5251 八省联考2018劈配(网络流)

  劈配,匹配,网络流。那么考虑怎么跑网络流。

  先看第一问。首先套路的建出超源超汇。不用想也知道导师向汇连容量为战队人数上限的边。特别地,给出局也建一个点,向汇连容量inf的边(似乎没有必要)。对于一个新学员,假设我们已经知道了之前的学员的最优选择,可以把之前的每名学员和可以选择的导师连边,并由源向学员连容量为1的边。然后对于该名新学员,先只连最优选择的边,如果此时跑出的最大流不等于学员数,则表明这名学员无法选择最优,那么删掉最优边(此时这些边里一定没有流量,可以通过容量改为0实现)并连上次优边,次优边还是不行的话继续,一直这样下去直到最大流等于学员数。第一问就做完了。至于复杂度,O(能过)。

  然后是第二问。对于每个人可以二分答案,然后跑最大流看其是否满足。似乎需要访问网络的历史状态?可持久化网络流!这玩意似乎没法可持久化啊……那就暴力记下来呗。由于有C的限制,这里面的边不会很多。于是就做完了,O(能过)。还有一种做法是先连上该学员可以选择的边,然后从第一名开始依次把最优选择加进去跑,直到无法满足,可能会快不少。

 

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
int read()
{
int x=,f=;char c=getchar();
while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
return x*f;
}
#define N 210
#define S 0
#define T 404
int test,c,n,m,b[N],s[N],p[N<<],cnt[N],t;
int cur[N<<],d[N<<],q[N<<],ans,cho[N];
vector<int> a[N][N];
struct data{int to,nxt,cap,flow;
}edge[N*N<<],hisedge[N][];
void addedge(int x,int y,int z)
{
t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].cap=z,edge[t].flow=,p[x]=t;
t++;edge[t].to=x,edge[t].nxt=p[y],edge[t].cap=,edge[t].flow=,p[y]=t;
}
void addedge(int x,int y,int cap,int flow)
{
t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].cap=cap,edge[t].flow=flow,p[x]=t;
t++;edge[t].to=x,edge[t].nxt=p[y],edge[t].cap=,edge[t].flow=-flow,p[y]=t;
}
void init()
{
for (int i=;i<=m;i++) b[i]=read();
for (int i=;i<=n;i++)
{
for (int j=;j<=m+;j++)
a[i][j].clear();
for (int j=;j<=m;j++)
{
int x=read();
if (x) a[i][x].push_back(j);
}
a[i][m+].push_back(m+);
}
for (int i=;i<=n;i++) s[i]=read();
}
bool bfs()
{
memset(d,,sizeof(d));d[S]=;
int head=,tail=;q[]=S;
do
{
int x=q[++head];
for (int i=p[x];~i;i=edge[i].nxt)
if (d[edge[i].to]==-&&edge[i].flow<edge[i].cap)
{
d[edge[i].to]=d[x]+;
q[++tail]=edge[i].to;
}
}while (head<tail);
return ~d[T];
}
int work(int k,int f)
{
if (k==T) return f;
int used=;
for (int i=cur[k];~i;i=edge[i].nxt)
if (d[k]+==d[edge[i].to])
{
int w=work(edge[i].to,min(f-used,edge[i].cap-edge[i].flow));
edge[i].flow+=w,edge[i^].flow-=w;
if (edge[i].flow<edge[i].cap) cur[k]=i;
used+=w;if (used==f) return f;
}
if (used==) d[k]=-;
return used;
}
void dinic()
{
while (bfs())
{
memcpy(cur,p,sizeof(p));
ans+=work(S,N);
}
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("bzoj5251.in","r",stdin);
freopen("bzoj5251.out","w",stdout);
const char LL[]="%I64d";
#else
const char LL[]="%lld";
#endif
test=read(),c=read();
while (test--)
{
n=read(),m=read();
t=-;
memset(p,,sizeof(p));
memset(cnt,,sizeof(cnt));
init();
for (int i=;i<=m;i++) addedge(n+i,T,b[i]);
addedge(n+m+,T,N);
for (int j=;j<=t;j+=)
if (edge[j].cap) hisedge[][++cnt[]]=edge[j],hisedge[][cnt[]].nxt=edge[j^].to;
ans=;
for (int i=;i<=n;i++)
{
addedge(S,i,);
for (int j=;j<=m+;j++)
{
int tmp=t;
for (int k=;k<a[i][j].size();k++)
addedge(i,n+a[i][j][k],);
dinic();
if (ans==i) {cho[i]=j;break;}
for (int k=tmp+;k<=t;k+=) edge[k].cap=;
}
for (int j=;j<=t;j+=)
if (edge[j].cap) hisedge[i][++cnt[i]]=edge[j],hisedge[i][cnt[i]].nxt=edge[j^].to;
}
/*for (int i=0;i<=n;i++)
{
for (int j=1;j<=cnt[i];j++)
cout<<hisedge[i][j].nxt<<' '<<hisedge[i][j].to<<' '<<hisedge[i][j].cap<<' '<<hisedge[i][j].flow<<endl;
cout<<endl;
}*/
for (int i=;i<=n;i++) printf("%d ",cho[i]);
cout<<endl;
for (int i=;i<=n;i++)
{
int l=,r=i-,add=-;
while (l<=r)
{
int mid=l+r>>;
memset(p,,sizeof(p));
t=-;
addedge(S,i,);
for (int j=;j<=s[i];j++)
for (int k=;k<a[i][j].size();k++)
addedge(i,n+a[i][j][k],);
for (int j=;j<=cnt[mid];j++)
addedge(hisedge[mid][j].nxt,hisedge[mid][j].to,hisedge[mid][j].cap,hisedge[mid][j].flow);
ans=;
dinic();
if (!ans) r=mid-;
else l=mid+,add=mid;
}
printf("%d ",i-add-);
}
cout<<endl;
}
return ;
}
上一篇:Vue项目中跨域的几种方式


下一篇:数据库技术之存储过程设计与实现(二)