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

最新推荐文章于 2019-04-23 19:36:28 发布
最新推荐文章于 2019-04-23 19:36:28 发布
阅读量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;
}
里阿奴摩西
关注
专栏目录
【2-SAT+前缀优化建图】BZOJ3495 PA2010 Riddle
linkfqy
07-28 1182
题面在这里很典型的2-SAT问题…… 可以对每个点是否为首都建图 考虑一条边(u,v)，如果u不是首都，则v必须是。反之亦然。 u′→v，v′→uu'\rightarrow v，v'\rightarrow u 然后对于一个国家里的点，有一个是首都则其他都不是首都但是这样边数是N2N^2级别的可以用前缀优化建图。 对于一个点u，建2个变量 分别表示 u是否为首都（u与u’） 和 u所在的国
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
【PA2009】【BZOJ3498】Cakes
CreationAugust is 14 years old forever
02-25 1580
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
BZOJ.3498.[PA2009]Cakes(三元环 枚举)
weixin_30341735的博客
07-13 161
题目链接 感觉我可能学的假的（复杂度没问题，但是常数巨大）。 一个比较真的说明见这儿：https://czyhe.me/blog/algorithm/3-mem-ring/3-mem-ring/。 \(Description\) n个点m条边的无向图，每个点有点权。对于任意一个三元环\((i,j,k),i<j<k\)，其贡献为\(max\{a_i,a_j,a_k\}\)。求所有三元环...
[BZOJ3498]-[PA2009]Cakes-阈值法
泉華子的OI足迹
05-23 356
说在前面 这个题真的是，卡的me心累 复杂度只要劣一点都过不去，T了一页最后学着别人的写法过了 题目 BZOJ3498传送门 题目大意 给出一张 nnn 个点 mmm 条边的图，点有点权 现在对于图中的每个三元环，定义其价值为 max(ai,aj,ak)max(ai,aj,ak)\max(a_i,a_j,a_k) 求出所有三元环的总价值 输入输出格式 输入格式：...
bzoj 3498 PA2009 Cakes
aoping9329的博客
08-17 108
N个点m条边，每个点有一个点权a。 对于任意一个三元环(i，j，k)（i<j<k），它的贡献为max(ai，aj，ak) 求所有三元环的贡献和。 N<100000 m<250000 模拟题意即可…… 1 #include <bits/stdc++.h> 2 using namespace std; 3 const...
bzoj 3498
weixin_30716725的博客
08-25 87
统计三元环 很多代码在bzoj都T诶 #include <iostream> #include <cstdio> #include <algorithm> #include <cmath> #include <cstring> #include <string> #define gc getchar()...
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
BZOJ题目镜像
09-14
【BZOJ题目镜像】是为了解决BZOJ（Blue Zoo Online Judge）系统因服务器问题时常出现访问不稳定或无法访问的情况而设立的一个备份资源。BZOJ是中国早期的在线编程竞赛平台，提供了大量的算法题目供编程爱好者进行...
bzoj problem
08-21
【标题】"bzoj problem" 指的是“八中OJ（Online Judge）问题集”，这是一个在线编程竞赛平台的题目合集。八中OJ可能是某所中学的在线编程训练系统，它提供了大量的编程题目供用户练习和挑战，以提升编程能力和算法...
【bzoj3498】PA2009 Cakes 判断三元环
wang_sj的博客
01-13 548
Description N个点m条边，每个点有一个点权a。 对于任意一个三元环(j，j，k)（i 为max(ai，aj，ak) 求所有三元环的贡献和。 N Input The first line of the standard input contains two integers n and m (1Bi) separated by a single space. They d
找三元环（bzoj 3498: PA2009 Cakes）
加载中...
08-29 1343
问题：给你一个无向图，问其中有多少个三元环 方法：暴力 复杂度：O(msqrt(m)) 步骤：按点的度数分成两类，分别暴力 ①统计每个点的度数 ②入度sqrt(m)的分为第二类 ③对于第一类，暴力每个点，然后暴力这个点的任意两条边，再判断这两条边的另一个端点是否连接 因为m条边最多每条边遍历一次，然后暴力的点的入度 ④对于第二类，直接暴力任意三个点，判断这三个点是否构成环，因为这一
[bzoj3498]Cakes
不来也不去的一只失忆蝴蝶
05-04 865
题目描述N个点m条边，每个点有一个点权a。 对于任意一个三元环(i，j，k)（i<j<k）(i，j，k)（i<j<k），它的贡献 为max(ai，aj，ak) 求所有三元环的贡献和。 N<100000，，m<250000。瞎做我们都知道一个定理就是一个图，枚举不比一个点度数小的点只有根号m个。 因此我们可以以度数来定义优先级，度数大优先级高，相同度数编号越大优先级越高。 对于一个三元环
单色三角形
maravilla_evol
04-14 2149
空间里有n个点，任意三点不共线。每两个点之间都用红色或者黑色线段链接。如果一个三角形的三条边同色，责成这个三角形是单色三角形。对于给定的红色线段列表，找出单色三角形的个数。分析：如果直接找需要枚举所有的三个点的组合，有C(n,3)种组合，当n比较大的时候，需要枚举的组合很多，复杂度为O(n3)；但是在遍历每个点的时候都可以枚举颜色不同的两条边。则每个点上不同的颜色的三角形有R*B，R为红色边的个数，B为黑色边的个数。所有点上的不同颜色边组之和为Q，则总共有C(n,3)-Q/2种相同颜色的边。时间复杂度变为了
单色三角形，同色三角形
Light
04-23 815
平面上有n个点（n≤8000），每两个点之间都有一条红色或者是黑色的线段，任意三点均不共线。 现在，已知哪些点之间连的线段是红色的，剩下的线段都是黑色的，要求计算这些点组成的三角 形中有多少是同色的（顶点编号从1到n）？ 输入描述 第一行是n, m（3≤n≤8000），n表示点的个数，m表示红色线段的条数。下面m行，每 一行都是两个整数a和b，表示点a和点b之间的线段是红色的（a<b）。 当组成异色三角形时， 必定有两个顶点的相邻两边为异色， 故可以算出所有异色三角形的个数
Bzoj 3498 Cakes（三元环）
weixin_30243533的博客
12-15 92
题面（权限题就不放题面了） 题解 三元环模板题，按题意模拟即可。 #include <cstdio> #include <cstring> #include <vector> using std::vector; const int N = 1e5 + 10, M = 2.5e5 + 10; int n, m, a[N], deg[N], u[M], v[M]...
BZOJ3498: PA2009 Cakes(三元环)
weixin_34203832的博客
10-25 101
题意 题目链接 Sol 按照套路把边转成无向图，我们采取的策略是从权值大的向权值小的连边 然后从按权值从小到大枚举每个点，再枚举他们连出去的点\(v\) 如果\(v\)的度数\(\leqslant M\)，那么就再暴力枚举\(v\)连出去的点\(t\)，看\(u\)与\(t\)是否联通(打标记) 否则暴力枚举\(u\)连出去的点\(t\)，看\(v\)与\(t\)是否联通(直接hash表) 复杂度...
《大白书》105页 单色三角形（组合计数）
细雨欣然
03-21 240
时间限制：1秒 内存限制：64M 【问题描述】 　　给定空间里的n个点，其中没有3个点共线。每两个点之间都用红色或黑色线段连接。求3条边同色的三角形个数。 【输入格式】 　　第一行一个整数n（n<=1000），表示点数。 　　接下来的n行表示一个n*n的矩阵，第i行第j列的数字表示连接点i与点j的线段的颜色，1表示红色，2表色黑色。当然i=j时为0。输入中i行j列与j行i列一定形同。【输出格
单色三角形问题
maravilla_evol
04-14 2786
空间里有n个点，任意三点不共线。每两个点之间都用红色或者黑色线段链接。如果一个三角形的三条边同色，责成这个三角形是单色三角形。对于给定的红色线段列表，找出单色三角形的个数。1，如果直接找需要枚举所有的三个点的组合，有C(n,3)种组合，当n比较大的时候，需要枚举的组合很多，复杂度为O(n3)；但是在遍历每个点的时候都可以枚举颜色不同的两条边。则每个点上不同的颜色的三角形有R*B，R为红色边的个数，B为黑色边的个数。所有点上的不同颜色边组之和为Q，则总共有C(n,3)-Q/2种相同颜色的边。时间复杂度变为了O
bzoj.org俄罗斯方块
最新发布
05-06
bzoj.org是一个在线的编程题库，其中包含了很多经典算法问题和实际应用中的编程问题。其中包括了很多有趣的游戏算法问题，比如Puzzle，Tetris等。 而bzoj.org俄罗斯方块就是在该平台上的一个经典的俄罗斯方块游戏。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值