建图一共建四层 第一层为N个歌手 第二层为{pi,li} 第三层为{si,li} 第四层为M首歌

除了S和第一层与第三层与T之间的边为[1,INF] 其他边均为[0,1]

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef int JQK;
const JQK INF = ;
struct node {
int p, s, l;
bool operator <(const node &rhs) const {
if (p != rhs.p)
return p > rhs.p;
else {
if (s != rhs.s)
return s > rhs.s;
else
return l > rhs.l;
}
}
} nodecnt;
int n, m, k, s1, t1, s2, t2;
map<pair<int, int>, int> pl, sl;
map<int, pair<int, int> > plb, slb;
map<node, int> psl;
vector<int> anser;
namespace dinic {
const int MAXN = ;
const int MAXM = ;
const JQK INF = ;
int Head[MAXN], cur[MAXN], lev[MAXN], to[MAXM << ], nxt[MAXM << ], ed = ;
int S, T, SS, TT, MAXP;
JQK f[MAXM << ], totflow[MAXN];
inline void addedge(int u, int v, JQK cap) {
to[++ed] = v;
nxt[ed] = Head[u];
Head[u] = ed;
f[ed] = cap;
to[++ed] = u;
nxt[ed] = Head[v];
Head[v] = ed;
f[ed] = ;
return;
}
inline void add(int u, int v, JQK caplow, JQK captop) {
addedge(u, v, captop - caplow);
totflow[u] -= caplow, totflow[v] += caplow;
return ;
}
inline bool BFS() {
int u;
for (int i = ; i <= MAXP + ; i++) {
lev[i] = -;
}
//memset(lev, -1, sizeof(lev));
queue<int>q;
lev[S] = ;
q.push(S);
while (q.size()) {
u = q.front();
q.pop();
for (int i = Head[u]; i; i = nxt[i])
if (f[i] && lev[to[i]] == -) {
lev[to[i]] = lev[u] + ;
q.push(to[i]);
/*
if (to[i] == T)
{
return 1;
}
magic one way optimize
*/
}
}
for (int i = ; i <= MAXP + ; i++) {
cur[i] = Head[i];
}
//memcpy(cur, Head, sizeof Head);
return lev[T] != -;
}
inline JQK DFS(int u, JQK maxf) {
if (u == T || !maxf) {
return maxf;
}
JQK cnt = , tem;
for (int &i = cur[u]; i; i = nxt[i])
if (f[i] && lev[to[i]] == lev[u] + ) {
tem = DFS(to[i], min(maxf, f[i]));
maxf -= tem;
f[i] -= tem;
f[i ^ ] += tem;
cnt += tem;
if (!maxf) {
break;
}
}
if (!cnt) {
lev[u] = -;
}
return cnt;
}
JQK Dinic() {
JQK ans = ;
while (BFS()) {
ans += DFS(S, INF);
}
return ans;
}
void del(int x) {
for (int i = Head[x]; i; i = nxt[i])
f[i] = f[i ^ ] = ;
}
void init(int S1, int T1, int S2, int T2) {
for (int i = ; i <= MAXP + ; i++) {
totflow[i] = Head[i] = ;
}
ed = ;
S = S1;
T = T1;
SS = S2;
TT = T2;
return;
}
inline JQK pushdownflow(int x) {
JQK sum = ;
for (int i = ; i <= x; i++) {
if (totflow[i] > ) {
addedge(SS, i, totflow[i]);
sum += totflow[i];
} else if (totflow[i] < ) {
addedge(i, TT, -totflow[i]);
}
}
return sum;//sum==dinic 循环流有解
}
JQK work() { //有源汇最大流
JQK ans = ;
JQK sum = pushdownflow(TT - );
addedge(T, S, INF);
swap(S, SS), swap(T, TT);
if (sum != Dinic()) { //Dinic(SS,TT) 求出一个有源汇有上下界可行流.此时的流不一定最大.
cout << - << endl;
return ;
}
ans = f[ed]; //可行流流量
f[ed] = f[ed ^ ] = ;
swap(S, SS), swap(T, TT); //Dinic(S,T)
ans += Dinic();//新增广出的s-t流量
cout << ans << endl;
for (int i = n + m + ; i <= n + m + k; i++) {
for (int j = Head[i]; j; j = nxt[j]) {
int v = to[j];
if (v >= n + m + k + && v <= n + m + * k && f[j] == ) {
nodecnt.p = plb[i].first;
nodecnt.s = slb[v].first;
nodecnt.l = slb[v].second;
anser.push_back(psl[nodecnt]);
}
}
}
sort(anser.begin(), anser.end());
for (int v : anser) {
cout << v << " ";
}
return ans;
}
}
int main() {
int u, v;
JQK c;
scanf("%d %d %d", &n, &m, &k);
dinic::MAXP = n + m + * k + ;
s1 = n + m + * k + ;
t1 = s1 + ;
s2 = t1 + ;
t2 = s2 + ;
dinic::init(s1, t1, s2, t2); //S T SS TT
for (int i = ; i <= n; i++)
dinic::add(s1, i, , INF);
for (int i = n + ; i <= n + m; i++)
dinic::add(i, t1, , INF);
for (int i = ; i <= k; i++) {
scanf("%d %d %d", &u, &v, &c);
nodecnt.p = u, nodecnt.s = v, nodecnt.l = c;
if (!psl[nodecnt])
psl[nodecnt] = i;
if (!pl[make_pair(u, c)]) {
pl[make_pair(u, c)] = n + m + i;
plb[n + m + i] = make_pair(u, c);
dinic::add(u, n + m + i, , );
}
if (!sl[make_pair(v, c)]) {
sl[make_pair(v, c)] = n + m + k + i;
slb[n + m + k + i] = make_pair(v, c);
dinic::add(n + m + k + i, n + v, , );
}
dinic::add(pl[make_pair(u, c)], sl[make_pair(v, c)], , );
}
dinic::work();
return ;
}

