BZOJ 4873 [Shoi2017]寿司餐厅 | 网络流 最大权闭合子图

链接

BZOJ 4873

题解

当年的省选题……还记得蒟蒻的我Day1 20分滚粗……

这道题是个最大权闭合子图的套路题。严重怀疑出题人就是先画好了图然后照着图编了个3000字的题面。和我喜欢的妹子当年给别人写的情书一样长……

最大权闭合子图

最大权闭合子图问题:一个有向图中,每个点带有一个权值(有正有负),有向边\(u \to v\)表示选\(u\)必须选\(v\),选出一些点使权值和最大,问权值和最大是多少。

最大权闭合子图的解法:网络流建图,对于每个点\(u\),设权值为\(w_u\),如果\(w_u\)为正,则连边\((S, u, w_u)\);若为负,连边\((u, T, |w_u|)\)。对于原图中每条有向边\(u \to v\),连边\((u, v, +\infty)\)。

答案即为所有正点权之和减去网络流的最小割。

这道题的依赖关系很明确:

  1. 选择了一个区间\([i, j](i < j)\),则必须选择它的两个子区间\([i, j - 1], [i + 1, j]\)。
  2. 选择了一个区间\([i, i]\),则必须选择里面的唯一寿司\(i\)。
  3. 选择了一个寿司\(i\),则必须选择它的代号\(a_i\)。

点的权值也很明确:

  1. “代号”\(i\)的权值为\(- m * i ^ 2\)。
  2. “代号”为\(i\)的寿司的权值为\(- i\)。
  3. 区间\([i, j]\)的权值就是输入数据中的“美味度”。

有了权值和依赖关系,按照上面的最大权闭合子图的建图套路无脑建图即可。

#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
typedef long long ll;
#define enter putchar('\n')
#define space putchar(' ')
template <class T>
void read(T &x){
    char c;
    bool op = 0;
    while(c = getchar(), c > '9' || c < '0')
    if(c == '-') op = 1;
    x = c - '0';
    while(c = getchar(), c >= '0' && c <= '9')
    x = x * 10 + c - '0';
    if(op) x = -x;
}
template <class T>
void write(T x){
    if(x < 0) putchar('-'), x = -x;
    if(x >= 10) write(x / 10);
    putchar('0' + x % 10);
}
const int N = 10005, M = 1000005, INF = 0x3f3f3f3f;
int n, m, src, des, sum;
int ecnt = 1, adj[N], cur[N], dis[N], nxt[M], go[M], cap[M];
int a[105], w[105][105], lst[1005], idx, type_num[1005], sushi_num[105], range_num[105][105];
void ADD(int u, int v, int _cap){
    go[++ecnt] = v;
    nxt[ecnt] = adj[u];
    adj[u] = ecnt;
    cap[ecnt] = _cap;
}
void add(int u, int v, int _cap){
    ADD(u, v, _cap);
    ADD(v, u, 0);
}
bool bfs(){
    static int que[N], qr;
    for(int i = 1; i <= des; i++)
    cur[i] = adj[i], dis[i] = -1;
    que[qr = 1] = src, dis[src] = 1;
    for(int ql = 1; ql <= qr; ql++){
    int u = que[ql];
    for(int e = adj[u], v; e; e = nxt[e])
        if(cap[e] && dis[v = go[e]] == -1){
        dis[v] = dis[u] + 1, que[++qr] = v;
        if(v == des) return 1;
        }
    }
    return 0;
}
int dfs(int u, int flow){
    if(u == des) return flow;
    int ret = 0, delta;
    for(int &e = cur[u], v; e; e = nxt[e])
    if(cap[e] && dis[v = go[e]] == dis[u] + 1){
        delta = dfs(v, min(cap[e], flow - ret));
        if(delta){
        cap[e] -= delta;
        cap[e ^ 1] += delta;
        ret += delta;
        if(ret == flow) return ret;
        }
    }
    dis[u] = -1;
    return ret;
}
int maxflow(){
    int ret = 0;
    while(bfs()) ret += dfs(src, INF);
    return ret;
}
int main(){
    read(n), read(m);
    for(int i = 1; i <= n; i++)
    read(a[i]), lst[++idx] = a[i];
    sort(lst + 1, lst + idx + 1);
    int tmp = idx = unique(lst + 1, lst + idx + 1) - lst - 1;
    for(int i = 1; i <= lst[idx]; i++)
    type_num[i] = lower_bound(lst + 1, lst + idx + 1, i) - lst;
    for(int i = 1; i <= n; i++)
    sushi_num[i] = ++idx;
    for(int i = 1; i <= n; i++)
    for(int j = i; j <= n; j++){
        range_num[i][j] = ++idx;
        read(w[i][j]);
    }
    src = ++idx, des = ++idx;
    for(int i = 1; i <= n; i++){
        add(sushi_num[i], des, a[i]);
    add(sushi_num[i], type_num[a[i]], INF);
    add(range_num[i][i], sushi_num[i], INF);
    }
    for(int i = 1; i <= n; i++)
    for(int j = i; j <= n; j++){
        if(j > i){
                add(range_num[i][j], range_num[i][j - 1], INF);
                add(range_num[i][j], range_num[i + 1][j], INF);
            }
            if(w[i][j] > 0) add(src, range_num[i][j], w[i][j]), sum += w[i][j];
            else add(range_num[i][j], des, -w[i][j]);
    }
    for(int i = 1; i <= tmp; i++)
        add(i, des, lst[i] * lst[i] * m);
    write(sum - maxflow()), enter;
    return 0;
}

