题面

背景


题目描述

从前有个

n

n

n 个方格排成一行,从左至右依此编号为

1

,

2

,

,

n

1,2,⋯,n

1,2,⋯,n。

有一天思考熊想给这

n

n

n 个方格染上黑白两色。

i

i

i 个方格上有

6

6

6 个属性:

a

i

,

b

i

,

w

i

,

l

i

,

r

i

,

p

i

a_i,b_i,w_i,l_i,r_i,p_i

ai​,bi​,wi​,li​,ri​,pi​。

如果方格

i

i

i 染成黑色就会获得

b

i

b_i

bi​ 的好看度。

如果方格

i

i

i 染成白色就会获得

w

i

w_i

wi​ 的好看度。

但是太多了黑色就不好看了。如果方格

i

i

i 是黑色,并且存在一个

j

j

j 使得

1

j

<

i

1≤j<i

1≤j<i 且

l

i

a

j

r

i

l_i≤a_j≤r_i

li​≤aj​≤ri​ 且方格

j

j

j 为白色,那么方格

i

i

i 就被称为奇怪的方格。

如果方格

i

i

i 是奇怪的方格,就会使总好看度减少

p

i

p_i

pi​。

也就是说对于一个染色方案,好看度为:

方格 i 为黑色

b

i

+

方格 i 为白色

w

i

方格 i 为奇怪的方格

p

i

∑_\text{方格 i 为黑色}b_i+∑_\text{方格 i 为白色}w_i−∑_\text{方格 i 为奇怪的方格}p_i

方格 i 为黑色∑​bi​+方格 i 为白色∑​wi​−方格 i 为奇怪的方格∑​pi​

现在给你

n

,

a

,

b

,

w

,

l

,

r

,

p

n,a,b,w,l,r,p

n,a,b,w,l,r,p,问所有染色方案中最大的好看度是多少。

数据范围

a

m

a

x

a_{max}

amax​ 为

a

,

l

,

r

a,l,r

a,l,r 中的最大值,

v

m

a

x

v_{max}

vmax​ 为

b

,

w

b,w

b,w 中的最大值,

p

m

a

x

p_{max}

pmax​ 为

p

p

p 中的最大值。

测试点编号

n

n

n

a

m

a

x

a_{max}

amax​

v

m

a

x

v_{max}

vmax​

p

m

a

x

p_{max}

pmax​

1

1

1

=

5

=5

=5

10

≤10

≤10

10

≤10

≤10

10

≤10

≤10

2

2

2

=

20

=20

=20

40

≤40

≤40

40

≤40

≤40

40

≤40

≤40

3

3

3

=

20

=20

=20

40

≤40

≤40

40

≤40

≤40

40

≤40

≤40

4

4

4

=

5000

=5000

=5000

10

≤10

≤10

200000

≤200000

≤200000

100000

≤100000

≤100000

5

5

5

=

5000

=5000

=5000

10

≤10

≤10

200000

≤200000

≤200000

300000

≤300000

≤300000

6

6

6

=

200

=200

=200

109

≤109

≤109

200000

≤200000

≤200000

200000

≤200000

≤200000

7

7

7

=

300

=300

=300

109

≤109

≤109

200000

≤200000

≤200000

220000

≤220000

≤220000

8

8

8

=

500

=500

=500

109

≤109

≤109

200000

≤200000

≤200000

400000

≤400000

≤400000

9

9

9

=

5000

=5000

=5000

5000

≤5000

≤5000

200000

≤200000

≤200000

150000

≤150000

≤150000

10

10

10

=

5000

=5000

=5000

109

≤109

≤109

200000

≤200000

≤200000

300000

≤300000

≤300000

2

s

,

48

M

B

\tt 2\,s~~,~~48\,MB

2s  ,  48MB

题解

最大的好看度,可以转化为:总好看度(含 b,w) - 最小损失好看度(含 b,w,p)

