「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=,f=;char ch=getchar();

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

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

     return x*f;

 }
 struct rec{ll x,y;}b[maxn];
 ll n,a[maxn],s[][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<<)-]=read(),a[i<<]=read();
     for1(i,n<<)b[i].x=a[i],b[i].y=i;
     sort(b+,b+*n+,cmp);
     int tot=;
     for1(i,n<<)
     {
         if(i==||b[i].x!=b[i-].x)tot++;
         a[b[i].y]=tot;c[tot]=b[i].x;
     }
     for1(i,n)
     {
        int l=a[(i<<)-],r=a[i<<];
        s[][]+=c[l];s[][l]-=c[l];
        s[][l]++;s[][r+]--;
     }
     int ans=;
     for1(i,tot+)
     {
         s[][i]+=s[][i-];s[][i]+=s[][i-];
         if(s[][i]+s[][i]*c[i]>s[][ans]+s[][ans]*c[ans])ans=i;
     }
     printf("%lld %lld\n",c[ans],s[][ans]+s[][ans]*c[ans]);
     return ;

 }