原文地址:https://www.cnblogs.com/RabbitHu/p/BZOJ4873.html

时间: 2024-10-24 05:24:04

BZOJ 4873 [Shoi2017]寿司餐厅 | 网络流 最大权闭合子图的相关文章

BZOJ_4873_[Shoi2017]寿司餐厅_最大权闭合子图

题意:http://www.lydsy.com/JudgeOnline/problem.php?id=4873 分析:我们发现分数正负都有,并且之间有依赖关系,很容易想到最大权闭合子图. 建图: 1.S向正点连边,负点向T连边. 2.选了[i~j]显然要选[i+1~j]和[i~j-1],分别连边. 3.对于i==j的点,向对应的寿司连边. 4.总花费m*x*x+c*x拆成两部分.对于每个代号x,向T连容量为m*x*x的边,c*x这部分我们考虑算f[i][i]时把f[i][i]的值减掉x,当然也可

BZOJ:4873: [Shoi2017]寿司餐厅

4873: [Shoi2017]寿司餐厅 首先很开心在膜你赛的时候做了出来. 看到数据范围,看到不能dp,看到贡献去重后计算,咦,流? 那就容易了,转最大权闭合子图,每个区间建一个点,取了就一定要取他的子区间(依赖关系),代价上也很容易用依赖关系搞. 提交完A了就没理,后来同学说我#1了??? #include<cstdio> #include<algorithm> #define MN 40001 using namespace std; int read_p,read_ca,r

【暖*墟】#网络流# 最大权闭合子图

[相关概念详解] 闭合图:有向图的一个点集,且这个点集的所有出边仍然指向该点集. 最大权闭合图:(每一个点有一个权值)在所有的合法闭合图中,点权之和最大的图. 处理问题:权值有正有负,重复选只算一次,选择有相互关联性 的问题. 首先有一个有向连通图(闭合图),每个点带有一个权值,例如: 造出一个超级源点S和一个超级汇点T,把S连边到所有带有正权的点上,每条边的容量是这个点的权: 把所有带负权的点连边到T,每条边的容量是这个点的权的相反数(正值).原来的边的容量设成无限大. 所有的点按权值的正负连

BZOJ 4873 寿司餐厅(最大权闭合图 网络流)

