[真正的单色三角形] BZOJ 3498 PA2009 Cakes

最新推荐文章于 2021-06-26 10:40:34 发布
最新推荐文章于 2021-06-26 10:40:34 发布
阅读量1.1k 收藏
点赞数
分类专栏： 图论 组合数学 暴力 杂题 Notes
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014609452/article/details/60368425
版权
杂题 同时被 3 个专栏收录
74 篇文章 0 订阅
订阅专栏
图论
66 篇文章 0 订阅
订阅专栏
组合数学
35 篇文章 0 订阅
订阅专栏

Description

给定一个无自环重边的无向图，求这个图的三元环的个数以及补图的三元环个数。
n≤105,m≤105

PS. 可以在这里交

Analysis

这个东西 首先我们可以用BZOJ 2916 [Poi1997]Monochromatic Triangles的方法算出原图三角形和补图三角形的个数之和 然后我们再算出原图三角形的个数就好了
以下的代码都是求原图三角形的个数 复杂度 O(mm−−√) 一般情况下跑得飞快

ps不知道哪里可以交啊写的代码没测过就放上来还说跑得快我也很不要脸啊

Solution1

首先我们可以考虑按阈值 m−−√ 把点分为两类
[奇葩方法]三元环题解

将所有点分成两类:度数 < sqrt(m)的和度数 > sqrt(m)的.
先求包含第一类点的三元环个数.
由于边很少,所以枚举2条边即可.由于一个点的度不超过sqrt(m),所以一条边最多被枚到(sqrt(m))次,最多枚M条边,所以这个操作时O(m*sqrt(m))的.
再求不包含第一类点的三元环个数.
由于每条边贡献2个度,所以二类点的数量是O(sqrt(m))级的.直接枚举三个点,复杂度O((sqrt(m))^3)=O(m*sqrt(m))
所以算法总的复杂度是O(m*sqrt(m))的.

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<cmath>
#define push_back pb
using namespace std;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=200005;

namespace HashMap{
  const int P=1000007,seed=2333;
  int u[N<<1],v[N<<1],next[N<<1];
  int head[P],inum;
  inline void add(int _u,int _v){
    int t=(_u*seed+_v)%P;
    u[++inum]=_u; v[inum]=_v; next[inum]=head[t]; head[t]=inum;
  }
  inline bool query(int _u,int _v){
    int t=(_u*seed+_v)%P;
    for (int p=head[t];p;p=next[p])
      if (u[p]==_u && v[p]==_v)
    return 1;
    return 0;
  }
}

struct edge{
  int u,v,next;
}G[N<<1];
int head[N],inum;
inline void add(int u,int v,int p){
  G[p].u=u; G[p].v=v; G[p].next=head[u]; head[u]=p;
}

int n,m;
int deg[N];
int s[N],pnt;
int Ans;
#define V1 G[p].v
#define V2 G[q].v
int main(){
  int iu,iv;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(n); read(m);
  for (int i=1;i<=m;i++){
    read(iu),read(iv);
    add(iu,iv,++inum); add(iv,iu,++inum);
    deg[iu]++,deg[iv]++;
    HashMap::add(iu,iv),HashMap::add(iv,iu);
  }
  int B=sqrt(m);
  for (int i=1;i<=n;i++)
    if (deg[i]<=B){
      for (int p=head[i];p;p=G[p].next)
    if (!(deg[V1]<=B && V1<i))
      for (int q=G[p].next;q;q=G[q].next)
        if (!(deg[V2]<=B && V2<i))
          if (HashMap::query(G[p].v,G[q].v))
        Ans++;
    }else
      s[++pnt]=i;
  for (int i=1;i<=pnt;i++)
    for (int j=i+1;j<=pnt;j++)
      if (HashMap::query(s[i],s[j]))
    for (int k=j+1;k<=pnt;k++)
      if (HashMap::query(s[j],s[k]) && HashMap::query(s[i],s[k]))
        Ans++;
  printf("%d\n",Ans);
  return 0;
}

Solution2

可以证明似乎三角形的个数就是 O(mm−−√) ？
求无向图三元环的个数，这种做法是O(m*sqrt(m))的吗？

给定n个点m条边的稀疏图，对于每一条边，枚举两个端点中度数较少的端点的邻接点，判断是否构成三元环。这种做法的判断次数在m条边形成一个完全图时达到最大，即O(m*sqrt(m))。不知道这个结论是不是对的？

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<cmath>
#define push_back pb
using namespace std;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=200005;

