HDU 6184&& 2017广西邀请赛 Counting Stars(三元环计数)
CountingStars
Time Limit: 4000/2000 MS(Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 309 Accepted Submission(s):82
Problem Description
Little A is an astronomy lover, and he hasfound that the sky was so beautiful!
So he is counting stars now!
There are n stars in the sky, and little A has connected them by mnon-directional edges.
It is guranteed that no edges connect one star with itself, and every two edgesconnect different pairs of stars.
Now little A wants to know that how many different "A-Structure"s arethere in the sky, can you help him?
An "A-structure" can be seen as a non-directional subgraph G, with aset of four nodes V and a set of five edges E.
If V=(A,B,C,D) and E=(AB,BC,CD,DA,AC), we call G as an "A-structure".
It is defined that "A-structure" G1=V1+E1 and G2=V2+E2 are same only in the condition that V1=V2 and E1=E2.
Input
There are no more than 300 test cases.
For each test case, there are 2 positive integers n and m in the first line.
2≤n≤105, 1≤m≤min(2×105,n(n−1)2)
And then m lines follow, in each line there are two positive integers u and v,describing that this edge connects node u and node v.
1≤u,v≤n
∑n≤3×105,∑m≤6×105
Output
For each test case, just output oneinteger--the number of different "A-structure"s in one line.
Sample Input
4 5
1 2
2 3
3 4
4 1
1 3
4 6
1 2
2 3
3 4
4 1
1 3
2 4
Sample Output
1
6
题目意思是说,在一个给定的无向图中,找"A-structure"的个数。
"A-structure"这个题目有解释,V=(A,B,C,D) and E=(AB,BC,CD,DA,AC)
这样的一个点集和边集就是符合条件的。
开始没太看懂题目,僵化地认为必须是一个一个由四个点组成的环+一条对角线
这样怎么也解释不了第二个样例。
只好转头在看一遍题目。
其实题目叫我们找的是两个三元环,这两个三元环共用一条边。
比如说V=(A,B,C,D) and E=(AB,BC,CD,DA,AC)这个里面是
三元环ABC和三元环ADC共用了AC。
样例二给的是一个完全图,一共有六条边。这六条边,每一条边都可以作为像上面AC一样的共边。,都能找到两个三元环。
这样,这道题目其实就用简单的方法就可以过。
我们每次枚举一条边,看它是否能作为共边,如果它是,那数有多少个点能够和这条边构成三元环。最后的答案就是C(cnt,2),(这符合条件的点中任选两个都符合A-structure)
比如说样例二,开始判断AB是否满足条件。发现它满足,C和D都可以与AB边构成三元环,所以以AB为共边的A-structure有C(2.2)个,也就是一个,其他5条边也是类似。
这样的确实能够解决问题,但是我们估计一下算法复杂度,发现时间不够。
所以我们需要优化一下。
我们有两种方法来枚举(具体看代码)
一种是枚举X,一种是枚举Y,
为了降低算法复杂度,我们需要根据两个点集的大小来选择哪种方式
这时候,我们可以选择以Sqrt(m)作为为一个分界点,具体可看代码
#include<iostream>
#include<cstring>
#include<stdio.h>
#include<math.h>
#include<string>
#include<stdio.h>
#include<map>
#include<vector>
#include<deque>
#include<algorithm>
#define maxn 100005
using namespace std;
map<pair<int,int>,int> maps;
vector<int> p[maxn];
int outdeg[maxn],link[maxn],vis[maxn];
int n,m;
long long ans;
void init()
{
ans=0;
maps.clear();
memset(outdeg,0,sizeof (outdeg));
memset(link,0,sizeof (link));
memset(vis,0,sizeof (vis));
for(int i=1;i<=n;i++)
p[i].clear();
}
int main()
{
int i,j,k;
int x,y,z;
int t,a,b,c;
long long sum;
pair<int,int> tmp1,tmp2;
while(scanf("%d%d",&n,&m)!=EOF)
{
init();//清空
t=sqrt(m);
for(i=1;i<=m;i++)
{
scanf("%d%d",&a,&b);
p[a].push_back(b);//存储双向边
p[b].push_back(a);
outdeg[a]++; outdeg[b]++;//统计每个点的度数
tmp1=make_pair(a,b);//把边加入Map里面
tmp2=make_pair(b,a);//方便后面查找
maps[tmp1]=1;//用set,Hash都可以,只要后面能找到边即可
maps[tmp2]=1;
}
for(i=1;i<=n;i++)//开始枚举
{
x=i;
vis[x]=1;
for(j=0;j<p[x].size();j++)
{
y=p[x][j];
link[y]=x;//与之相连的边标记出来
}
for(j=0;j<p[x].size();j++)
{
sum=0;
y=p[x][j];//这里比较重要,此时我们选择的边是XY
if(vis[y])//每个点枚举一次
continue;
if(outdeg[y]<=t)//分界点
{
for(k=0;k<p[y].size();k++)
{
z=p[y][k];//此时选择YZ边
if(link[z]==x)//如果XZ相连
sum++;//可构成三元环,计数
}
}
else
{
for(k=0;k<p[x].size();k++)
{
z=p[x][k];//此时XZ相连
pair<int,int> tmp=make_pair(z,y);
if(maps[tmp]!=0)//如果ZY相连
sum++;//可构成三元环,计数
}
}
ans+=(sum*(sum-1))/2;//最后答案是C(cnt,2);
}
}
printf("%lld\n",ans);
}
return 0;
}