HDU - 2435

题目链接:HDU - 2435


似乎数据太水,很多可以卡死的代码都能秒过。

思路:
跑完一次最小割之后,整个图被分成了属于点1和属于点n的两块。
我们连边肯定是从点1的快,连向点n的块。于是我们可以枚举加的这条边的两个顶点。

这还不TLE?

我们考虑到第一次最小割才是最花时间的,对于后面连的边,我们相当于是在残量网络上面跑,所以是很快的。但是很多人没有记录残量网络也能过。

注意加边之后需要删掉。


AC代码:

#pragma GCC optimize("-Ofast","-funroll-all-loops")
#include<bits/stdc++.h>
//#define int long long
using namespace std;
const int inf=0x3f3f3f3f;
const int N=110,M=1e6+10;
int n,m,s,t,h[N],vis[N];
int head[N],nex[M],to[M],w[M],tot;
struct node{int to,w,nex;}cl[M];
inline void ade(int a,int b,int c){
	to[++tot]=b; nex[tot]=head[a]; w[tot]=c; head[a]=tot;
}
inline void add(int a,int b,int c){ade(a,b,c);	ade(b,a,0);}
inline int bfs(){
	queue<int> q;	q.push(s);	memset(h,0,sizeof h);	h[s]=1;
	while(q.size()){
		int u=q.front();	q.pop();
		for(int i=head[u];i;i=nex[i]){
			if(w[i]&&!h[to[i]]){
				h[to[i]]=h[u]+1;	q.push(to[i]);
			}
		}
	}
	return h[t];
}
int dfs(int x,int f){
	if(x==t)	return f;	int fl=0;
	for(int i=head[x];i&&f;i=nex[i]){
		if(w[i]&&h[to[i]]==h[x]+1){
			int mi=dfs(to[i],min(w[i],f));
			w[i]-=mi,w[i^1]+=mi,fl+=mi,f-=mi;
		}
	}
	if(!fl)	h[x]=-1;
	return fl;
}
inline int dinic(){
	int res=0;
	while(bfs())	res+=dfs(s,inf);
	return res;
}
void dfs(int x){
	vis[x]=1;
	for(int i=head[x];i;i=nex[i]){
		if(!vis[to[i]]&&w[i])	dfs(to[i]);
	}
}
inline void solve(){
	tot=1;	memset(head,0,sizeof head);	memset(vis,0,sizeof vis);
	cin>>n>>m;	s=1,t=n;
	for(int i=1,a,b,c;i<=m;i++)	scanf("%d %d %d",&a,&b,&c),add(a,b,c);
	int res=dinic(),mx=0;	dfs(1);
	for(int i=2;i<=tot;i++)	cl[i]={to[i],w[i],nex[i]};
	for(int i=2;i<n;i++)	if(vis[i]){
		for(int j=2;j<n;j++)	if(i!=j&&!vis[j]){
			int hi=head[i],hj=head[j];	add(i,j,inf);
			mx=max(mx,dinic());	tot-=2;
			head[i]=hi,head[j]=hj;
			for(int i=1;i<=tot;i++)	to[i]=cl[i].to,w[i]=cl[i].w,nex[i]=cl[i].nex; 
		}
	}
	printf("%d\n",res+mx);
}
signed main(){
	int T;	cin>>T;
	while(T--)	solve();
	return 0;
}
发布了579 篇原创文章 · 获赞 242 · 访问量 4万+

猜你喜欢

转载自blog.csdn.net/weixin_43826249/article/details/104300214
hdu