codeforces 962 F. Simple Cycles Edges(点双连通分量||tarjan缩点)

http://codeforces.com/contest/962/problem/F

题意：
给出一个无向图，问有哪些边只属于一个简单环。

思路：
找到每个点双连通分量，如果一个x个点的点双连通分量正好由x条边构成，那么这些边就都只属于一个简单环

#include
#include
#include
#include
#include
#define N 100005
#define M 200005
using namespace std;


//___________tarjan__________________

struct edge
{
    int to,next;
    edge(){}
    edge(int to,int next)
    {
        this->to=to;
        this->next=next;
    }
}ed[M];
int head[N],lnum,esum,index,top,bccnum;
int dfn[N],low[N],mark[M],bj[N],st[M],belong[M],isCut[N],num[N];
pair<int,int>E[M];
void addline(int from,int to)
{
    ed[lnum]=edge(to,head[from]);
    head[from]=lnum++;
}
void tarjan(int root,int fa){//有自环时不加自环的边
                            //点双连通缩点方法：清空路径，枚举E[]数组中存储的路径，建立双向边。
    dfn[root]=low[root]=++index; //新点初始化
    int child=0;
//初始节点需要两个以上儿子且dfn[root]<=low[v] 才是割点
    for(int i=head[root];~i;i=ed[i].next){ //遍历root指出去的边
        int v=ed[i].to;
        if(mark[i]) continue;
        mark[i]=mark[i^1]=1;
        st[++top]=i;//边入栈，需注意此语句要放在判continue之后
         if(!dfn[v]){ //如果v节点未去过，搜索v节点
            child++;
            tarjan(v,root);
            low[root]=min(low[root],low[v]); //更新low值
            if(dfn[root]<=low[v]){
                isCut[root]=1;          //此点是割点，需注意初始节点要有两个儿子
                bccnum++;//注意这里是N++，建数组时要注意开至少两倍大
                for(;;){
                    int j=st[top--];
            //bj[]数组用来标记节点所属的bcc，割点会改变，无意义
            //E[]存新图的边，esum是其数量，tarjan结束后建双向边
                    if(bj[ed[j].to]!=bccnum){
                        bj[ed[j].to]=bccnum;
                        num[bccnum]++;
                        E[++esum]=make_pair(ed[j].to,bccnum);
                    }
                    if(bj[ed[j^1].to]!=bccnum){
                        bj[ed[j^1].to]=bccnum;
                        num[bccnum]++;
                        E[++esum]=make_pair(ed[j^1].to,bccnum);
                    }
                    belong[(j>>1)+1]=bccnum;//标记边所属的bcc
                    if(i==j)break;
                }
            }
        }
        else low[root]=min(low[root],dfn[v]);
    //与有向图区分，此处else不需要判别v节点是否在栈内
	 }
	if(root==fa && child<2)isCut[root]=0;
	//如果初始节点没有2个以上儿子，标记清0
}

void init()
{
    memset(head,-1,sizeof(head));
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(mark,0,sizeof(mark));
    memset(belong,0,sizeof(belong));
    memset(isCut,0,sizeof(isCut));
	top=0;
	lnum=0;
	index=0;
	bccnum=0;
	esum=0;
}


//___________tarjan__________________

vector<int>ans;
vector<int> v[N];
int n,m,x,y;
int main()
{

	ios::sync_with_stdio(false);
	cin.tie(0);
	init();
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		cin>>x>>y;
		addline(x,y);
		addline(y,x);
	}
	for(int i=1;i<=n;i++)
	{
		if(!dfn[i]) tarjan(i,i);
	}
	for(int i=1;i<=m;i++)
	{
		v[belong[i]].push_back(i);
	}
	for(int i=1;i<=bccnum;i++)
	{
		if(v[i].size()==num[i])//边数等于点数
		{
			for(auto x:v[i])
			{
				ans.push_back(x);
			}
		}
	}
	sort(ans.begin(),ans.end());
	cout<<ans.size()<<endl;
	for(auto x:ans) cout<<x<<" ";

}