#77. A+B Problem
题目名称是吸引你点进来的。
从前有个 nn 个方格排成一行,从左至右依此编号为 1,2,⋯,n1,2,⋯,n。
有一天思考熊想给这 nn 个方格染上黑白两色。
第 ii 个方格上有 66 个属性:ai,bi,wi,li,ri,piai,bi,wi,li,ri,pi。
如果方格 ii 染成黑色就会获得 bibi 的好看度。
如果方格 ii 染成白色就会获得 wiwi 的好看度。
但是太多了黑色就不好看了。如果方格 ii 是黑色,并且存在一个 jj 使得 1≤j<i1≤j<i 且 li≤aj≤rili≤aj≤ri 且方格 jj 为白色,那么方格 ii 就被称为奇怪的方格。
如果方格 ii 是奇怪的方格,就会使总好看度减少 pipi。
也就是说对于一个染色方案,好看度为:
∑方格i为黑色bi+∑方格i为白色wi−∑方格i为奇怪的方格pi∑方格i为黑色bi+∑方格i为白色wi−∑方格i为奇怪的方格pi
现在给你 n,a,b,w,l,r,pn,a,b,w,l,r,p,问所有染色方案中最大的好看度是多少。
输入格式
第一行一个正整数 nn。
接下来 nn 行中第 ii 行有用空格隔开的 66 个非负整数依次表示 ai,bi,wi,li,ri,piai,bi,wi,li,ri,pi。
保证 li≤rili≤ri。
输出格式
一个非负整数表示所有染色方案中最大的好看度。
样例一
input
10
0 1 7 3 9 2
7 4 0 9 10 5
1 0 4 2 10 2
7 9 1 5 7 2
6 3 5 3 6 2
6 6 4 1 8 1
6 1 6 0 6 5
2 2 5 0 9 3
5 1 3 0 2 5
5 6 7 1 1 2
output
55
explanation
最优染色方案为:白 黑 白 黑 白 黑 白 白 白 白
可以发现只有方格 66 为奇怪的方格。
所以好看度为:
==w1+b2+w3+b4+w5+b6+w7+w8+w9+w10−p67+4+4+9+5+6+6+5+3+7−155w1+b2+w3+b4+w5+b6+w7+w8+w9+w10−p6=7+4+4+9+5+6+6+5+3+7−1=55
限制与约定
设 amaxamax 为 a,l,ra,l,r 中的最大值,vmaxvmax 为 b,wb,w 中的最大值, pmaxpmax 为 pp 中的最大值。
测试点编号 | nn | amaxamax | vmaxvmax | pmaxpmax |
---|---|---|---|---|
1 | =5=5 | ≤10≤10 | ≤10≤10 | ≤10≤10 |
2 | =20=20 | ≤40≤40 | ≤40≤40 | ≤40≤40 |
3 | =20=20 | ≤40≤40 | ≤40≤40 | ≤40≤40 |
4 | =5000=5000 | ≤10≤10 | ≤200000≤200000 | ≤100000≤100000 |
5 | =5000=5000 | ≤10≤10 | ≤200000≤200000 | ≤300000≤300000 |
6 | =200=200 | ≤109≤109 | ≤200000≤200000 | ≤200000≤200000 |
7 | =300=300 | ≤109≤109 | ≤200000≤200000 | ≤220000≤220000 |
8 | =500=500 | ≤109≤109 | ≤200000≤200000 | ≤400000≤400000 |
9 | =5000=5000 | ≤5000≤5000 | ≤200000≤200000 | ≤150000≤150000 |
10 | =5000=5000 | ≤109≤109 | ≤200000≤200000 | ≤300000≤300000 |
时间限制:2s2s
空间限制:48MB48MB
来源
VFleaKing
傻到了一种境界。。。。。
线段树insert的操作的本质是修改,所以当你修改叶子节点时如果这个叶子节点本来就有值,你就是把它覆盖了
所以你这个叶子要向原来的叶子连边。
#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#define maxn 300005
#define maxm 1000005
#define LL long long
#define inf 999999999
using namespace std;
int n;
int info[maxn],Prev[maxm],to[maxm],cap[maxm],cnt_e=1;
inline void Node(int u,int v,int c){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v,cap[cnt_e]=c; }
inline void Line(int u,int v,int c){ Node(u,v,c),Node(v,u,0); }
int s,t,stm,vd[maxn],d[maxn];
int sa[maxm];
int lc[maxn],rc[maxn],tot,rt[maxn];
inline int NewNode()
{
return ++tot;
}
void Insert(int &now,int l,int r,int pos,int id)
{
int tmp=now;
now=NewNode();
lc[now]=lc[tmp],rc[now]=rc[tmp];
if(l==r)
{
Line(now,id,inf);
if(tmp) Line(now,tmp,inf);
return;
}
int mid=(l+r)>>1;
if(pos<=mid) Insert(lc[now],l,mid,pos,id);
else Insert(rc[now],mid+1,r,pos,id);
if(lc[now]) Line(now,lc[now],inf);
if(rc[now]) Line(now,rc[now],inf);
}
void connect(int now,int l,int r,int ql,int qr,int id)
{
if( qr<l || ql>r) return;
if(ql<=l && r<=qr)
{
Line(id,now,inf);
return;
}
int mid=(l+r)>>1;
if(lc[now])connect(lc[now],l,mid,ql,qr,id);
if(rc[now])connect(rc[now],mid+1,r,ql,qr,id);
}
int aug(int now,int Max)
{
if(now==t)
return Max;
int inc,st=Max,Min=tot-1;
for(int i=info[now];i;i=Prev[i])
if(cap[i])
{
if(d[to[i]]+1==d[now])
{
inc=aug(to[i],min(st,cap[i]));
st-=inc,cap[i]-=inc,cap[i^1]+=inc;
if(!st || d[s]>=tot) return Max-st;
}
Min=min(Min,d[to[i]]);
}
if(Max==st)
{
(!--vd[d[now]]) && (d[s]=tot);
vd[d[now]=Min+1]++;
}
return Max-st;
}
int sap()
{
memset(vd,0,sizeof vd);
memset(d,0,sizeof d);
for(vd[0]=tot,stm=0;d[s]<tot;)
stm+=aug(s,inf);
return stm;
}
int a[5005],b[5005],w[5005],l[5005],r[5005],p[5005];
int main()
{
#ifndef ONLINE_JUDGE
freopen("1.in","r",stdin);
#endif
int ans=0;
scanf("%d",&n);
tot=n*2;
s=NewNode(),t=NewNode();
for(int i=1;i<=n;i++)
{
scanf("%d%d%d%d%d%d",&a[i],&b[i],&w[i],&l[i],&r[i],&p[i]);
sa[++sa[0]]=a[i];
sa[++sa[0]]=l[i];
sa[++sa[0]]=r[i];
}
sort(sa+1,sa+sa[0]+1);
sa[0]=unique(sa+1,sa+1+sa[0])-sa-1;
for(int i=1;i<=n;i++)
{
a[i]=lower_bound(sa+1,sa+1+sa[0],a[i])-sa;
l[i]=lower_bound(sa+1,sa+1+sa[0],l[i])-sa;
r[i]=lower_bound(sa+1,sa+1+sa[0],r[i])-sa;
Line(s,i,b[i]);
Line(i,t,w[i]);
Line(i,i+n,p[i]);
ans+=b[i]+w[i];
}
for(int i=1;i<=n;i++)
{
if(i>1)connect(rt[i-1],1,sa[0],l[i],r[i],i+n);
rt[i]=rt[i-1];
Insert(rt[i],1,sa[0],a[i],i);
//ans+=p;
}
printf("%d\n",ans-sap());
}
话说这个主席树只维护线段间的二分关系,感觉和LCT维护LCA一样大材小用。。。。。。
20190314upd:
又错在了同一个地方,代码变好(难)看了。
#include<bits/stdc++.h>
#define maxn 200005
#define maxm maxn*6
#define inf 0x3f3f3f3f
using namespace std;
int n;
int buf[maxn],info[maxn],Prev[maxm],to[maxm],cap[maxm],cnt_e=1;
void Node(int u,int v,int c){ Prev[++cnt_e]=info[u],info[u]=cnt_e,cap[cnt_e]=c,to[cnt_e]=v; }
void Line(int u,int v,int c){ Node(u,v,c),Node(v,u,0); }
int S,T,tot,h[maxn],gap[maxn];
int rt[5005],lc[maxn],rc[maxn];
void TLine(int now,int l,int r,int ql,int qr,int p,int pos)
{
if(!now || r<ql || l>qr) return;
if(ql<=l && r<=qr)
{
Line(pos,now,inf);
return;
}
int mid = (l+r) >> 1;
TLine(lc[now],l,mid,ql,qr,p,pos);
TLine(rc[now],mid+1,r,ql,qr,p,pos);
}
void Insert(int &now,int l,int r,int pos,int id)
{
int tmp = now;
++tot , lc[tot] = lc[now] , rc[tot] = rc[now];
now = tot;
if(l == r)
{
Line(now,id,inf);
if(tmp)Line(now,tmp,inf);
return;
}
int mid = (l+r) >> 1;
if(pos<=mid) Insert(lc[now],l,mid,pos,id);
else Insert(rc[now],mid+1,r,pos,id);
if(lc[now]) Line(now,lc[now],inf);
if(rc[now]) Line(now,rc[now],inf);
}
int aug(int now,int Max)
{
if(now == T) return Max;
int inc , st = Max;
for(int &i=info[now];i;i=Prev[i])
if(cap[i] && h[to[i]]+1==h[now])
{
inc = aug(to[i],min(cap[i],st));
st -= inc, cap[i] -= inc , cap[i^1] += inc;
if(!st || h[S]>tot) return Max-st;
}
if(!--gap[h[now]]) h[S] = tot+1;
++gap[++h[now]],info[now]=buf[now];
return Max-st;
}
int main()
{
scanf("%d",&n) , tot = n;
S = ++tot , T = ++tot;
int ans = 0 , Mx = 1e9;
for(int i=1;i<=n;i++)
{
int a,b,w,l,r,p;
scanf("%d%d%d%d%d%d",&a,&b,&w,&l,&r,&p);
Line(S,i,b),Line(i,T,w),ans+=b+w;
TLine(rt[i-1],0,Mx,l,r,p,++tot),Line(i,tot,p);
Insert(rt[i]=rt[i-1],0,Mx,a,i);
}
memcpy(buf,info,sizeof info);
for(gap[0]=tot;h[S]<=tot;)
ans -= aug(S,inf);
printf("%d\n",ans);
}
UPD 20191207
#include<bits/stdc++.h>
#define maxn 5005
#define maxp maxn * 17
#define maxm maxp * 20
#define inf 0x3f3f3f3f
#define Clear(a,b) memset(a,b,sizeof a)
#define Copy(a,b) memcpy(a,b,sizeof a)
using namespace std;
int n,a[maxn],b[maxn],w[maxn],l[maxn],r[maxn],p[maxn],sb[maxn*3];
int info[maxp],Prev[maxm],to[maxm],cap[maxm],cnt_e=1,tot=0;
void Node(int u,int v,int c){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v,cap[cnt_e]=c; }
void Line(int u,int v,int c){ Node(u,v,c),Node(v,u,0); }
int rt,lc[maxp],rc[maxp];
int S,T,h[maxp],g[maxp],buf[maxp];
void ade(int u,int l,int r,int ql,int qr,int p){
if(l>qr||ql>r||!u) return;
if(ql<=l && r<=qr) return (void)(Line(p,u,inf));
int mid=(l+r)>>1;
ade(lc[u],l,mid,ql,qr,p),ade(rc[u],mid+1,r,ql,qr,p);
}
void ins(int &u,int l,int r,int p,int d){
lc[++tot]=lc[u],rc[tot]=rc[u];
if(l==r) return (void)((u)&&((Line(tot,u,inf)),0),Line(tot,d,inf),u=tot);
u=tot;
int mid=(l+r)>>1;
if(p<=mid) ins(lc[u],l,mid,p,d);
else ins(rc[u],mid+1,r,p,d);
if(lc[u]) Line(u,lc[u],inf);
if(rc[u]) Line(u,rc[u],inf);
}
int aug(int u,int mx){
if(u==T) return mx;
int st=mx,inc;
for(int &i=info[u],v;i;i=Prev[i])
if(cap[i] && h[v=to[i]]+1==h[u] && (inc=aug(v,min(cap[i],st)))){
cap[i]-=inc,cap[i^1]+=inc;
if(!(st-=inc) || h[0]) return mx-st;
}
if(!--g[h[u]]) h[0]=1;
++g[++h[u]];info[u]=buf[u];
return mx-st;
}
int main(){
scanf("%d",&n);tot = n;S=++tot,T=++tot;
for(int i=1;i<=n;i++) scanf("%d%d%d%d%d%d",&a[i],&b[i],&w[i],&l[i],&r[i],&p[i]),sb[++sb[0]]=l[i],sb[++sb[0]]=r[i],sb[++sb[0]]=a[i];
sort(sb+1,sb+1+sb[0]),sb[0]=unique(sb+1,sb+1+sb[0])-sb-1;
int stm=0;
for(int i=1;i<=n;i++){
Line(S,i,b[i]),Line(i,T,w[i]);stm+=b[i]+w[i];
a[i] = lower_bound(sb+1,sb+1+sb[0],a[i])-sb;
l[i] = lower_bound(sb+1,sb+1+sb[0],l[i])-sb;
r[i] = lower_bound(sb+1,sb+1+sb[0],r[i])-sb;
ade(rt,1,sb[0],l[i],r[i],++tot);
Line(i,tot,p[i]);
ins(rt,1,sb[0],a[i],i);
}
for(Clear(h,0),Clear(g,0),Copy(buf,info);!h[0];) stm-=aug(S,inf);
printf("%d\n",stm);
}