[最大密度子图] 2017 计蒜之道 初赛 第三场 腾讯狼人杀

最新推荐文章于 2021-11-04 15:32:01 发布
最新推荐文章于 2021-11-04 15:32:01 发布
阅读量565 收藏
点赞数
分类专栏： 网络流
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014609452/article/details/72863722
版权

这里写图片描述

最大密度子图的建图 详见论文
这个东西是支持强制选的

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define cl(x) memset(x,0,sizeof(x))
using namespace std;
typedef double ld;

inline char nc(){
  static char buf[1000000],*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 ld eps=1e-7;

namespace D{
  const int N=505;
  const int M=N*N*2;
  struct edge{
    int u,v,next; ld f;
  }G[M<<1];
  int head[N],inum=1;
  inline void add(int u,int v,ld f,int p){
    G[p].u=u; G[p].v=v; G[p].f=f; G[p].next=head[u]; head[u]=p;
  }
  inline void link(int u,int v,ld f){
    add(u,v,f,++inum); add(v,u,0,++inum);
  }
#define V G[p].v
  int S,T;
  int Q[N],l,r,dis[N];
  inline bool bfs(){
    for (int i=1;i<=T;i++) dis[i]=-1;
    l=r=-1; Q[++r]=S; dis[S]=0;
    while (l<r){
      int u=Q[++l];
      for (int p=head[u];p;p=G[p].next)
    if (G[p].f>eps && dis[V]==-1){
      dis[V]=dis[u]+1; Q[++r]=V;
      if (V==T) return 1;
    }
    }
    return 0;
  }
  int cur[N];
  inline ld dfs(int u,ld flow){
    if (u==T) return flow;
    ld used=0;
    for (int p=cur[u];p;p=G[p].next){
      cur[u]=p;
      if (G[p].f>eps && dis[V]==dis[u]+1){
    ld d=dfs(V,min(G[p].f,flow-used));
    G[p].f-=d; G[p^1].f+=d;
    used+=d; if (fabs(flow-used)<eps) break;
      }
    }
    if (used<eps) dis[u]=-1;
    return used;
  }
  inline ld Dinic(){
    ld ans=0;
    while (bfs()){
      memcpy(cur,head,sizeof(int)*(T+5));
      ans+=dfs(S,1e10);
    }
    return ans;
  }
}

const int N=405;

int n,m;
int w[N][N]; ld d[N];
int a[N];

const int U=400*(100+2000);

inline bool check(ld g){
  using namespace D; cl(head); inum=1;
  S=n+1,T=n+2;
  for (int i=1;i<=n;i++){
    d[i]=0;
    for (int j=1;j<=n;j++)
      if (i^j)
    d[i]+=w[i][j]+2*g;
  }
  ld ret=0;
  for (int i=1;i<=n;i++)
    if (!a[i]){
      link(S,i,U);
      link(i,T,U+2.0*(2*n-1)*g-d[i]);
    }else{
      link(S,i,1e10);
      ret+=U+2.0*(2*n-1)*g-d[i];
    }
  for (int i=1;i<=n;i++)
    for (int j=i+1;j<=n;j++)
      link(i,j,w[i][j]+2*g),link(j,i,w[i][j]+2*g);
  ret+=Dinic();
  ret-=U*n;
  ret=-ret;
  return ret>=-eps;
}

int main(){
  int iu,iv,iw;
  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),read(iw),w[iu][iv]=w[iv][iu]=iw;
  for (int i=1;i<=n;i++) read(a[i]);
  ld L=0,R=400,MID;
  while (R-L>eps)
    if (check(MID=(L+R)/2))
      L=MID;
    else
      R=MID;
  printf("%.6lf\n",(L+R)/2);
  return 0;
}

一开始写的最大权闭合子图被卡了

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define cl(x) memset(x,0,sizeof(x))
using namespace std;
typedef double ld;