namespace HashMap{
  const int P=1000007,seed=2333;
  int u[N<<1],v[N<<1],next[N<<1];
  int head[P],inum;
  inline void add(int _u,int _v){
    int t=(_u*seed+_v)%P;
    u[++inum]=_u; v[inum]=_v; next[inum]=head[t]; head[t]=inum;
  }
  inline bool query(int _u,int _v){
    int t=(_u*seed+_v)%P;
    for (int p=head[t];p;p=next[p])
      if (u[p]==_u && v[p]==_v)
    return 1;
    return 0;
  }
}

struct edge{
  int u,v,next;
}G[N<<1];
int head[N],inum;
inline void add(int u,int v,int p){
  G[p].u=u; G[p].v=v; G[p].next=head[u]; head[u]=p;
}

int n,m;
int u[N],v[N];
int deg[N];
int Ans;
#define V G[p].v
int main(){
  freopen("t.in","r",stdin);
  freopen("t2.out","w",stdout);
  read(n); read(m);
  for (int i=1;i<=m;i++){
    read(u[i]),read(v[i]);
    add(u[i],v[i],++inum); add(v[i],u[i],++inum);
    deg[u[i]]++,deg[v[i]]++;
    HashMap::add(u[i],v[i]),HashMap::add(v[i],u[i]);
  }
  for (int i=1;i<=m;i++){
    int u=::u[i],v=::v[i];
    if (deg[u]>deg[v]) swap(u,v);
    for (int p=head[u];p;p=G[p].next)
      if (v!=V && HashMap::query(v,V))
    Ans++;
  }
  printf("%d\n",Ans/3);
  return 0;
}

学习了这里的写法 其实是一样的 不过常数小很多

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<cmath>
#define push_back pb
using namespace std;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=200005;

namespace HashMap{
  const int P=1000007,seed=2333;
  int u[N<<1],v[N<<1],next[N<<1];
  int head[P],inum;
  inline void add(int _u,int _v){
    int t=(_u*seed+_v)%P;
    u[++inum]=_u; v[inum]=_v; next[inum]=head[t]; head[t]=inum;
  }
  inline bool query(int _u,int _v){
    int t=(_u*seed+_v)%P;
    for (int p=head[t];p;p=next[p])
      if (u[p]==_u && v[p]==_v)
    return 1;
    return 0;
  }
}

struct edge{
  int u,v,next;
}G[N<<1];
int head[N],inum;
inline void add(int u,int v,int p){
  G[p].u=u; G[p].v=v; G[p].next=head[u]; head[u]=p;
}

int n,m;
int deg[N];
int Ans;
#define V G[p].v
int s[N],pnt;
int main(){
  int iu,iv;
  freopen("t.in","r",stdin);
  freopen("t3.out","w",stdout);
  read(n); read(m);
  for (int i=1;i<=m;i++){
    read(iu),read(iv);
    add(iu,iv,++inum); add(iv,iu,++inum);
    deg[iu]++,deg[iv]++;
    HashMap::add(iu,iv),HashMap::add(iv,iu);
  }
  for (int u=1;u<=n;u++){
    pnt=0;
    for (int p=head[u];p;p=G[p].next)
      if (deg[u]<deg[V] || (deg[u]==deg[V] && u<V))
    s[++pnt]=V;
    for (int i=1;i<=pnt;i++)
      for (int j=1;j<i;j++)
    if (HashMap::query(s[i],s[j]))
      Ans++;
  }
  printf("%d\n",Ans);
  return 0;
}

UPD.找到题了 跑得并不快 555 好像大家各有奇技淫巧

BZOJ 3498: PA2009 Cakes

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<cmath>
#define push_back pb
using namespace std;
typedef long long ll;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=250005;

namespace HashMap{
  const int P=1000007,seed=2333;
  int u[N<<1],v[N<<1],next[N<<1];
  int head[P],inum;
  inline void add(int _u,int _v){
    int t=(_u*seed+_v)%P;
    u[++inum]=_u; v[inum]=_v; next[inum]=head[t]; head[t]=inum;
  }
  inline bool query(int _u,int _v){
    int t=(_u*seed+_v)%P;
    for (int p=head[t];p;p=next[p])
      if (u[p]==_u && v[p]==_v)
    return 1;
    return 0;
  }
}

struct edge{
  int u,v,next;
}G[N<<1];
int head[N],inum;
inline void add(int u,int v,int p){
  G[p].u=u; G[p].v=v; G[p].next=head[u]; head[u]=p;
}

