【NOIP2016提高A组集训第12场11.10】灵知的太阳信仰

Description

在炽热的核熔炉中，居住着一位少女，名为灵乌路空。
据说，从来没有人敢踏入过那个熔炉，因为人们畏缩于空所持有的力量——核能。
核焰，可融真金。

咳咳。
每次核融的时候，空都会选取一些原子，排成一列。然后，她会将原子序列分成一些段，并将每段进行一次核融。
一个原子有两个属性：质子数和中子数。
每一段需要满足以下条件：
1、同种元素会发生相互排斥，因此，同一段中不能存在两个质子数相同的原子。
2、核融时，空需要对一段原子加以防护，防护罩的数值等于这段中最大的中子数。换句话说，如果这段原子的中子数最大为x，那么空需要付出x的代价建立防护罩。求核融整个原子序列的最小代价和。

Input

第一行一个正整数N，表示原子的个数。
接下来N行，每行两个正整数pi和ni，表示第i个原子的质子数和中子数。

Output

输出一行一个整数，表示最小代价和。

Sample Input

5
3 11
2 13
1 12
2 9
3 13

Sample Output

26

Data Constraint

对于20%的数据，1<=n<=100
对于40%的数据，1<=n<=1000
对于100%的数据，1<=n<=10^5，1<=pi<=n，1<=ni<=2*10^4

Solution

这题我做过原题，而且这题还更简单，然而我居然没看出来…………
有一道题叫做军训，题解请点进去看

Code

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
#define fo(i,a,b) for(ll i=a;i<=b;i++)
#define N 101000
#define INF 214748367
using namespace std;
struct node{
    ll f,a,l,b;
};
node t[N*20];
int bz[N],n,up[N],a[N],k;
ll f[N],ans;
void get(ll v)
{
    t[v].a=max(t[v*2+1].a,t[v*2].a);t[v].f=min(t[v*2+1].f,t[v*2].f);
    ll x=t[v*2].b,y=t[v*2+1].b;
    x=x==0?INF:x;y=y==0?INF:y;
    t[v].b=min(x,y);
}
void down(ll v)
{
    if(t[v].l==0) return;
    t[v*2].a=max(t[v*2].a,t[v].l);t[v*2+1].a=max(t[v*2+1].a,t[v].l);
    t[v*2].l=max(t[v*2].l,t[v].l);t[v*2+1].l=max(t[v*2+1].l,t[v].l);
    t[v*2].b=t[v*2].a+t[v*2].f;t[v*2+1].b=t[v*2+1].a+t[v*2+1].f;
    t[v].l=0;
}
int fing(ll v,ll l,ll r,ll x)
{
    if(l==r) return l;
    ll m=(l+r)/2;down(v);
    if(t[v*2+1].a>x) return fing(v*2+1,m+1,r,x);
    else if(t[v*2].a>x) return fing(v*2,l,m,x);
    return 0;
}
void finh(ll v,ll l,ll r,ll x,ll y,ll z)
{
    if(l==x&&r==y)
    {
        if(t[v].a>=z)
        {
            int kk=fing(v,l,r,z);
            k=max(k,kk);
        }       
        return;
    }
    ll m=(l+r)/2;down(v);
    if(y<=m) finh(v*2,l,m,x,y,z);
    else if(x>m) finh(v*2+1,m+1,r,x,y,z);
         else finh(v*2,l,m,x,m,z),finh(v*2+1,m+1,r,m+1,y,z);
}
void insert(ll v,ll l,ll r,ll x,ll y,ll z)
{
    if(l==x&&r==y) {t[v].a=z;t[v].l=z;t[v].b=t[v].f+t[v].a;return;}
    ll m=(l+r)/2;down(v);
    if(y<=m) insert(v*2,l,m,x,y,z);
    else if(x>m) insert(v*2+1,m+1,r,x,y,z);
         else insert(v*2,l,m,x,m,z),insert(v*2+1,m+1,r,m+1,y,z);
    get(v);
}
void find(ll v,ll l,ll r,ll x,ll y)
{
    if(l==x&&r==y) {ans=min(ans,t[v].b);return;}
    ll m=(l+r)/2;down(v);
    if(y<=m) find(v*2,l,m,x,y);
    else if(x>m) find(v*2+1,m+1,r,x,y);
         else find(v*2,l,m,x,m),find(v*2+1,m+1,r,m+1,y);
}
void change(ll v,ll l,ll r,ll x,ll z)
{
    if(l==r) {t[v].f=z;t[v].b=t[v].f+t[v].a;return;}
    ll m=(l+r)/2;down(v);
    if(x<=m) change(v*2,l,m,x,z);
    else change(v*2+1,m+1,r,x,z);
    get(v);
}
int main()
{
    freopen("array.in","r",stdin);freopen("array.out","w",stdout);
    scanf("%lld",&n);
    fo(i,1,n)
    {
        ll x;scanf("%lld %lld",&x,&a[i]);
        up[i]=max(up[i-1],bz[x]+1);bz[x]=i;
    }
    fo(i,1,n)
    {
        k=0;finh(1,1,n,up[i],i,a[i]);
        if(k+1<=i)insert(1,1,n,max(k+1,up[i]),i,a[i]);
        ans=INF;find(1,1,n,up[i],i);f[i]=ans;
        change(1,1,n,i+1,f[i]);
    }
    printf("%lld",f[n]);
    fclose(stdin);fclose(stdout);
}