inline char nc(){
  static char buf[1000000],*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 ld eps=1e-7;

namespace D{
  const int N=100005;
  const int M=N<<2;
  struct edge{
    int u,v,next; ld f;
  }G[M<<1];
  int head[N],inum=1;
  inline void add(int u,int v,ld f,int p){
    G[p].u=u; G[p].v=v; G[p].f=f; G[p].next=head[u]; head[u]=p;
  }
  inline void link(int u,int v,ld f){
    add(u,v,f,++inum); add(v,u,0,++inum);
  }
#define V G[p].v
  int S,T;
  int Q[N],l,r,dis[N];
  inline bool bfs(){
    for (int i=1;i<=T;i++) dis[i]=-1;
    l=r=-1; Q[++r]=S; dis[S]=0;
    while (l<r){
      int u=Q[++l];
      for (int p=head[u];p;p=G[p].next)
    if (G[p].f>eps && dis[V]==-1){
      dis[V]=dis[u]+1; Q[++r]=V;
      if (V==T) return 1;
    }
    }
    return 0;
  }
  int cur[N];
  inline ld dfs(int u,ld flow){
    if (u==T) return flow;
    ld used=0;
    for (int p=cur[u];p;p=G[p].next){
      cur[u]=p;
      if (G[p].f>eps && dis[V]==dis[u]+1){
    ld d=dfs(V,min(G[p].f,flow-used));
    G[p].f-=d; G[p^1].f+=d;
    used+=d; if (fabs(flow-used)<eps) break;
      }
    }
    if (used<eps) dis[u]=-1;
    return used;
  }
  inline ld Dinic(){
    ld ans=0;
    while (bfs()){
      memcpy(cur,head,sizeof(int)*(T+5));
      ans+=dfs(S,1e10);
    }
    return ans;
  }
}

const int N=405;

int n,m;
int w[N][N];
int a[N];

inline bool check(ld ans){
  using namespace D; cl(head); inum=1;
  S=n+n*(n-1)/2+1; T=n+n*(n-1)/2+2;
  int cur=0; ld ret=0;
  for (int i=1;i<=n;i++)
    for (int j=i+1;j<=n;j++){
      ++cur;
      link(S,n+cur,w[i][j]+2*ans);
      if (!a[i]) link(n+cur,i,1e8);
      if (!a[j]) link(n+cur,j,1e8);
      ret+=w[i][j]+2*ans;
    }
  int c=0;
  for (int i=1;i<=n;i++)
    if (!a[i])
      link(i,T,ans*(2*n-1));
    else
      ret+=-ans*(2*n-1);
  ret-=Dinic();
  return ret>=-eps;
}

int main(){
  int iu,iv,iw;
  freopen("t.in","r",stdin);
  freopen("t1.out","w",stdout);
  read(n); read(m);
  for (int i=1;i<=m;i++)
    read(iu),read(iv),read(iw),w[iu][iv]=w[iv][iu]=iw;
  for (int i=1;i<=n;i++) read(a[i]);
  ld L=0,R=1e3,MID;
  while (R-L>eps)
    if (check(MID=(L+R)/2))
      L=MID;
    else
      R=MID;
  printf("%.6lf\n",(L+R)/2);
  return 0;
}
里阿奴摩西
关注
专栏目录
poj 3155 (最大密度子图)
菜鸟起航的专栏
09-27 2081
题意：一个公司有n个人，给出了一些有冲突的人的对数(u,v)，公司决定裁人，那么总裁现在要裁掉冲突率最高的那些人（冲突率=在这些人中存在的冲突数/人数）。就是求出一些点，这些点之间的边数/点数最大。最大密度子图。 思路：胡伯涛的论文《最小割模型在信息学竞赛中的应用》介绍了两种方法： 第一种：转换为最大权闭合图的模型来求解： 设max g = f(x)= |E‘|/|V’| ，找一个子图的
最大密度子图
harry12138小白的博客
02-24 2162
最大密度子图 首先要了解什么是最大密度子图，顾名思义，就是所有子图中密度最大的那一个。那么什么是一个子图的密度呢，这里规定子图的密度就是子图中边数和点数的比值。
[最大密度子图] 2017 计蒜之道 初赛 第三场. 腾讯狼人杀
CE玩家
12-19 435
边带权的最大密度子图带权的最大密度子图有两种做法。1. 二分答案 xx，然后判断 ∑wi,j−xk(2n−k)\sum w_{i,j} -xk(2n-k) 的最大值是否大于零 ∑wi,j−xk(2n−k)\sum w_{i,j} -xk(2n-k) 这东西可以拆成 ∑(wi,j+2x)−k(2n−1)x\sum (w_{i,j}+2x)-k(2n-1)x 相当于每条边权加上 2x2x，点权设为
2017 计蒜之道 初赛 第三场 腾讯狼人杀（简单)(暴搜)
twh233的算法blog
06-20 425
传送门 不久前，腾讯推出了一款改进后的狼人杀游戏。如下是这款游戏的背景和详细规则： 很久很久以前，在莱茵河畔一座岸壁陡峭的山顶上，有一个名叫“杜斯特瓦德”的小村庄。不久前，这个小村庄每晚都会受到狼人的侵袭。每个夜晚，狼人都会在村中进行抢劫，并且会有一个村民成为这群狼人的牺牲品。然而村民们不会坐以待毙，他们试图在白天找到这些狼人并处决他们。 作为杜斯特瓦德村的长老，你需要组织村民来对抗狼人，村
poj 3155(Hard Life)分数规划/最大密度子图
08-20 863
题目链接：http://poj.org/problem?id=3155 题意简述：一个公司有n个人，给出了一些有冲突的人的对数(u,v)，所有为了公司更好的发展，公司的总裁决定裁人，那么总裁现在裁要裁掉冲突率最高的哪些（冲突率=人数/在这些人中存在的冲突数） （题意有所
2017计蒜之道初赛第四场-商汤科技的行人检测（简单）
algzjh的博客
05-28 999
商汤科技的行人检测（简单） 商汤科技近日推出的 SenseVideo 能够对视频监控中的对象进行识别与分析，包括行人检测等。在行人检测问题中，最重要的就是对行人移动的检测。由于往往是在视频监控数据中检测行人，我们将图像上的行人抽象为二维平面上若干个的点。那么，行人的移动就相当于二维平面上的变换。 在这道题中，我们将行人的移动过程抽象为 平移，有两个 移动参数：d_xd​x​​ 和 d_yd​y
腾讯2013编程马拉松初赛第三场
03-22
【腾讯2013编程马拉松初赛第三场】是一场由腾讯主办的编程竞赛，其中包含两道题目，难度相对较高，因为据描述只有一个参赛者完成了全部题目。比赛的标签为“腾讯马拉松”，通常这类比赛旨在挑战程序员的编程技能、...
2017年第二十三届NOIP（C语言）普及组初赛试题及详细答案
03-12
第三轮，他向前走3单位距离，然后右转……他一直这么走下去。请问第2017轮后，他的坐标是：（___,____）。** - **解析:** 该问题可以通过观察规律解决。每四轮，这个人会回到原点，但是离原点更远了。2017除以4的...
【POJ3155】【网络流】【最大密度子图】Hard Life 题解
It's Maverick.
06-07 633
John is a Chief Executive Officer at a privately owned medium size company. The owner of the company has decided to make his son Scott a manager in the company. John fears that the owner will ultimatel
2017 计蒜之道 初赛 第三场 A+B
妍恐龙的博客....^~^眼熟我
05-27 592
在腾讯课堂的物理课上，进行了一个有趣的物理实验。在一个长度为 LL 米的光滑轨道上，小车 A 在 00 时刻以 1\mathrm{m/s}1m/s 的速度从左端出发向右运动，小车 B 在 tt 时刻以 1\mathrm{m/s}1m/s 的速度从右端出发向左运动，两个小车的质量相等。假设所有碰撞都是弹性碰撞，也就是当两个小车相向碰撞时，他们各自会以原来的速度向相反的方向运动；小车和轨道两端发生碰撞时
【学习笔记】最大密度子图
七代目鸽影
11-04 944
内容挺单薄的……
POJ 3155 Hard Life（最大密度子图）
yasolx的博客
08-11 501
Hard Life Time Limit: 8000MS Memory Limit: 65536K Total Submissions: 9111 Accepted: 2648 Case Time Limit: 2000MS Special Judge DescriptionJohn is a Chief Executive Officer at a private
UvaLive 3709 Hard Life(最大密度子图)
静静地码
09-03 711
题目大意：求出最大的密度子图解题思路：详见胡伯涛:算法合集之《最小割模型在信息学竞赛中的应用》第一种解法，最大权闭合的解法#include <cstdio> #include <cstring> #include <algorithm> #include <vector> #include <queue> using namespace std; #define N 1110 const doubl
计蒜客-腾讯狼人杀（网络流+最大密度子图）
十分残念的博客
07-14 570
传送门：https://nanti.jisuanke.com/t/15549 一道论文题，具体求解过程就不详将了（毕竟本人太菜了） 推荐看《最小割模型在信息学竞赛中的应用》这本论文，很好的诠释了最大密度子图在网络流中的应用 题解：网络流+最大密度子图 c(s,v)=U c(v,t)=U+4ng-2g-du c(u,v)=mp[u][v]+2g 其中du为与点u相连的边的总权值，g为二
UVALive 7037 (最大密度子图 网络流)
morejarphone~
03-22 779
每个数都看成一个节点,每个逆序对之间的节点连边,于是只需要求出这个图的最大密度子图,可以用最小割模型解决. #include #include #include #include #include #include #include using namespace std; const double INF = 1e8; #define maxn 1511 #define maxm
NOIP2017初赛普及组C++试题解析
"NOIP2017初赛普及组C++试题" 全国青少年信息学奥林匹克联赛（NOIP）是中国计算机学会举办的一项旨在普及计算机科学知识，提高青少年计算机程序设计能力的比赛。2017年的初赛普及组C++试题主要测试参赛者对于计算机...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值