线段树维护最后一个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<<);
 const int inf=(<<);
 const double ESP=1e-;
 const int mod=(int)1e9+;
 const int N=(int)2e5+;

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

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

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

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

     Build(l,mid,rt<<);
     Build(mid+,r,rt<<|);
     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)>>;
     dn(rt);
     if(pos<=mid)
         update_dot(pos,V,l,mid,rt<<);
     else
         update_dot(pos,V,mid+,r,rt<<|);
     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)>>;
     dn(rt);
     if(L<=mid)
         update_qu(L,R,V,l,mid,rt<<);
     if(R>mid)
         update_qu(L,R,V,mid+,r,rt<<|);
     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)>>;

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

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

 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-);
 }

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

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

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