int n,m;
int deg[N];
ll Ans;
int val[N];
#define V G[p].v
int s[N],pnt;
int main(){
  int iu,iv;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(n); read(m);
  for (int i=1;i<=n;i++) read(val[i]);
  for (int i=1;i<=m;i++){
    read(iu),read(iv);
    add(iu,iv,++inum); add(iv,iu,++inum);
    deg[iu]++,deg[iv]++;
    HashMap::add(iu,iv),HashMap::add(iv,iu);
  }
  for (int u=1;u<=n;u++){
    pnt=0;
    for (int p=head[u];p;p=G[p].next)
      if (deg[u]<deg[V] || (deg[u]==deg[V] && u<V))
    s[++pnt]=V;
    for (int i=1;i<=pnt;i++)
      for (int j=1;j<i;j++)
    if (HashMap::query(s[i],s[j]))
      Ans+=max(val[u],max(val[s[i]],val[s[j]]));
  }
  printf("%lld\n",Ans);
  return 0;
}
里阿奴摩西
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
无向图中三角形个数
yang20141109的专栏
12-22 6350
给定一个无向图，求出此图中包含三角形的个数。 图中三角形指的是：在图中任意选择三个顶点，且这三个顶点之间任意两个点都有一条来自原图中的边。求出图中三角形个数有许多应用(1)在qq中两个好友之间的共同好友个数(2)好友推荐应用。(3)图的聚集系数。(4)小世界网络。 本文中所给出的算法思想： 1.对图中的所有顶点按照度数由小到大进行排序。 2.计算图中每个顶点
简单环（状压 dp，详细注释）
百般回首曾经的相遇，留下的只是一缕残忆
10-06 3476
简单环 题目描述 给定一张n个点m条边的无向图，求出图中所有简单环的数量。(简单环：简单环又称简单回路，图的顶点序列中，除了第一个顶点和最后一个顶点相同外，其余顶点不重复出现的回路叫简单回路。或者说，若通路或回路不重复地包含相同的边，则它是简单的) 输入描述: 第一行三个数n m k (k在输出描述中提到) 接下来m行每行两个数u，v表示u到v之间存在一条边（无重边） 输出描述: 输出k行每行1个...
无向图判断三元环
weixin_33980459的博客
04-10 450
在无向图中，如果三个不同的顶点之间都有边，则称他们组成了一个三角形。 在一张无向图 G 中，有且仅有一个三角形。现在你的任务是找到它。 输入格式 第一行两个数n,mn,m，表示 G 的顶点个数和边的条数。 接下来mm行，每行两个数i,ji,j表示点ii和jj之间有一条边。题目保证没有重边和自环。 输出格式 输出一行，三个整数，i<j<ki<j&...
单色三角形问题（容斥定理）
llin-黎辰
03-22 1289
空间里有n个点，任意三点不共线，每两个点之间都用黑色或红色线段连接，如果一个三角形的三条边同色，那么称这个三角形为单色三角型，给出红色线段的列表 问题的解决：如果直接枚举点的话，复杂度是n^3，虽然也是多项式复杂度，但是适用的数据量太小，但是如果我们换个角度思想，从补集的角度切入，因为三角形的个数是C(n,3),也就是全集已知，对于每个点，只要与它相连的两个边颜色不同，就不能构成单色三角形，而恰
无向图求三元环
codeswarrior的博客
07-14 1536
无向图是一个稀疏图，点n&amp;lt;105n&amp;lt;105n &lt; 10^5 边数m&amp;lt;min(105,n×(n−1)2)m&amp;lt;min(105,n×(n−1)2)m &lt; min(10^5,\frac{n\times(n-1)}{2}),现在给出n,m,E(边集合)，求出这个三元环的个数，两个三元环不同是：当于对两个三元环的边，某个三元环存在一条边在另一个三元环中无法找到。 解法：暴力的方法可...
单色三角形问题
maravilla_evol
04-14 2739
空间里有n个点，任意三点不共线。每两个点之间都用红色或者黑色线段链接。如果一个三角形的三条边同色，责成这个三角形是单色三角形。对于给定的红色线段列表，找出单色三角形的个数。1，如果直接找需要枚举所有的三个点的组合，有C(n,3)种组合，当n比较大的时候，需要枚举的组合很多，复杂度为O(n3)；但是在遍历每个点的时候都可以枚举颜色不同的两条边。则每个点上不同的颜色的三角形有R*B，R为红色边的个数，B为黑色边的个数。所有点上的不同颜色边组之和为Q，则总共有C(n,3)-Q/2种相同颜色的边。时间复杂度变为了O
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
BZOJ3230相似子串测试数据
08-19
标题"BZOJ3230相似子串测试数据"指的是一个与在线判题平台BZOJ（BestCoder ZOJ，中国大学算法竞赛在线评测系统）相关的编程问题，具体问题编号为3230，涉及的主要知识点是"相似子串"的搜索和处理。在编程竞赛中，...
bzoj problem
08-21
【标题】"bzoj problem" 指的是“八中OJ（Online Judge）问题集”，这是一个在线编程竞赛平台的题目合集。八中OJ可能是某所中学的在线编程训练系统，它提供了大量的编程题目供用户练习和挑战，以提升编程能力和算法...
BZOJ题目镜像
09-14
【BZOJ题目镜像】是为了解决BZOJ（Blue Zoo Online Judge）系统因服务器问题时常出现访问不稳定或无法访问的情况而设立的一个备份资源。BZOJ是中国早期的在线编程竞赛平台，提供了大量的算法题目供编程爱好者进行...
统计无向图中三角形的个数，复杂度m*sqrt(m).
weixin_33973600的博客
10-17 1055
统计无向图中三角形的个数，复杂度m*sqrt(m). #include<stdio.h> #include<vector> #include<set> #include<math.h> #include<algorithm> using namespace std; #define LL long...
无向图三元环计数 (思维，图论)
chineseherofeng的博客
06-26 587
题目链接 思路： 看完这题第一眼思路是直接枚举两条邻边(m^2)，但一看数据范围2e5便知道这个办法不可行，于是考虑将图简化。 简化方法：将原图中度数小的点向度数大的点连一条有向边，如果度数相同便从编号小的向编号大的连边，可以发现这样简化完后的新图是一个DAG(因为如果要成还的话环里的边必然有一条是从度数大的连向度数小的，或者从编号大的连向编号小的，与规则矛盾)，而原图中的每个三元环(u,v,w)在新图中一定有一个(u->v,u->w,v->w)的子图与之相对应，我们只需要在新图中枚举每个
求无向图中的三元环个数
FlushHip
04-16 6917
无向图是个稀疏图，点数n<105n<10^5, 边数m<min(105,n∗(n−1)2)m<min(10^5, \frac{n*(n - 1)}{2}).现在给出n,m,E(边的集合)n, m,E(边的集合)。求出这个无向图中三元环的个数，两个三元环不同是：当对于两个三元环的边，某个三元环存在一条边在另一个三元环的边中无法找到！首先想到的方法只能是暴力求解，很直接的暴力就是枚举每一条边，以一条边度
[奇葩方法]三元环题解
qw
10-16 2810
[题目描述] 给你一个包含n个点,m条边的无向图,求三元环个数. [数据范围] n,m [题解]         这个方法是lypYY出的一个奇葩方法.ORZlyp,YY都能YY出这种方法.         将所有点分成两类:度数sqrt(m)的.         先求包含第一类点的三元环个数.         由于边很少,所以枚举2条边即可.由于一个点的度不超过sqrt(m),所
HDU-6184 (无向图三元环计数)
yo_bc的博客
09-06 1867
Little A is an astronomy lover, and he has found 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 m non-directional edges. It is guranteed that no edges connect one star with it
[竞赛图三元环存在性判断] HDU 4324 Triangle LOVE
雯舞
03-04 1853
给一个竞赛图 判断是否有三元环 首先 如果我们任意一个环 就可以找到一个三元环 证明： 假设有环上三个相邻的点a-> b-> c 如果c->a有边 就已经形成了一个三元环 如果c->a没边 那么a->c肯定有边 这样就形成了一个n-1元环 递归即可 我们可以考虑Tarjan求出强连通分量即可一个强连通分量必定包含一个环 不考虑读入 判断单次 O(n+m)O(n
【PA2009】【BZOJ3498】Cakes
CreationAugust is 14 years old forever
02-25 1542
DescriptionN个点m条边，每个点有一个点权a。 对于任意一个三元环(i，j，k)（i< j< k），它的贡献 为max(ai，aj，ak) 求所有三元环的贡献和。 N< 100000，，m< 250000。InputThe first line of the standard input contains two integers n and m (1<=N<=100000,1
HDU 6184 Counting Stars（无向图三元环计数）
yasolx的博客
09-05 1746
Counting StarsTime Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 459 Accepted Submission(s): 115Problem Description Little A is an astronomy l
bzoj.org俄罗斯方块
最新发布
05-06
bzoj.org是一个在线的编程题库，其中包含了很多经典算法问题和实际应用中的编程问题。其中包括了很多有趣的游戏算法问题，比如Puzzle，Tetris等。 而bzoj.org俄罗斯方块就是在该平台上的一个经典的俄罗斯方块游戏。在该游戏中，玩家需要控制方块不断向下移动，直到方块落到底部或者与其他方块重叠。游戏过程中，玩家可以通过旋转和移动方块来调整方块的位置，使其能够合理地落下。当一行方块被填满时，该行会被清除，并获得相应的分数。 bzoj.org俄罗斯方块不仅仅是一款游戏，更是一个很好的算法练习平台，你可以通过该游戏来学习和研究游戏算法问题。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值