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