2019 Petrozavodsk Winter Camp, Yandex Cup C. Diverse Singing 上下界网络流的更多相关文章

  1. 2018 Petrozavodsk Winter Camp, Yandex Cup

    A. Ability Draft solved by RDC 60min start, 148 min AC, 1Y 题意:两只 Dota 队伍,每队 \(n\) 个英雄,英雄一开始无技能,他们需要按 ...

  2. 2019 wannafly winter camp day 3

    2019 wannafly winter camp day 3 J 操作S等价于将S串取反,然后依次遍历取反后的串,每次加入新字符a,当前的串是T,那么这次操作之后的串就是TaT.这是第一次转化. 涉 ...

  3. 2019 wannafly winter camp

    2019 wannafly winter camp Name Rank Solved A B C D E F G H I J K day1 9 5/11 O O O O O day2 5 3/11 O ...

  4. Petrozavodsk Winter Camp, Warsaw U, 2014, A The Carpet

    一个地图上有若干障碍,问允许出现一个障碍的最大子矩形为多大? 最大子矩形改编 #include<bits/stdc++.h> using namespace std; #define re ...

  5. Petrozavodsk Winter Camp, Day 8, 2014, Ship

    $dp(i,j)$表示i-j这段还没运走时的状态,包括 运输了多少次,还剩多少空间 每次枚举运输左边还是右边转移 #include <bits/stdc++.h> #define rep( ...

  6. Petrozavodsk Winter Camp, Day 8, 2014, Fine Brochures

    1的个数-块的个数+多减的个数+flag 多减的只会在一个循环末尾出现 #include <bits/stdc++.h> using namespace std; #define rep( ...

  7. Petrozavodsk Winter Camp, Day 8, 2014, Second Trip

    给你一棵树,每次询问一个(a,b),问有多少有路径与a-b没有交集 找lca #include <bits/stdc++.h> using namespace std; #define r ...

  8. Petrozavodsk Winter Camp, Day 8, 2014, Mosaic

    给你三个数字a,b,c,让你用1-m的数字凑出来 结论:有2个1和2个2肯定凑不出来,然后就搜索 #include <bits/stdc++.h> using namespace std; ...

  9. Petrozavodsk Winter Camp, Day 8, 2014, Rectangle Count

    给一个n*m的格点图,问其中有多少个矩形? $ \sum_{x=1}^{nm} \sum_{ab=x} [a + b \leq n](n - a - b + 1)\sum_{cd=x} [c + d ...

随机推荐

  1. Python数据结构与语法

    字典:Python字典是另一种可变容器模型,且可存储任意类型对象,如字符串.数字.元组.字典等其他容器模型:值可以取任何数据类型,但键必须是不可变的,如字符串,数字或元组:遍历字典时遍历的是键:访问v ...

  2. windows下exfat无法写入怎么修复?

    为了能够实现mac与windows文件共享,把移动硬盘格式化为exfat了,可是在osx中放入文件后,在windows上紧进行读取写入时出现错误,提示使用chkdsk进行修正,下面是修正步骤. 方法/ ...

  3. K/3 Cloud 单据关联查询

    销售出库单 下推 销售退货单,如何获知他们的关联关系?T_SAL_OUTSTOCKENTRY 是销售出库单分录T_SAL_RETURNSTOCKENTRY 是销售退货单分录T_SAL_RETURNST ...

  4. Java学习笔记-对象与垃圾回收

    Java存在垃圾回收机制,JVM会去回收垃圾,释放资源,而不是像C++一样有程序员去完成 垃圾回收机制的特点 垃圾回收机制只负责回收堆内存中的对象,不会回收任何物理资源(例如数据库连接.网络IO等资源 ...

  5. 监听input框变化,即时搜索 compositionstart, compositionend

    前话: 如果直接通过input监听, 它是键盘输入按键按下了就触发时间,这样可能导致一些问题,比如在谷歌浏览器想输入中文输入不了: 解决方案: 用 compositionstart 和 composi ...

  6. Mybatis插件之Mybatis-Plus的CRUD方法

    使用Mybatis-plus进行基本的CRUD(增查改删)操作. 实体类(User)代码: import com.baomidou.mybatisplus.annotation.IdType; imp ...

  7. [转帖]中国新超算彻底告别进口CPU 国产芯片已可与国外抗衡

    中国新超算彻底告别进口CPU 国产芯片已可与国外抗衡 蓝天·2017-10-17·本土IC 来源: 观察者网 https://www.laoyaoba.com/html/news/newsdetail ...

  8. 关于quartz定期,起服务时不新增配置文件中的定期问题

    关于quartz定期,起服务时不新增配置文件中的定期问题 问题原因:生产环境中起服务,未加载配置文件信息,且quartz连接超时 查找原因发现 由于别人新建了一个定期文件 并将 quartz工厂类的i ...

  9. WinForm打包

    首先要在想要打包的项目下创建一个新的项目, 创建好setup项目,之后点击属性,去修改打包软件的名字,ProductName....可以选填 到此已经创建好了setup工程了,那么下面开始将要打包的d ...

  10. AtCoder AISing Programming Contest 2019 Task D. Nearest Card Game

    题目分析在代码注释里. int main() { #if defined LOCAL && !defined DUIPAI ifstream in("main.in" ...