于是,可以用网络流最小割解决“最小损失好看度”的问题。

  • S

    ,

    T

    S,T

    S,T 为源点、汇点,

    S

    S

    S 部表示涂黑色,

    T

    T

    T 部表示涂白色。

  • 我们把每个点

    i

    i

    i 连两条边:

    S

    i

    S\rightarrow i

    S→i(边权为

    b

    i

    b_i

    bi​,若保留,则表示

    i

    i

    i 涂黑色),

    i

    T

    i\rightarrow T

    i→T(边权为

    w

    i

    w_i

    wi​ ,若保留,则表示

    i

    i

    i 涂白色)。

  • 对于每个

    i

    i

    i ,新建一个点

    i

    i'

    i′,与每个满足

    1

    j

    <

    i

    ,

    l

    i

    a

    j

    r

    i

    1\leq j<i,l_i\leq a_j\leq r_i

    1≤j<i,li​≤aj​≤ri​ 的

    j

    j

    j 连一条边

    i

    j

    i'\rightarrow j

    i′→j (边权为

    +

    +\infty

    +∞),然后连一条边

    i

    i

    i\rightarrow i'

    i→i′(边权为

    p

    i

    p_i

    pi​),这样,如果所有的

    j

    j

    j 有任意一个是白色(

    j

    T

    j\rightarrow T

    j→T 保留),为了使

    S

    ,

    T

    S,T

    S,T 不可达,就必须保留

    i

    T

    i\rightarrow T

    i→T 断掉

    S

    i

    S\rightarrow i

    S→i(选白色),或者保留

    S

    i

    S\rightarrow i

    S→i 断掉

    i

    i

    i\rightarrow i'

    i→i′ (选黑色,但是多付出

    p

    i

    p_i

    pi​ 的代价)。

虽然建图可以随便完成,但是边的数量太多,不仅空间存不下(48 Mb),而且最大流会超时。

于是我们可以优化建图

我们建立新点

i

i'

i′ 的目的,其实是要使得存在这么一个点,通过若干容量为无穷的边,可以到达所有符合

1

j

<

i

,

l

i

a

j

r

i

1\leq j<i,l_i\leq a_j\leq r_i

1≤j<i,li​≤aj​≤ri​ 的

j

j

j。但是,不同的

i

i'

i′ ,可能连向一大批相同的

j

j

j ,

j

j

j 和

a

j

a_j

aj​ 都是连续在一个范围内的

线段树优化建图

瞎扯了这么多,我也不绕弯了,用线段树优化建图。每次的

i

i'

i′ 只需要连向线段树上的某些(log的数量级)区间点,每个父亲向左右儿子连容量无穷的边,每个叶子向

a

j

a_j

aj​ 等于某一值的所有

j

j

j 连容量无穷的边。

有两个条件,

1

j

<

i

1\leq j<i

1≤j<i 和

l

i

a

j

r

i

l_i\leq a_j\leq r_i

li​≤aj​≤ri​ ,难道要树套树?

其实只用可持久化线段树就行了。

分块优化建图

n

5000

n\leq 5000

n≤5000 ,貌似有点小?根号和 log 差不太多啊。

可持久化分块,何如?每次只用新建 1~2 个点,边数也不用如此繁多。

实际上比线段树优秀得多:

CODE

线段树

