hdu 6184 Counting Stars(三元环计数)
题意:
给一张n个点m条边的无向图,问有多少个\(A-structure\)
其中\(A-structure\)满足\(V=(A,B,C,D)\) && \(E=(AB,BC,CD,DA,AC)\)
显然\(A-structure\)是由两个有公共边的三元环构成的
\(1 <=n <= 1e5\)
\(1 <= m <= min(2e5,n*(n-1)/2)\)
思路:
三元环计数
做法1、
①统计每个点的度数
②入度\(<=sqrt(m)\)的分为第一类,入度\(>sqrt(m)\)的分为第二类
③对于第一类,暴力每个点,然后暴力这个点的任意两条边,再判断这两条边的另一个端点是否连接
因为\(m\)条边最多每条边遍历一次,然后暴力的点的入度\(<=sqrt(m)\),所以复杂度约为\(O(msqrt(m))\)
④对于第二类,直接暴力任意三个点,判断这三个点是否构成环,因为这一类点的个数不会超过\(sqrt(m)\)个,所以复杂度约为\(O(sqrt(m)^3)=O(msqrt(m))\)
⑤判断两个点是否连接可以用set,map和Hash都行,根据具体情况而行
这种做法建的是双向边,常数很大
更优的做法2、建有向边 复杂度为\(O(msqrt(n))\)
对所有边按照两端点的度数定向,度数小的往度数大的走,度数相同则按编号小到大走,这样定向后
可以保证是个有向无环图。
为什么呢,要想定向存在环,则这个环上的点度数必须相同,由于保证了编号从小到大走
所以是不可能存在环的。
这样定向同时还保证了每个点的出度不超过\(sqrt(n)\),乍一眼看觉得应该是\(sqrt(m)\)的,仔细一想确实比这个低,不过不知道怎么证。
对于这道题,我们在求三元环的时候,统计一下每条边有多少对点能构成三元环,\(C(cnt,2)\)累计一下即可
做法1、
#include<bits/stdc++.h>
#define LL long long
using namespace std;
void read(int &x){
x = 0;
char c = getchar();
while(c < '0' || c > '9') c = getchar();
while(c >= '0' && c <= '9') x = x * 10 + c - '0',c = getchar();
}
const int N = 1e5 + 10;
set<LL> g;
int deg[N];
vector<int> G[N];
int vis[N],vi[N];
int main(){
int n, m, u, v, Sz;
while(scanf("%d%d",&n,&m) != EOF){
Sz = sqrt(m + 0.5);
g.clear();
for(int i = 1;i <= n;i++){
vis[i] = vi[i] = deg[i] = 0;
G[i].clear();
}
for(int i = 0;i < m;i++){
scanf("%d%d",&u,&v);
g.insert(u + 1LL * v * n);
g.insert(v + 1LL * u * n);
deg[u]++,deg[v]++;
G[u].push_back(v);
G[v].push_back(u);
}
LL ans = 0;
for(int u = 1;u <= n;u++){
vis[u] = 1;
for(auto v:G[u]) vi[v] = u;
for(auto v:G[u]){
int cnt = 0;
if(vis[v]) continue;
if(deg[v] <= Sz){
for(auto vv:G[v]){
if(vi[vv] == u) cnt++;
}
}else{
for(auto vv:G[u]){
if(g.find(1LL * v * n + vv) != g.end()) cnt++;
}
}
ans += 1LL * cnt * (cnt - 1) / 2;
}
}
printf("%lld\n",ans);
}
return 0;
}
做法二、
#include<bits/stdc++.h>
#define LL long long
#define P pair<int,int>
using namespace std;
void read(int &x){
x = 0;
char c = getchar();
while(c < '0' || c > '9') c = getchar();
while(c >= '0' && c <= '9') x = x * 10 + c - '0',c = getchar();
}
const int N = 1e5 + 10;
set<LL> g;
int deg[N];
vector<pair<int,int> > G[N];
int vi[N];
int X[N * 2],Y[N * 2],cnt[N],pos[N];
int main(){
int n, m, u, v, Sz;
while(scanf("%d%d",&n,&m) != EOF){
Sz = sqrt(m);
for(int i = 1;i <= n;i++){
vi[i] = deg[i] = pos[i] = 0;
G[i].clear();
}
g.clear();
int tot = 0;
for(int i = 0;i < m;i++){
scanf("%d%d",&X[i],&Y[i]);
u = X[i],v = Y[i];
deg[u]++,deg[v]++;
}
for(int i = 0;i < m;i++){
cnt[i] = 0;
if(deg[X[i]] < deg[Y[i]]) G[X[i]].push_back(make_pair(Y[i],i));
else if(deg[Y[i]] < deg[X[i]]) G[Y[i]].push_back(P(X[i],i));
else{
if(X[i] < Y[i]) G[X[i]].push_back(P(Y[i],i));
else G[Y[i]].push_back(P(X[i],i));
}
}
LL ans = 0;
for(int i = 0;i < m;i++){
u = X[i],v = Y[i];
for(auto vp:G[u]) pos[vp.first] = vp.second,vi[vp.first] = i + 1;
for(auto vp:G[v]){
int vv = vp.first;
if(vi[vv] == i + 1){
cnt[i]++;
cnt[pos[vv]]++;
cnt[vp.second]++;
}
}
}
for(int i = 0;i < m;i++) ans += 1LL * cnt[i] * (cnt[i] - 1) / 2;
printf("%lld\n",ans);
}
return 0;
}