HDU 6184 && 2017广西邀请赛:Counting Stars(三元环)




题意:

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);
	}
}

  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值