#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<ctime>
#include<queue>
#include<bitset>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 5005
#define LL long long
#define ULL unsigned long long
#define UI unsigned int
#define DB double
#define ENDL putchar('\n')
#define lowbit(x) (-(x) & (x))
#define FI first
#define SE second
#define eps (1e-4)
#define SI(x) set<x>::iterator
#define MI map<int,int>::iterator
#define SQ 51
//#pragma GCC optimize(2)
//#pragma GCC optimize(3)
//#pragma GCC optimize("Ofast")
LL read() {
LL f=1,x=0;char s = getchar();
while(s < '0' || s > '9') {if(s=='-')f = -f;s = getchar();}
while(s >= '0' && s <= '9') {x=x*10+(s-'0');s = getchar();}
return f*x;
}
void putpos(LL x) {
if(!x) return ;
putpos(x/10); putchar('0'+(x%10));
}
void putnum(LL x) {
if(!x) putchar('0');
else if(x < 0) putchar('-'),putpos(-x);
else putpos(x);
}
void AIput(LL x,char c) {putnum(x);putchar(c);} int n,m,s,o,k;
int cnp,hd[MAXN*200],S,T;
int v[MAXN*400],nx[MAXN*400],rev[MAXN*400],cne;
LL w[MAXN*400];
int ins(int x,int y,LL k) {
if(!x || !y) return 0;
nx[++ cne] = hd[x];v[cne] = y;w[cne] = k;hd[x] = cne;
nx[++ cne] = hd[y];v[cne] = x;w[cne] = 0;hd[y] = cne;
rev[cne] = cne-1; rev[cne-1] = cne; return cne-1;
}
int hd2[MAXN*200],d[MAXN*200];
bool bfs(int S,int T) {
for(int i = 1;i <= cnp;i ++) {
hd2[i] = hd[i]; d[i] = -1;
}
queue<int> b;
d[S] = 0;
b.push(S);
while(!b.empty()) {
int t = b.front();b.pop();
if(t == T) return 1;
for(int i = hd[t];i;i = nx[i]) {
if(d[v[i]] < 0 && w[i] > 0) {
d[v[i]] = d[t] + 1;
b.push(v[i]);
}
}
}
return 0;
}
LL dfs(int x,LL flow) {
if(x == T) return flow;
LL res = 0;
for(int i = hd2[x];i;i = nx[i]) {
if(d[v[i]] == d[x] + 1 && w[i] > 0) {
LL mx = dfs(v[i],min(flow-res,w[i]));
res += mx; w[i] -= mx; w[rev[i]] += mx;
if(res == flow) break;
}
hd2[x] = nx[i];
}
return res;
}
LL dinic(int S,int T) {
LL ans = 0;
while(bfs(S,T)) {
ans += dfs(S,(LL)1e18);
}return ans;
} int A[MAXN],bk[MAXN],we[MAXN],ll[MAXN],rr[MAXN],P[MAXN];
int ar[MAXN];
int tre[MAXN<<2],M;
void maketree(int n) {
M=1;while(M<n+2)M<<=1;
}
void addtree(int x,int y) {
int s=M+x;ins(++ cnp,tre[s],(LL)1e18);
tre[s] = cnp;ins(tre[s],y,(LL)1e18);s >>= 1;
while(s) {
tre[s] = ++ cnp;
ins(tre[s],tre[s<<1],(LL)1e18);
ins(tre[s],tre[s<<1|1],(LL)1e18);
s >>= 1;
}return ;
}
int buildp(int l,int r) {
int p = ++ cnp;
int s = M+l-1,t = M+r+1;
while(s || t) {
if((s>>1) ^ (t>>1)) {
if(!(s&1)) ins(p,tre[s^1],(LL)1e18);
if(t & 1) ins(p,tre[t^1],(LL)1e18);
}else break;
s >>= 1;t >>= 1;
}return p;
}
int main() {
n = read();
cnp = n+2; S = n+1; T = n+2;
LL ans = 0;
for(int i = 1;i <= n;i ++) {
A[i] = read(); ar[i] = A[i];
bk[i] = read(); we[i] = read();
ll[i] = read(); rr[i] = read();
P[i] = read();
ins(S,i,bk[i]); ins(i,T,we[i]);
ans += we[i]+bk[i];
}
sort(ar + 1,ar + 1 + n);
for(int i = 1;i <= n;i ++) {
A[i] = lower_bound(ar + 1,ar + 1 + n,A[i]) - ar;
ll[i] = lower_bound(ar + 1,ar + 1 + n,ll[i]) - ar;
rr[i] = lower_bound(ar + 1,ar + 1 + n,rr[i]+1) - ar - 1;
}
maketree(n);
for(int i = 1;i <= n;i ++) {
int p = buildp(ll[i],rr[i]);
ins(i,p,P[i]); addtree(A[i],i);
}
ans -= dinic(S,T);
AIput(ans,'\n');
// for(int i = 1;i <= n;i ++) {
// printf("%d(%d): ",i,A[i]);
// if(w[fr[i]] > 0) {
// printf("black");
// }
// else printf("white");
// printf(" [%d,%d]\n",ll[i],rr[i]);
// }
return 0;
}

