当前位置:   article > 正文

黑科技——三元环计数

三元环计数

如果给你一个无向图,问其中有多少个三元环

暴力?一个菊花图就没了

我们考虑如果一条边(u,v)由度数大的向度数小的连边。如果度数一样,那么编号小的向编号大的连边(这个编号大的向小的连也可以)。

这时候整张图是个DAG,所以只要枚举每个点,给每个点的出点打上标记,再枚举出点的出点,如果有标记就计数,显然每个三元环只会被算到一次。

我们考虑一下时间复杂度

\sum_{i=1}^{m}degout(u_{i})

乍一看时间复杂度不太对劲,然而,由于我们每次是由度数大的向度数小的连边,所以:

1.如果degout(u)>\sqrt{m},这样的u只会被算根号次,因为只有度数大于等于根号m的点才能连向他,而这样的点最多只有根号m个

2.degout(u)<=\sqrt{m}

综上时间复杂度为m\sqrt{m}

例题:

hdu6184

  1. #include<bits/stdc++.h>
  2. #define ll long long
  3. using namespace std;
  4. const int N=2e5+10;
  5. struct node{
  6. int u,to,w;
  7. };
  8. node edge[N<<1];
  9. int head[N],s[N],t[N],x[N],y[N],deg[N];
  10. int k,n,m;
  11. ll ans,la[N];
  12. void add(int x,int y,int z){
  13. edge[k].u=y; edge[k].to=head[x]; edge[k].w=z; head[x]=k++;
  14. }
  15. int main(){
  16. while (~scanf("%d%d",&n,&m)) {
  17. for (int i=1;i<=n;i++) head[i]=-1; k=0;
  18. ans=0;
  19. for (int i=1;i<=n;i++) deg[i]=0;
  20. for (int i=1;i<=m;i++) {
  21. scanf("%d%d",&x[i],&y[i]); deg[x[i]]++; deg[y[i]]++;
  22. la[i]=0;
  23. }
  24. for (int i=1;i<=n;i++) s[i]=0,t[i]=0;
  25. for (int i=1;i<=m;i++) if (deg[x[i]]>deg[y[i]]||((deg[x[i]]==deg[y[i]])&&x[i]<y[i]))
  26. add(x[i],y[i],i); else add(y[i],x[i],i);
  27. for (int i=1;i<=n;i++) {
  28. for (int j=head[i];j!=-1;j=edge[j].to){
  29. int u=edge[j].u; s[u]=i; t[u]=edge[j].w;
  30. }
  31. for (int j=head[i];j!=-1;j=edge[j].to){
  32. int u=edge[j].u;
  33. for (int p=head[u];p!=-1;p=edge[p].to){
  34. int v=edge[p].u;
  35. if (s[v]==i) {
  36. ans+=la[edge[j].w]+la[edge[p].w]+la[t[v]];
  37. la[edge[j].w]++; la[edge[p].w]++; la[t[v]]++;
  38. }
  39. }
  40. }
  41. }
  42. printf("%lld\n",ans);
  43. }
  44. }

 

声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:【wpsshop博客】
推荐阅读
相关标签
  

闽ICP备14008679号