猴子拆房 题解
题目描述
输入
输出
样例输入
【样例输入1】 2 2 3 4 5 【样例输入2】 3 2 4 2 5 1 3 【样例输入3】 6 3 5 3 4 1 7 1 7 4 2 4 1样例输出
【样例输出1】 3 【样例输出2】 0 【样例输出3】 10数据范围限制
提示
这个是我的，是我的QWQ，我没有转载，只是把以前的博客搬运过来了。
然后就是因为markdown丢了，所以我截图：
不过代码还在OJ上存着，你们放心：
#include<cstdio>
#include<algorithm>
#define ll long long
#define INF 1e12
using namespace std;
ll n;
struct node{
friend bool operator <(const node &x,const node &y){ // 懒得改了
return x.v<y.v;
}
friend bool operator >(const node &x,const node &y){
return x.v>y.v;
}
ll h,v; // 小根堆 大根堆
} a[200010],b[200010],t1[200010],t2[200010];
ll num[200010];
ll cnt1,cnt2;
ll big,sum,mn=INF;
ll spend[200010];
ll s[200010];
// 下面是手写堆
// 下面的肯定对了！
void insert1(node k){
t1[++cnt1]=k;
ll x=cnt1;
while(x>1 && t1[x/2]>t1[x]){
swap(t1[x/2],t1[x]);
x/=2;
}
}
void pop1(){
swap(t1[1],t1[cnt1]);
cnt1--;
ll x=1;
while(2*x<=cnt1&&((2*x+1<=cnt1 && t1[x*2+1]<t1[x]) || t1[x*2]<t1[x])){
if(2*x+1>cnt1){
swap(t1[x*2],t1[x]);
x*=2;
}
else if(t1[x*2+1]<t1[x*2]){
swap(t1[x*2+1],t1[x]);
x=x*2+1;
}
else{
swap(t1[x*2],t1[x]);
x*=2;
}
}
}
void insert2(node k){
t2[++cnt2]=k;
ll x = cnt2;
while(x>1 && t2[x/2]<t2[x]){
swap(t2[x/2],t2[x]);
x/=2;
}
}
void pop2(){
swap(t2[1],t2[cnt2]);
cnt2--;
ll x=1;
while(2*x<=cnt2&&((2*x+1<=cnt2 && t2[x*2+1]>t2[x]) || t2[x*2]>t2[x])){
if(2*x+1>cnt2){
swap(t2[x*2],t2[x]);
x*=2;
}
else if(t2[x*2+1]>t2[x*2]){
swap(t2[x*2+1],t2[x]);
x=x*2+1;
}
else{
swap(t2[x*2],t2[x]);
x*=2;
}
}
}
void solve(ll l,ll r){
if(l==r) return;
ll mid=(l+r)>>1;
solve(l,mid);
solve(mid+1,r);
ll pos1=l,pos2=mid+1;
for(ll i=l;i<=r;i++){
if(pos2>r||(pos1<=mid && a[pos1].h<a[pos2].h)){
b[i]=a[pos1++];
}
else if(pos2<=r){
b[i]=a[pos2++];
}
}
for(ll i=l;i<=r;i++){
a[i]=b[i];
}
}
// 上面的肯定对了！
// 代码主要部分
void run(ll h){
ll k=(n-s[h])-(num[h]-1); // 最少拆多少个房子
while(k<cnt2 && cnt2){ // 多了
insert1(t2[1]);
sum-=t2[1].v;
pop2();
}
while(cnt2<k){ // 不够
if(!cnt1){
return;
}
insert2(t1[1]);
sum+=t1[1].v;
pop1();
}
big-=spend[h];
mn=min(mn,sum+big);
}
int main(){
freopen("house.in","r",stdin);
freopen("house.out","w",stdout);
scanf("%lld",&n);
for(ll i=1;i<=n;i++){
scanf("%lld %lld",&a[i].h,&a[i].v);
num[a[i].h]++;
spend[a[i].h]+=a[i].v;
}
solve(1,n);
for(ll i=100000;i>=1;i--){
s[i]=s[i+1]+num[i];
big+=spend[i];
}
ll i=1;
for(ll h=1;h<=100000;h++){
if(num[h]){
while(i<=n){
if(a[i].h>=h) break;
insert2(a[i]);
sum+=a[i].v;
insert1(t2[1]);
sum-=t2[1].v;
pop2();
i++;
}
run(h);
}
}
printf("%lld",mn);
}
另外再贴HLY大佬的题解一篇（单堆，不理解我的可以试着理解他的）：
虽然格式有点……但是思路很新奇。
代码帮你们从OJ上copy过来了，码风和我的不太一样：
#include<bits/stdc++.h>
using namespace std;
int n,b[100001],t[100001],tt[100001],cnt,sum,h=1,p[1000001],f[100001],k[100001],g,ans=23407234;
struct st{
int h;
int c;
}a[100001];
void up(long long x)
{
p[++cnt]=x;
long long wz=cnt;
while(wz>1&&p[wz]>p[wz/2])
{
swap(p[wz],p[wz/2]);
wz/=2;
}
}
long long down()
{
long long maxn=p[1],wz=1;
p[1]=p[cnt--];
while(wz*2<=cnt&&p[wz]<p[wz*2]
||wz*2+1<=cnt&&p[wz]<p[wz*2+1])
{
long long mx=p[wz*2],k=wz*2;
if(wz*2+1<=cnt&&p[wz*2]<p[wz*2+1])
{
mx=p[wz*2+1],k=wz*2+1;
}
swap(p[wz],p[k]);
wz=k;
}
return maxn;
}
void hb(int x,int y,int p,int q)
{
int i=x,j=p,k=x;
st b[200001];
while(i<=y&&j<=q)
{
if(a[i].h<a[j].h) b[k].c=a[i].c,b[k++].h=a[i++].h;
else b[k].c=a[j].c,b[k++].h=a[j++].h;
}
while(i<=y) b[k].c=a[i].c,b[k++].h=a[i++].h;
while(j<=q) b[k].c=a[j].c,b[k++].h=a[j++].h;
for(int l=x;l<=q;l++)
{
a[l].h=b[l].h,a[l].c=b[l].c;
}
}
void gb(int l,int r)
{
if(l==r) return;
int mid=(l+r)/2;
gb(l,mid);
gb(mid+1,r);
hb(l,mid,mid+1,r);
}
int main()
{
freopen("house.in","r",stdin);
freopen("house.out","w",stdout);
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%d%d",&a[i].h,&a[i].c);
}
gb(1,n);
for(int i=1;i<=n;i++)
{
if(t[a[i].h]==0)
{
b[++b[0]]=a[i].h;
}
t[a[i].h]++;
tt[a[i].h]+=a[i].c;
sum+=a[i].c;
}
for(int i=1;i<=b[0];i++)
{
int uu=b[i],g1=0;
sum-=tt[uu];
f[i]=sum;
while(a[h].h<uu) g+=a[h].c,up(a[h++].c);
if(t[uu]<=cnt)
{
for(int j=1;j<=t[uu]-1;j++)
{
k[++k[0]]=down();
g1+=k[k[0]];
}
f[i]+=g-g1;
for(int j=1;j<=k[0];j++)
{
up(k[j]);
}
k[0]=0;
}
ans=min(ans,f[i]);
}
printf("%d",ans);
}
另外手动膜拜HLY大佬，这里贴上他的博客