分块

#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<ctime>
#include<queue>
#include<bitset>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 5005
#define LL long long
#define ULL unsigned long long
#define UI unsigned int
#define DB double
#define ENDL putchar('\n')
#define lowbit(x) (-(x) & (x))
#define FI first
#define SE second
#define eps (1e-4)
#define SI(x) set<x>::iterator
#define MI map<int,int>::iterator
#define SQ 51
//#pragma GCC optimize(2)
//#pragma GCC optimize(3)
//#pragma GCC optimize("Ofast")
LL read() {
LL f=1,x=0;char s = getchar();
while(s < '0' || s > '9') {if(s=='-')f = -f;s = getchar();}
while(s >= '0' && s <= '9') {x=x*10+(s-'0');s = getchar();}
return f*x;
}
void putpos(LL x) {
if(!x) return ;
putpos(x/10); putchar('0'+(x%10));
}
void putnum(LL x) {
if(!x) putchar('0');
else if(x < 0) putchar('-'),putpos(-x);
else putpos(x);
}
void AIput(LL x,char c) {putnum(x);putchar(c);} int n,m,s,o,k;
int cnp,hd[MAXN*200],S,T;
int v[MAXN*400],nx[MAXN*400],rev[MAXN*400],cne;
LL w[MAXN*400];
int ins(int x,int y,LL k) {
if(!x || !y) return 0;
nx[++ cne] = hd[x];v[cne] = y;w[cne] = k;hd[x] = cne;
nx[++ cne] = hd[y];v[cne] = x;w[cne] = 0;hd[y] = cne;
rev[cne] = cne-1; rev[cne-1] = cne; return cne-1;
}
int hd2[MAXN*200],d[MAXN*200];
bool bfs(int S,int T) {
for(int i = 1;i <= cnp;i ++) {
hd2[i] = hd[i]; d[i] = -1;
}
queue<int> b;
d[S] = 0;
b.push(S);
while(!b.empty()) {
int t = b.front();b.pop();
if(t == T) return 1;
for(int i = hd[t];i;i = nx[i]) {
if(d[v[i]] < 0 && w[i] > 0) {
d[v[i]] = d[t] + 1;
b.push(v[i]);
}
}
}
return 0;
}
LL dfs(int x,LL flow) {
if(x == T) return flow;
LL res = 0;
for(int i = hd2[x];i;i = nx[i]) {
if(d[v[i]] == d[x] + 1 && w[i] > 0) {
LL mx = dfs(v[i],min(flow-res,w[i]));
res += mx; w[i] -= mx; w[rev[i]] += mx;
if(res == flow) break;
}
hd2[x] = nx[i];
}
return res;
}
LL dinic(int S,int T) {
LL ans = 0;
while(bfs(S,T)) {
ans += dfs(S,(LL)1e18);
}return ans;
} int bl[MAXN],cnb,tt[MAXN];
int A[MAXN],bk[MAXN],we[MAXN],ll[MAXN],rr[MAXN],P[MAXN];
int ar[MAXN];
int buildp(int l,int r) {
int s = l/SQ+1,t = r/SQ+1;
int p = ++ cnp;
if(s == t) {
for(int i = l;i <= r;i ++) ins(p,tt[i],(LL)1e18);
}
else {
for(int i = l;(i/SQ+1) == s;i ++) ins(p,tt[i],(LL)1e18);
for(int i = r;(i/SQ+1) == t;i --) ins(p,tt[i],(LL)1e18);
for(int i = s+1;i < t;i ++) {
ins(p,bl[i],(LL)1e18);
}
}return p;
}
int main() {
n = read();
cnp = n+2; S = n+1; T = n+2;
LL ans = 0;
for(int i = 1;i <= n;i ++) {
A[i] = read(); ar[i] = A[i];
bk[i] = read(); we[i] = read();
ll[i] = read(); rr[i] = read();
P[i] = read();
ins(S,i,bk[i]); ins(i,T,we[i]);
ans += we[i]+bk[i];
}
sort(ar + 1,ar + 1 + n);
for(int i = 1;i <= n;i ++) {
A[i] = lower_bound(ar + 1,ar + 1 + n,A[i]) - ar;
ll[i] = lower_bound(ar + 1,ar + 1 + n,ll[i]) - ar;
rr[i] = lower_bound(ar + 1,ar + 1 + n,rr[i]+1) - ar - 1;
}
for(int i = 1;i <= n;i ++) {
int p = buildp(ll[i],rr[i]);
ins(i,p,P[i]); int B = A[i]/SQ+1; cnb = max(cnb,B);
ins(++ cnp,bl[B],(LL)1e18);
bl[B] = cnp;
ins(bl[B],i,(LL)1e18);
ins(++ cnp,tt[A[i]],(LL)1e18);
tt[A[i]] = cnp;
ins(tt[A[i]],i,(LL)1e18);
}
ans -= dinic(S,T);
AIput(ans,'\n');
// for(int i = 1;i <= n;i ++) {
// printf("%d(%d): ",i,A[i]);
// if(w[fr[i]] > 0) {
// printf("black");
// }
// else printf("white");
// printf(" [%d,%d]\n",ll[i],rr[i]);
// }
return 0;
}

