HDU 6184 && 2017广西邀请赛:Counting Stars(三元环)
程序员文章站
2022-06-04 15:49:32
...
n个点m条边的无向图,问有多少个A-structure
其中A-structure满足V=(A,B,C,D) && E=(AB,BC,CD,DA,AC)
可以看出A-structure是由两个有公共边的三元环构成的
msqrt(m)暴力三元环就好了
#include<stdio.h>
#include<vector>
#include<set>
#include<math.h>
#include<algorithm>
using namespace std;
#define LL long long
vector<int> G[100005];
set<LL> st;
int vis[100005], link[100005], out[100005];
int main(void)
{
LL ans, sum;
int n, i, j, m, k, x, y, z, B;
while(scanf("%d%d", &n, &m)!=EOF)
{
B = sqrt(m);
st.clear();
for(i=1;i<=n;i++)
{
G[i].clear();
vis[i] = out[i] = link[i] = 0;
}
for(i=1;i<=m;i++)
{
scanf("%d%d", &x, &y);
G[x].push_back(y), out[x]++;
G[y].push_back(x), out[y]++;
st.insert((LL)x*n+y);
st.insert((LL)y*n+x);
}
ans = 0;
for(i=1;i<=n;i++)
{
x = i;
vis[x] = 1;
for(j=0;j<G[x].size();j++)
link[G[x][j]] = x;
for(j=0;j<G[x].size();j++)
{
sum = 0;
y = G[x][j];
if(vis[y])
continue;
if(out[y]<=B)
{
for(k=0;k<G[y].size();k++)
{
z = G[y][k];
if(link[z]==x)
sum++;
}
}
else
{
for(k=0;k<G[x].size();k++)
{
z = G[x][k];
if(st.find((LL)z*n+y)!=st.end())
sum++;
}
}
ans += sum*(sum-1)/2;
}
}
printf("%lld\n", ans);
}
}
上一篇: 盘点一下蒜末茼蒿怎么做以及它的功效
下一篇: 白领早餐吃什么对身体好 推荐7种食物