Codeforces Round #532 (Div. 2) E. Andrew and Taxi(拓扑/dfs判环)

 

题意

给定一个有环图，把这个图变成DAG

可以将一些边反向，反向的cost为边权

多条边反向为这些边边权的最大值

求如何反向使cost最小，输出cost、边数和边集

题解

二分一个值ans，先不管比这个小的边权的边，

取只有比这个大的边权的边和所有顶点的生成子图，

二分出生成子图无环的最小值ans，

dfs判环/拓扑判环均可//各写了一个

然后对生成子图拓扑排序，

再考察其余比ans边权小的边，若拓扑逆序则加入rev边集，代表需要反转

各路神仙dfs、手写队列的拓扑、手撸scc的代码看得我真是头皮发麻

思路来源

各路cf神仙的提交代码

代码

#include <iostream>
#include <algorithm> 
#include <cstring>
#include <cstdio>
#include <cmath>
#include <set>
#include <map>
#include <vector>
#include <stack>
#include <queue>
#include <functional>
const int INF=0x3f3f3f3f;
const int maxn=1e5+10; 
const int mod=1e9+7;
const int MOD=998244353;
const double eps=1e-7;
typedef long long ll;
#define vi vector<int> 
#define si set<int>
#define pii pair<int,int> 
#define pi acos(-1.0)
#define pb push_back
#define mp make_pair
#define lowbit(x) (x&(-x))
#define sci(x) scanf("%d",&(x))
#define scll(x) scanf("%I64d",&(x))
#define sclf(x) scanf("%lf",&(x))
#define pri(x) printf("%d",(x))
#define rep(i,j,k) for(int i=j;i<=k;++i)
#define per(i,j,k) for(int i=j;i>=k;--i)
#define mem(a,b) memset(a,b,sizeof(a)) 
using namespace std;
struct edge
{
  int from,to,nex,w,id;	
}e[maxn];
int head[maxn],cnt,n,m,vis[maxn];
vector<int>res; 
int pp;
int que[maxn],s=0,t=0,num,color[maxn];
int in[maxn];
void add(int u,int v,int w,int i)
{
	e[cnt].from=u;
	e[cnt].to=v;
	e[cnt].nex=head[u];
	e[cnt].w=w;
	e[cnt].id=i;
	head[u]=cnt++;
}
void init()
{
	mem(head,-1);
	cnt=0;
}
void dfs(int u,int maxw,bool &flag)
{
	vis[u]=1;
	for(int i=head[u];~i;i=e[i].nex)
	{
	   int w=e[i].w,v=e[i].to;
	   if(w<=maxw)continue;
	   if(vis[v]==1)flag=1;//在栈中且已访问	
	   else if(!vis[v])dfs(v,maxw,flag);//未访问 
	}
	vis[u]=2;//不在栈中且已访问
}
bool cycle(int mid)
{
	bool flag=0;
	mem(vis,0); 
	for(int i=0;i<n&&!flag;++i)
	if(!vis[i])dfs(i,mid,flag);
	return flag;
}

int erfen(int l,int r)
{
	while(l<r)
	{
		int mid=(l+r)/2;
		if(cycle(mid))l=mid+1;
		else r=mid;
	}
	return l;
}
int main()
{
	sci(n),sci(m);
	init();
	rep(i,0,m-1)
	{
		int u,v,w;
		sci(u),sci(v),sci(w);
		u--,v--;
		add(u,v,w,i+1);
	}
	int ans=erfen(0,1e9);//无环的最大权值 
	rep(i,0,n-1)
	{
		for(int j=head[i];~j;j=e[j].nex)
		{
			int v=e[j].to,w=e[j].w;
			if(w<=ans)continue;
			in[v]++;
		}
	}
	rep(i,0,n-1)if(in[i]==0)que[t++]=i;
	while(s<t)
	{
		int tmp=que[s++];
		color[tmp]=++num;//可能会被忽略一个点 
		for(int j=head[tmp];~j;j=e[j].nex)
		{
			int v=e[j].to,w=e[j].w;
			if(w<=ans)continue;
			in[v]--;
			if(in[v]==0)que[t++]=v;
		}
	}
	rep(i,0,m-1)
	{
		int u=e[i].from,v=e[i].to,w=e[i].w,id=e[i].id;
		if(w>ans)continue;
		if(color[u]>color[v])res.push_back(id);
	}
	printf("%d %d\n",ans,pp=res.size());
	rep(i,0,pp-1)
	printf("%d%c",res[i],i==pp-1?'\n':' ');
	return 0;
}