UOJ#XX A+B Problem (罔烙硫)的更多相关文章

  1. UOJ#77. A+B Problem [可持久化线段树优化建边 最小割]

    UOJ#77. A+B Problem 题意:自己看 接触过线段树优化建图后思路不难想,细节要处理好 乱建图无果后想到最小割 白色和黑色只能选一个,割掉一个就行了 之前选白色必须额外割掉一个p[i], ...

  2. ACM程序设计选修课——1041: XX's easy problem(神烦的多次字符串重定向处理)

    1041: XX's easy problem Time Limit: 1 Sec  Memory Limit: 128 MB Submit: 41  Solved: 7 [Submit][Statu ...

  3. BZOJ 3218 UOJ #77 A+B Problem (主席树、最小割)

    大名鼎鼎的A+B Problem, 主席树优化最小割-- 调题死活调不对,一怒之下改了一种写法交上去A了,但是改写法之后第4,5个点常数变大很多,于是喜提UOJ全站倒数第三 目前还不知道原来的写法为什 ...

  4. BZOJ3218 UOJ#77 A+B Problem(最小割+主席树)

    竟然在BZOJ上拿了Rank1太给力啦. p.s.:汗,一发这个就被一堆人在2月27号强势打脸-- 传送门(BZOJ) 传送门(UOJ) 说说这道题目吧: 首先是说说这个构图吧.因为有选择关系,我们很 ...

  5. UOJ#77. A+B Problem

    题目名称是吸引你点进来的. 从前有个 n 个方格排成一行,从左至右依此编号为 1,2,⋯,n. 有一天思考熊想给这 n 个方格染上黑白两色. 第 i 个方格上有 6 个属性:ai,bi,wi,li,r ...

  6. 【UOJ #246】【UER #7】套路

    http://uoj.ac/contest/35/problem/246 神奇!我这辈子是想不出这样的算法了. 对区间长度分类讨论:题解很好的~ 我已经弱到爆了,看完题解后还想了一晚上. 题解中&qu ...

  7. 【UOJ #244】【UER #7】短路

    http://uoj.ac/contest/35/problem/244 对其他人来说好简单的一道题,我当时却不会做TWT 注定滚粗啊 题解很好的~ #include<cstdio> #i ...

  8. UOJ 180【UR #12】实验室外的攻防战

    http://uoj.ac/contest/25/problem/180 从前往后对比串A,B 当$A_i,B_i$不相同时找到$B_i$在A中的位置j 若$min{A_1,A_2,A_3...... ...

  9. UOJ 58 (树上带修改的莫队)

    UOJ 58 糖果公园 Problem : 给一棵n个点的树,每个点上有一种颜色,对于一条路径上的点,若 i 颜色第 j 次出现对该路径权值的贡献为 w[i] * c[j], 每次询问一条路径的权值, ...

