Newcoder 83 E.无向图的弦（分块）

Description

对于一个无向图G 上的简单环$(simplecycle)C$ 以及$G$ 上的一条边$e$，我们称$e$ 是C 的「弦」当且仅当$e$ 的两端点都在$C$ 上且$e$ 本身不是存在于$C$ 的边。

现在给你一个无向图以及好多个询问，每个询问会给你此图上的其中一个简单环，请回答此简单环有多少条「弦」。

Input

输入共有$1+m+1+q$ 行。

第一行有两个正整数$n,m$，分别代表图的点数和边数。

接下来$M$ 行中的第$ i$ 行有两个整数 $x_i,y_i$，代表此图第$ i$ 条边的两端点。

接着有一个正整数$q$，代表有多少个询问。

最后$q$行中的第$i$行有$k_i+1$个整数，第$1$个整数是$k_i$ 代表第$i $个询问的环的长度，接下来的$k_i $个整数$v_{i,0},v_{i,1},…,v_{i,k_i −1}$ 代表此环的的$k_i$ 个点，对于所有$0\le j<k_i$，点$v_{i,j}$ 和点$v_{i,(j+1)mod{k}_i}$ 是此环中相邻的两个点。

$(3\le n,m\le 3\cdot 10^5,1\le q\le 10^5,x_i<y_i,3\le k_i\le m,\sum _{1\le i\le q}{k}_i\le 3\cdot 10^5)$

Output

总共要输出$Q$行，每个询问个输出一行包含一个整数，代表第 $i$ 个询问所给的环共有多少条「弦」。

Sample Input

7 10
0 1
1 2
2 3
3 4
4 5
5 6
0 6
0 3
3 6
2 5
4
7 0 1 2 3 4 5 6
5 0 3 2 5 6
5 0 1 2 5 6
3 0 3 6

Sample Output

3
1
0
0

Solution

对环分块，假设分块界限为$S=\lfloor \sqrt{\frac{m}{2}}\rfloor$，直接统计两点在环上的边，去掉该环的边数即为该环的弦数

1.若$k_i>S$，该类环至多$O(\sqrt{m})$个，直接枚举每条边判断其是否在环上即可，时间复杂度$O(m\sqrt{m})$

2.若$k_i\le S$，对于每个点$u$，记录包含该点环的编号，枚举$u$的所有邻接点$v$，给$v$点打上$u$的标记，之后依次枚举$u$所在环上的点，看其是否有$u$的标记即可，枚举$v$的累加复杂度$O(m)$，这类环上每个点最多被遍历$O(\sqrt{m})$次，而环上点的个数和是$O(m)$规模，时间复杂度$O(m\sqrt{m})$

注意到一条弦不要重复统计，记录$u$在环上的编号后，寻找的是编号大于$u$编号且被打上$u$标记的$v$点

同时$S$没有取$\lfloor \sqrt{m}\rfloor$的原因是第二部分的常数更大，故降低环上每个点被遍历次数

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
typedef pair<int,int>P;
const int maxn=300005,maxq=100005,S=775;
int n,m,q,x[maxn],y[maxn],ans[maxn],vis[maxn];
vector<int>g[maxn],cir[maxq];
vector<P>pos[maxn];
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=0;i<m;i++)
	{
		scanf("%d%d",&x[i],&y[i]);
		g[x[i]].push_back(y[i]),g[y[i]].push_back(x[i]);
	}
	scanf("%d",&q);
	memset(vis,-1,sizeof(vis));
	for(int l=0;l<q;l++)
	{
		int k;
		scanf("%d",&k);
		ans[l]=-k;
		if(k>S)
		{
			for(int i=0;i<k;i++)
			{
				int temp;
				scanf("%d",&temp);
				vis[temp]=l;
			}
			for(int i=0;i<m;i++)
				if(vis[x[i]]==l&&vis[y[i]]==l)ans[l]++;
		}
		else 
		{
			cir[l].resize(k);
			for(int i=0;i<k;i++)
			{
				scanf("%d",&cir[l][i]);
				pos[cir[l][i]].push_back(P(l,i));
			}
		}
	}
	memset(vis,-1,sizeof(vis));
	for(int u=0;u<n;u++)
	{
		for(int i=0;i<g[u].size();i++)vis[g[u][i]]=u;
		for(int i=0;i<pos[u].size();i++)
		{
			int id=pos[u][i].first;
			for(int j=pos[u][i].second+1;j<cir[id].size();j++)
				if(vis[cir[id][j]]==u)ans[id]++;
		}
	}
	for(int i=0;i<q;i++)printf("%d\n",ans[i]);
	return 0;
}