寿司餐厅 时间限制: 1 Sec  内存限制: 512 MB提交: 6  解决: 3[提交][状态][讨论版] 题目描述 Kiana 最近喜欢到一家非常美味的寿司餐厅用餐.每天晚上,这家餐厅都会按顺序提供n种寿司,第i种寿司有一个代号 ai和美味度 di,i,不同种类的寿司有可能使用相同的代号.每种寿司的份数都是无限的,Kiana 也可以无限次取寿司来吃,但每种寿司每次只能取一份,且每次取走的寿司必须是按餐厅提供寿司的顺序连续的一段,即 Kiana 可以一次取走第 1,2种寿司各一份,也可以一次

bzoj 1497 [NOI2006]最大获利【最大权闭合子图+最小割】

不要被5s时限和50000点数吓倒!大胆网络流!我一个5w级别的dinic只跑了1s+! 看起来没有最大权闭合子图的特征--限制,实际上还是有的. 我们需要把中转站看成负权点,把p看成点权,把客户看成正权点,把c看成点权,然后把中转站点a.b作为客户点的依赖点 s点向所有正权点连边,流量为点权:所有负权点向t连边,流量为负点权(即正数!) 对于所有有依赖关系的点,由客户点向中转站点连边,流量为inf,也就是最大权闭合子图中的向其依赖点连边 连边的意义详见:http://www.cnblogs.c

BZOJ 1565 植物大战僵尸(拓扑排序+最大权闭合子图)

图中的保护关系就类似于最大权闭合子图.即你想杀x,你就一定要杀掉保护x的点,那么把x向保护它的点连边.那么题目就转化成了最大权闭合子图的问题. 但是这个图有点特殊啊... 考虑有环的情况,显然这个环以及指向这个环的点都不能选. 所以还要把这个图的反图进行一遍拓扑排序,这样忽略掉了这些点了... # include <cstdio> # include <cstring> # include <cstdlib> # include <iostream> # i

FZU2295 Human life:网络流-最大权闭合子图-二进制优化-第九届福建省大学生程序设计竞赛

目录 Catalog Solution: (有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 Catalog Problem:Portal传送门 ?原题目描述在最下面. ?题意就是很裸的最大权闭合子图. ?推荐阅读:胡伯涛<最小割模型在信息学竞赛中的应用> ?完完全全的模板题:新疆大学五月月赛-D-勤奋的杨老师 ?本题题意:m(50)个任务,n个技能.完成每个任务由相应的收益,完成每个任务前必须学一些技能.有些技能由先修技能. ?有些任务不能同时完成. Solution: ?训练

我和最大权闭合子图

第一次接触最大权闭合子图大概是2017年3月27号星期一,那段时间有5个同学(ZJC/LKQ/LWD/WJJ/...)去了湖南师大附中听PTY.BK他们讲课,因为我没有去Hfu一直很怪罪.后来安排我和LXY在70去电子科大之前去考一周的试,上午LXY的Mom接送,下午我们坐metro到世纪城再让LXY的Mom送回学校上晚自习.那一天我去70,T1是分治Floyd,T2是轮廓线DP,T3是一道最大权闭合子图(后来知道是BZOJ 3774 最优选择).那时我只会Dinic的模板,我的初高中学长Ana

【最大权闭合子图】bzoj4873: [Shoi2017]寿司餐厅

为什么跑得这么慢 Description Kiana最近喜欢到一家非常美味的寿司餐厅用餐.每天晚上,这家餐厅都会按顺序提供n种寿司,第i种寿司有一个 代号ai和美味度di,i,不同种类的寿司有可能使用相同的代号.每种寿司的份数都是无限的,Kiana也可以无限次 取寿司来吃,但每种寿司每次只能取一份,且每次取走的寿司必须是按餐厅提供寿司的顺序连续的一段,即Kiana 可以一次取走第1,2种寿司各一份,也可以一次取走第2,3种寿司各一份,但不可以一次取走第1,3种寿司.由于餐 厅提供的寿司种类繁多,