2015年普及组第四题

来来，题目在这里：https://www.luogu.org/problemnew/show/P2672

至于这道题嘛，说实话个人感觉有点恶心。
还是一样先说60分的做法，枚举距离最远的一家住户的位置，然后在前面选出前i-1个（如果说是要到i家的话）最大的劳累值，然后求出最大的和就好了。值得一提的是，在最后一家住户越来越远的时候，如果这家住户从i移到了i+1的话，无非就是将i的值添加到原来维护的那个数组里面，就是添加一项，然后使用插入排序。但是，如果说记录每一个数在前i个中出现了多少次的话，其实就更快一些，毕竟插入枚举的是n,这种方法枚举的是a[n]的最大值，就这道题的数据来讲，明显是后者较优。至于为什么，能大就取大嘛~~~
具体实现请看代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<utility>
#include<cstring>
using namespace std;
const int maxN=100100;
int n,a[maxN],dis[maxN],ans[maxN];
int num[1010],maxA;//f[x][y]表示以x结尾的第y大的数 

int main()
{
    freopen("salesman.in","r",stdin);
    freopen("salesman.out","w",stdout);
    cin>>n;
    for(int i=1;i<=n;++i)
        cin>>dis[i];
    for(int i=1;i<=n;++i){
        cin>>a[i];
        maxA=max(maxA,a[i]);
    }
    for(int i=1;i<=n;++i){//i是最后一个的位置 
        num[a[i-1]]++;
        int sum=a[i]+2*dis[i];
        ans[1]=max(ans[1],sum);
        int k=maxA;
        int temp=num[k];
        for(int j=1;j<=i-1;++j){//j是访问的个数 
            while(temp==0){
                k--;
                temp=num[k];
            }
            sum+=k;
            temp--; 
            ans[j+1]=max(ans[j+1],sum);
        }
    }
    for(int i=1;i<=n;++i)
        cout<<ans[i]<<endl;
    return 0;
}

接下来是满分的做法：
通过一系列的严谨的证明，我们发现：如果在某种情况下，如果选k这家住户是最优方案的话，那么在之后的方案中都会保留原来的选择。假如说第一次选的是第i个住家，那么剩余的住家分为两个部分：（1）i以前的住家，这些住家不需要考虑距离问题，因为最远是到i住家这个位置。 （2）i之后的住家，这些住家不仅要考虑自己的值，还要考虑距离的问题，应为它比现在最远的远了

由此我们可以将两个范围分成两个大根堆，分别维护最大值，然后比较两边的最大值所能创造出的疲劳值哪个大。如果是左边的大，那么直接在答案上加上这个值。如果是右边的大，那么不光要加上这个值，还要把最后选的这个数之前的数加到左边的堆里面，然后在右边的堆里面去除。由于我们不知道去除的数在队中的位置，所以可以用堆的那个数组来维护序号（因为所有住户的值都是知道的，求一下就好了）

大概思路就是这样的，具体实现还是看代码吧。。。

#include<iostream>
#include<cstdio>
using namespace std;
const int maxN=100100;
const int d[2][2]={{0,1},{2,1}};
int n,a[maxN],value[maxN],ans,last=0;
int len[2],h[2][maxN],p[maxN];

int getvalue(int x,int y)
{
    return d[x][0]*a[h[x][y]]+d[x][1]*value[h[x][y]];
}

void down(int x,int y)
{
    int z=y+y;
    while(z<=len[x]){
        if(z+1<=len[x]&&getvalue(x,z+1)>getvalue(x,z))
            z++;
        if(getvalue(x,z)>getvalue(x,y)){
            swap(h[x][y],h[x][z]);
            p[h[x][y]]=y;
            p[h[x][z]]=z;
            y=z;
            z=y+y;
        }       
        else
            break;
    }
}

void up(int x,int y)
{
    while(y>1&&getvalue(x,y)>getvalue(x,y/2)){
        swap(h[x][y],h[x][y/2]);
        p[h[x][y]]=y;
        p[h[x][y/2]]=y/2;
        y>>=1; 
    }
}
void build(int x)
{
    for(int i=len[1]/2;i>=1;i--){
        down(1,i);
    }
}

void init()
{
    len[1]=n;
    for(int i=1;i<=n;i++){
        h[1][i]=i;
        p[i]=i;
    }
    build(1); 
}


void insert(int x,int y){
    h[x][++len[x]]=y;
    p[y]=len[x];
    up(x,len[x]);
}

void del(int x,int y)
{
    h[x][y]=h[x][len[x]--];
    p[h[x][y]]=y;
    up(x,y);
    down(x,y);
}

void work(int x,int y)
{
    del(1,1);
    for(int i=x;i<=y-1;i++){
        del(1,p[i]);
        insert(0,i);
    }
}

int main()
{
    cin>>n;
    for(int i=1;i<=n;++i){
        cin>>a[i];
    }
    for(int j=1;j<=n;++j){
        cin>>value[j];
    }
    init();
    for(int i=1;i<=n;i++)
    {
        if(len[0]==0){
            ans+=2*(a[h[1][1]]-a[last])+value[h[1][1]];
            int temp=h[1][1];
            work(last+1,h[1][1]);
            last=temp;
            cout<<ans<<endl;
            continue;
        }
        if(len[1]==0){
            ans+=value[h[0][1]];
            del(0,1);
            cout<<ans<<endl;
            continue;
        }
        if(value[h[0][1]]>2*(a[h[1][1]]-a[last])+value[h[1][1]]){
            ans+=value[h[0][1]];
            del(0,1);
        }
        else{
            ans+=2*(a[h[1][1]]-a[last])+value[h[1][1]];
            int temp=h[1][1];
            work(last+1,h[1][1]);
            last=temp;
        }
        cout<<ans<<endl;
    }
    return 0;
}

终于写完了，撒花~~~