「Poetize10」封印一击

描述 Description

Nescafe由n种元素组成（编号为1~n）， 第i种元素有一个封印区间[ai,bi]。当封印力度E小于ai时，该元素将获得ai的封印能量；当封印力度E在ai到bi之间时，该元素将获得E的封印 能量；而当封印力度E大于bi时，该元素将被破坏从而不能获得任何封印能量。现在圣主applepi想选择恰当的E，使得封印获得的总能量尽可能高。为了 封印的最后一击尽量完美，就请你写个程序帮他计算一下吧！

题解：
首先必须离散化，然后差分可以算出在区间外且<a[i]的总得分，再差分一次可以算出每个点被区间覆盖了几次，然后就是统计了。
代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn 1000000
#define maxm 500+100

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define mod 1000000007

using namespace std;

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}
struct rec{ll x,y;}b[maxn];
ll n,a[maxn],s[3][maxn],c[maxn];
inline bool cmp(rec a,rec b){return a.x<b.x;}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);

    n=read();
    for1(i,n)a[(i<<1)-1]=read(),a[i<<1]=read();
    for1(i,n<<1)b[i].x=a[i],b[i].y=i;
    sort(b+1,b+2*n+1,cmp);
    int tot=0;
    for1(i,n<<1)
    {
        if(i==1||b[i].x!=b[i-1].x)tot++;
        a[b[i].y]=tot;c[tot]=b[i].x;
    }
    for1(i,n)
    {
       int l=a[(i<<1)-1],r=a[i<<1];
       s[1][1]+=c[l];s[1][l]-=c[l];
       s[2][l]++;s[2][r+1]--;
    }
    int ans=0;
    for1(i,tot+1)
    {
        s[1][i]+=s[1][i-1];s[2][i]+=s[2][i-1];
        if(s[1][i]+s[2][i]*c[i]>s[1][ans]+s[2][ans]*c[ans])ans=i;
    }
    printf("%lld %lld
",c[ans],s[1][ans]+s[2][ans]*c[ans]);
    return 0;

}