随机推荐

  1. 方法(method)

    方法是可以完成某个特定的功能,并且可以重复利用的代码片段...C中叫为函数 方法定义在类体中,不可定义在主方法下. 一个方法执行完就会被释放, 提高代码的复用性 相同的业务逻辑就可以不用重复,,,,因 ...

  2. rosbag遍历数据出错:(unicode error) 'utf-8' codec can't decode byte 0xcd in position 31: invalid continuation byte

    主题: 前言 针对ros系统记录的bag文件,可以使用python的rosbag包,按照不同起止时间和topic进行提取. 然而,有的topic可以使用rosbag读取,但是不能遍历,存在解码错误.原 ...

  3. 「非软文」零基础学习TypeScript(源码开源)

    今天,这篇文章篇幅很短,主要开放我最近学习整理TypeScript源码. 源码地址 https://github.com/maomincoding/typeScript_study 更多内容请见原文, ...

  4. hadoop集群搭建——单节点(伪分布式)

    1. 准备工作: 前提:需要电脑安装VM,且VM上安装一个Linux系统 注意:本人是在学习完尚学堂视频后,结合自己的理解,在这里做的总结.学习的视频是:大数据. 为了区分是在哪一台机器做的操作,eg ...

  5. uipath 如何利用函数split切割换行符?

    uipath 如何利用函数split切割换行符? 答案在这 https://rpazj.com/thread-178-1-1.html

  6. archlinux-小米pro15_2020款使用archlinux(MX350显卡驱动安装)

    1.官网下载archlinux ISO镜像 https://archlinux.org/download/   使用磁力链接下载 2.使用软碟通将镜像写入U盘,制作成U盘启动盘 3.进入BIOS 关掉 ...

  7. 我是如何将一个老系统的kafka消费者服务的性能提升近百倍的

    ☞☞☞ 我是如何将一个老系统的kafka消费者服务的性能提升近百倍的 ☜☜☜ ○○○○○○○○○○○○○○○ 大家好,又见面了~ kafka作为一种高吞吐量的分布式发布订阅消息系统,在业务系统中被广泛 ...

  8. Docker安装NextCloud使用MySQL

    安装 1.拉取并启动MySQL,最好把数据可目录挂载到宿主机,以便容器被误删后恢复: docker run --name=nextcloud_db \ -e MYSQL_ROOT_PASSWORD=X ...

  9. 数据类型 简单扩展(Java)

    public class HelloWorld { public static void main(String[] args) { //整数拓展 进制 二进制0b 十进制 八进制0 十六进制0x i ...

  10. docker容器管理操作

    Docker容器的四种状态: 运行 已暂停 重新启动 已退出 1.容器的创建 容器创建:就是将镜像加载到容器的过程. 创建容器时如果没有指定容器名称,系统会自动创建一个名称. 新创建的容器默认处于停止 ...