全局最小割 Stoer-Wagner (SW算法)优化

优化吃藕了,感谢放宽时限,感谢平板电视 (pb_ds)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <ext/pb_ds/priority_queue.hpp> typedef long long ll;
typedef __gnu_pbds::priority_queue<std::pair<int, int>, std::less<std::pair<int, int> >, __gnu_pbds::binomial_heap_tag >pq; struct Edge {int from, to, nxt, val;}; const int maxn = 3e3 + ;
const int maxm = 4e5 + ;
const int inf = ;
int n, r;
bool vis[maxn], bin[maxn];
Edge e[maxm * + ];
int head[maxn], id[maxn][maxn], esz; inline ll read() {
int x = , f = ; char ch = getchar();
while(ch < '' || ch > '') {if(ch == '-')f = -; ch = getchar();}
while(ch >= '' && ch <= '') {x = x * + ch - ''; ch = getchar();}
return x * f;
} inline void add_edge(int from, int to, int val) {
if(id[from][to] == -) {
e[esz].from = from, e[esz].to = to, e[esz].val = val, e[esz].nxt = head[from];
id[from][to] = esz; head[from] = esz++;
} else {
e[id[from][to]].val += val;
}
} inline void init() {
memset(bin, false, sizeof(bin));
memset(head, -, sizeof(head));
memset(id, -, sizeof(id));
esz = ;
} pq que;
pq::point_iterator it[maxn];
inline int contract(int &s, int &t) { // Find s,t
memset(vis, false, sizeof(vis));
int i, j, k, mincut, maxc;
while(!que.empty()) que.pop();
for(i = ; i <= n; ++i)
if(!bin[i])it[i] = que.push(std::make_pair(, i));
else it[i] = nullptr;
for(i = ; i <= n; i++) {
k = -; maxc = -;
if(que.empty()) return mincut;
k = que.top().second, maxc = que.top().first; que.pop();
s = t; t = k; mincut = maxc; vis[k] = true;
for(j = head[k]; ~j; j = e[j].nxt) {
int v = e[j].to, w = e[j].val;
if(!bin[v] && !vis[v]) {
que.modify(it[v], std::make_pair((it[v]->first) + w, v));
}
}
}
return mincut;
} inline int Stoer_Wagner() {
int mincut, i, j, s, t, ans;
for(mincut = inf, i = ; i < n; i++) {
ans = contract(s, t);
bin[t] = true;
if(mincut > ans) mincut = ans;
if(mincut == )return ;
for(j = ; j <= n; j++) if(!bin[j]) {
if(id[j][t] == -)continue;
if(id[s][j] == -)add_edge(s, j, e[id[j][t]].val);
else e[id[s][j]].val += e[id[j][t]].val;
if(id[j][s] == -)add_edge(j, s, e[id[j][t]].val);
else e[id[j][s]].val += e[id[j][t]].val;
}
}
return mincut;
} int main() {
int m;
while(scanf("%d%d", &n, &m) != -) {
init();
for(int i = ; i <= m; i++) {
int a, b, c;
a = read(), b = read(), c = read();
add_edge(a, b, c); add_edge(b, a, c);
}
printf("%d\n", Stoer_Wagner());
}
return ;
}

