线段树维护最后一个0的位置（Restore Permutation）Manthan, Codefest 19 (open for everyone, rated, Div. 1 + Div. 2)

题意：https://codeforc.es/contest/1208/problem/D

给你长度为n的序列，s[i]的值为p[1]到p[i-1]中比p[i]小的数的和，让你求出p序列。

思路：

首先我们要想到，最后一个0的位置一定就是当前剩余还没用确定的数里的最小值，所以从1，2，3 。。。一直下去就行了。

选过了就在线段树上改为INF，同时pushup维护的最右边为0的位置。

一开始单点修改忘记down下去了，改了好久，学到了学到了。

#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
#include <cstdio>//sprintf islower isupper
#include <cstdlib>//malloc  exit strcat itoa system("cls")
#include <iostream>//pair
#include <fstream>//freopen("C:\Users\13606\Desktop\草稿.txt","r",stdin);
#include <bitset>
//#include <map>
//#include<unordered_map>
#include <vector>
#include <stack>
#include <set>
#include <string.h>//strstr substr
#include <string>
#include <time.h>//srand(((unsigned)time(NULL))); Seed n=rand()%10 - 0~9;
#include <cmath>
#include <deque>
#include <queue>//priority_queue<int, vector<int>, greater<int> > q;//less
#include <vector>//emplace_back
//#include <math.h>
//#include <windows.h>//reverse(a,a+len);// ~ ! ~ ! floor
#include <algorithm>//sort + unique : sz=unique(b+1,b+n+1)-(b+1);+nth_element(first, nth, last, compare)
using namespace std;//next_permutation(a+1,a+1+n);//prev_permutation
//******************
int abss(int a);
int lowbit(int n);
int Del_bit_1(int n);
int maxx(int a,int b);
int minn(int a,int b);
double fabss(double a);
void swapp(int &a,int &b);
clock_t __STRAT,__END;
double __TOTALTIME;
void _MS(){__STRAT=clock();}
void _ME(){__END=clock();__TOTALTIME=(double)(__END-__STRAT)/CLOCKS_PER_SEC;cout<<"Time: "<<__TOTALTIME<<" s"<<endl;}
//***********************
#define rint register int
#define fo(a,b,c) for(rint a=b;a<=c;++a)
#define fr(a,b,c) for(rint a=b;a>=c;--a)
#define mem(a,b) memset(a,b,sizeof(a))
#define pr printf
#define sc scanf
#define ls rt<<1
#define rs rt<<1|1
typedef long long ll;
const double E=2.718281828;
const double PI=acos(-1.0);
const ll INF=(1LL<<60);
const int inf=(1<<30);
const double ESP=1e-9;
const int mod=(int)1e9+7;
const int N=(int)2e5+10;

ll a[N];
ll add[N<<2],min_[N<<2];
int POS[N<<2];

void up(int rt,int l,int r)
{
    min_[rt]=min(min_[rt<<1],min_[rt<<1|1]);
    if(min_[ls]<min_[rs]) POS[rt]=POS[ls];
    else POS[rt]=POS[rs];
}

void dn(int rt)
{
    if(add[rt]!=0)
    {
        add[ls]+=add[rt];
        add[rs]+=add[rt];
        min_[rt<<1]+=add[rt];
        min_[rt<<1|1]+=add[rt];
        add[rt]=0;
    }
}

void Build(int l,int r,int rt)
{
    if(l==r)
    {
        min_[rt]=a[l];
        POS[rt]=l;
        add[rt]=0;
        return;
    }
    int mid=(l+r)>>1;

    Build(l,mid,rt<<1);
    Build(mid+1,r,rt<<1|1);
    up(rt,l,r);
}

void update_dot(int pos,ll V,int l,int r,int rt)
{
    if(l==r)
    {
        min_[rt]=V;
        return;
    }

    int mid=(l+r)>>1;
    dn(rt);
    if(pos<=mid)
        update_dot(pos,V,l,mid,rt<<1);
    else
        update_dot(pos,V,mid+1,r,rt<<1|1);
    up(rt,l,r);
}
void update_qu(int L,int R,int V,int l,int r,int rt)
{
    if(L>R)return;
    if(L<=l&&r<=R)
    {
        min_[rt]+=V;
        add[rt]+=V;
        return;
    }

    int mid=(l+r)>>1;
    dn(rt);
    if(L<=mid)
        update_qu(L,R,V,l,mid,rt<<1);
    if(R>mid)
        update_qu(L,R,V,mid+1,r,rt<<1|1);
    up(rt,l,r);
}
void check(int pos,int l,int r,int rt)
{
    if(l==r)
    {
        cout<<min_[rt]<<' ';
        return ;
    }
    int mid=(l+r)>>1;

    dn(rt);
    if(pos<=mid)
        check(pos,l,mid,rt<<1);
    else
        check(pos,mid+1,r,rt<<1|1);
}
int ans[N];
int main()
{
//    fill(min_,min_+N-3,INF);
    int n;
    sc("%d",&n);
    for(int i=1;i<=n;++i)
        sc("%lld",&a[i]);
    Build(1,n,1);
    for(int i=1;i<=n;++i)
    {
        int p=POS[1];
    //    pr("%d: %d
",i,p);
        ans[p]=i;
        update_dot(p,INF,1,n,1);
    //    for(int j=1;j<=n;++j)check(j,1,n,1);cout<<endl;
        update_qu(p+1,n,-i,1,n,1);
    //    for(int j=1;j<=n;++j)check(j,1,n,1);cout<<endl;
    }
    for(int i=1;i<=n;++i)
        pr("%d ",ans[i]);
    return 0;
}

/**************************************************************************************/

int maxx(int a,int b)
{
    return a>b?a:b;
}

void swapp(int &a,int &b)
{
    a^=b^=a^=b;
}

int lowbit(int n)
{
    return n&(-n);
}

int Del_bit_1(int n)
{
    return n&(n-1);
}

int abss(int a)
{
    return a>0?a:-a;
}

double fabss(double a)
{
    return a>0?a:-a;
}

int minn(int a,int b)
{
    return a<b?a:b;
}