bzoj 1178: [Apio2009]CONVENTION会议中心(少见做法掉落!)【贪心+二分】
数组若干+手动二分一个的算法,bzoj rank8
=废话分割线=====
我我我我我!一定要说一下我的心路历程!我只用了几个数组和一个手动二分!在洛谷和bzoj都过了所以应该是对的!(跑的挺快甚至和学长合了影
事情是这样的,我首先在洛谷瞎贪心贪了55,然后调不出来去看正解,发现好麻烦啊,同时觉得我的贪心还挺对的,于是搞了一份标程拍,拍着拍着就拍出错了,发现判断字典序不太对,我忘记排过序了所以不能直接取尾部,于是顺着pr往前跳取min,于是A了!bzoj rank15欸!看到了两个学长欸!卡卡常卡卡常,于是升到rank8并且为了全场最短改了代码风格
=废话分割线=====
先打上id顺序按r排序,f存当前最多能接待的公司,s存前缀f最大值(值相同取字典序最小),p存s最大值位置,pr存这个点是从哪个点转移过来的。然后对第i个l,二分找出最右的r小于当前l,设为w,用s[w]更新当前数组f[i],连上前缀pr[i]=p[w],然后分三种情况修改s和p数组。
1、s[i-1]<f[i],直接s[i]=s[i-1],p[i]=p[i-1]即可;
2、s[i-1]>f[i],即当前前缀和最大值为当前i点,所以s[i]=f[i],p[i]=i;
3、对于相等的情况对于这两个点i和p[i-1],沿着pr往前跳,直到pr相等(因为pr相等之后两个点经过的点就相同了),对于经过的点的id取min,于是这两个min就是升序排序后最早不同的方案位置,所以比较这两个来决定转到情况1或2.
然后就没了。是不是超简单?
但是还是去学一下思路神奇的正解吧
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=200005;
int n,s[N],p[N],q[N],top,pr[N];
struct qwe {
int l,r,id;
} a[N];
bool cmp(const qwe &a,const qwe &b) {
return a.r<b.r;
}
int read() {
int r=0,f=1;
char p=getchar();
while(p>'9'||p<'0') {
if(p=='-') f=-1;
p=getchar();
}
while(p>='0'&&p<='9') {
r=r*10+p-48;
p=getchar();
}
return r*f;
}
int main() {
n=read();
for(int i=1; i<=n; i++) a[i].l=read(),a[i].r=read(),a[i].id=i;
sort(a+1,a+1+n,cmp);
for(int i=1; i<=n; i++) {
int w=0,l=0,r=i-1,k=a[i].l;
while(l<=r) {
int mid=(l+r)>>1;
if(a[mid].r<k) w=mid,l=mid+1;
else r=mid-1;
}
int now=s[w]+1;
pr[i]=p[w];
if(s[i-1]>now) s[i]=s[i-1],p[i]=p[i-1];
else if(s[i-1]==now) {
int p1=p[i-1],p2=i,mn1=1e9,mn2=1e9;
while(pr[p1]!=pr[p2]) {
if(a[p1].id<mn1) mn1=a[p1].id;
if(a[p2].id<mn2) mn2=a[p2].id;
p1=pr[p1],p2=pr[p2];
}
if(a[p1].id<mn1) mn1=a[p1].id;
if(a[p2].id<mn2) mn2=a[p2].id;
if(mn1<mn2) s[i]=s[i-1],p[i]=p[i-1];
else s[i]=now,p[i]=i;
}
else s[i]=now,p[i]=i;
}
printf("%d\n",s[n]);
int now=p[n];
while(now)
q[++top]=a[now].id,now=pr[now];
sort(q+1,q+1+top);
for(int i=1; i<=top; i++) printf("%d ",q[i]);
return 0;
}