2017"百度之星"程序设计大赛 - 资格赛 1002 度度熊的王国战略的更多相关文章

  1. hdu 6082 度度熊与邪恶大魔王(2017"百度之星"程序设计大赛 - 资格赛 )

    度度熊与邪恶大魔王 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total S ...

  2. 2017"百度之星"程序设计大赛 - 资格赛【1001 Floyd求最小环 1002 歪解(并查集),1003 完全背包 1004 01背包 1005 打表找规律+卡特兰数】

    度度熊保护村庄 Accepts: 13 Submissions: 488 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/3276 ...

  3. [SinGuLaRiTy] 2017 百度之星程序设计大赛-资格赛

    [SinGuLaRiTy-1034] Copyright (c) SinGuLaRiTy 2017. All Rights Reserved. 度度熊保护村庄  Time Limit: 2000/10 ...

  4. 2017"百度之星"程序设计大赛 - 资格赛

    度度熊与邪恶大魔王  Accepts: 3666  Submissions: 22474  Time Limit: 2000/1000 MS (Java/Others)  Memory Limit: ...

  5. 2014年百度之星程序设计大赛 - 资格赛 1002 Disk Schedule(双调欧几里得旅行商问题)

    Problem Description 有非常多从磁盘读取数据的需求,包含顺序读取.随机读取.为了提高效率,须要人为安排磁盘读取.然而,在现实中,这样的做法非常复杂.我们考虑一个相对简单的场景.磁盘有 ...

  6. 2017"百度之星"程序设计大赛 - 资格赛-度度熊与邪恶大魔王(dp+后缀最小值)

    度度熊与邪恶大魔王 思路:由于防御和血量的范围很小,所以暴力枚举出对于每种防御造成的每种伤害所需的最小花费,最后只需在伤害大于等于血量的情况下再找到最小花费(这个只需要后缀最小值预处理一下就可以了) ...

  7. 2018"百度之星"程序设计大赛 - 资格赛 1002 子串查询

    题面又是万能的毒毒熊... 实在不想写了,就只写了这题 记26个前缀和查询枚举最小值直接算 实在是氵的死 而且我忘记输出Case #%d 想了很久 >_< #include<bits ...

  8. 2017"百度之星"程序设计大赛 - 资格赛 寻找母串

    Problem Description 对于一个串S,当它同时满足如下条件时,它就是一个01偏串: 1.只由0和1两种符组成: 2.在S的每一个前缀中,0的个数不超过1的个数: 3.S中0的个数和1的 ...

  9. 2017"百度之星"程序设计大赛 - 资格赛 度度熊的王国战略

    度度熊的王国战略 度度熊国王率领着喵哈哈族的勇士,准备进攻哗啦啦族. 哗啦啦族是一个强悍的民族,里面有充满智慧的谋士,拥有无穷力量的战士. 所以这一场战争,将会十分艰难. 为了更好的进攻哗啦啦族,度度 ...

随机推荐

  1. 蓝牙4.0模块,AT指令集

    一,LED状态 二,蓝牙模块有两种通信模式 1,AT指令模式 2,数据透传模式 三.AT指令程序设计 1.设置模块的名字 void usart3_send_str(char *pbuf) { whil ...

  2. 7月清北学堂培训 Day 2

    今天是林永迪老师的讲授~ 继续昨日的贪心内容. 我们继续看例题: 分析样例的过河方法: 首先1和2先过河,总时间为2: 然后1回来,总时间为3: 然后5和10过河,总时间为13: 然后2回来,总时间为 ...

  3. Spring Boot教程(七)通过springboot 去创建和提交一个表单

    创建工程 涉及了 web,加上spring-boot-starter-web和spring-boot-starter-thymeleaf的起步依赖. <dependencies> < ...

  4. vue 路由跳转记住当前页面位置

    从列表页面跳去详情页面, 在列表页面的生命周期:deactivated  中把当前的scrollTop位置存下来,可以存在localstorage中,也可以存在vuex中, 从详情页面返回列表页面:a ...

  5. OpenDayLight Beryllium版本 下发流表实现hardtimeout

    1. 实验拓扑 2. 创建拓扑 from mininet.topo import Topo class MyTopo(Topo): def __init__(self): # initilaize t ...

  6. Linux系统下查找最近修改过的文件

    Linux的终端上,没有windows的搜索那样好用的图形界面工具,但find命令确是很强大的. 比如按名字查找一个文件,可以用 find / -name targetfilename . 唉,如果只 ...

  7. kotlin 类的委托

    fun main(arg: Array<String>) { val baseImpl = baseImpl() demo(baseImpl).printL() } interface b ...

  8. MySQL标准化、自动化部署

    机器标准化:服务器的硬件配置(CPU.硬盘.内存.RAID)要保持一致 目录标准化: mysql的安装目录,存放数据的目录.日志文件等,多台机器上配置要统一 参数标准化: 多台服务器上,存放的my.c ...

  9. 阶段5 3.微服务项目【学成在线】_day05 消息中间件RabbitMQ_7.RabbitMQ研究-工作模式-工作队列模式

    RabbitMQ有以下几种工作模式 : 1.Work queues 2.Publish/Subscribe 3.Routing 4.Topics 5.Header 6.RPC 1.Work queue ...

  10. double,float,BigDecimal类型数值的操作

    float四舍五入保留两位小数 /** * float四舍五入保留两位小数 * */ public static float formatDecimal(float n) { return (Math ...