荷马史诗（NOI2015）提高组

1. Problem Description

追逐影子的人，自己就是影子 ——荷马

Allison 最近迷上了文学。她喜欢在一个慵懒的午后，细细地品上一杯卡布奇诺，静静地阅读她爱不释手的《荷马史诗》。但是由《奥德赛》和《伊利亚特》 组成的鸿篇巨制《荷马史诗》实在是太长了，Allison 想通过一种编码方式使得它变得短一些。

一部《荷马史诗》中有n种不同的单词，从1到n进行编号。其中第i种单 词出现的总次数为wi。Allison 想要用k进制串si来替换第i种单词，使得其满足如下要求：

对于任意的 1 ≤ i, j ≤ n ， i ≠ j ，都有：si不是sj的前缀。

现在 Allison 想要知道，如何选择si，才能使替换以后得到的新的《荷马史诗》长度最小。在确保总长度最小的情况下，Allison 还想知道最长的si的最短长度是多少？

一个字符串被称为k进制字符串，当且仅当它的每个字符是 0 到 k − 1 之间（包括 0 和 k − 1 ）的整数。

字符串 str1 被称为字符串 str2 的前缀，当且仅当：存在 1 ≤ t ≤ m ，使得str1 = str2[1..t]。其中，m是字符串str2的长度，str2[1..t] 表示str2的前t个字符组成的字符串。

2. Input

输入的第 1 行包含 2 个正整数 n, k ，中间用单个空格隔开，表示共有 n种单词，需要使用k进制字符串进行替换。

接下来n行，第 i + 1 行包含 1 个非负整数wi ，表示第 i 种单词的出现次数。

3. Output

输出包括 2 行。

第 1 行输出 1 个整数，为《荷马史诗》经过重新编码以后的最短长度。

第 2 行输出 1 个整数，为保证最短总长度的情况下，最长字符串 si 的最短长度。

4.详解

#include <iostream>
#include <cstdio>
#define ll long long

using namespace std;

int k;
const int M=200001;
struct node{
    ll id,data;
    ll lenth;
}temp,s[M],x;
int r,l;

ll n,all;
ll asize,size,w[100001],fa[200001][10];
bool cmp(node a,node b)
{ return a.data<b.data||(a.data==b.data&&a.lenth<b.lenth);}

void swapx(int i, int j){
    node tmp;
    tmp=s[i]; s[i]=s[j]; s[j]=tmp;
}
/*void out (){
    for(int i=0;i<=r;i++){
        printf("%lld ",s[i].data);
    }
    cout<<endl;
}*///调试
void insert(node x){    //将数x插入堆中
    r++;s[r]=x; l=r;
    while ((l>1)&&(cmp(s[l],s[l>>1]))){
        swapx(l,l>>1); l=l>>1;
    }
}
void change(node x){  //调整小根堆
    int rc,lc;
    node min=x; lc=2*x.id; rc=2*x.id+1;
    if ((lc<=r)&&(cmp(s[lc],s[min.id]))) min.id=lc;
    if ((rc<=r)&&(cmp(s[rc],s[min.id]))) min.id=rc;
    if (min.id!=x.id){
        swapx(x.id,min.id); change(min);
    }
}
node get(){    // 取堆顶元素并从堆中删除
    node res,a1;
    a1.data=1,a1.id=1,a1.lenth=1;
    res=s[1];  s[1]=s[r]; r--;
    if (r>1) change(a1);
    return res;
}
void init(){
    cin>>n>>k;
    temp.lenth=0;
    for(int i=1;i<=n;i++){
        cin>>w[i];
        fa[i][k]=1;
        temp.id=i;
        temp.data=w[i];
        l++;
        insert(temp);
        //out();
    }
    if(k!=2){
        while(n%(k-1)!=1){
            n++;
            temp.id=n;
            temp.data=0;
            fa[n][k]=1;
            insert(temp);
        }
    }
    all=n+1;
    while(r){
        temp.data=0;
        temp.lenth=0;
        temp.id=all++;
        for(int i=0;i<k;i++){
            temp.lenth=max((s[1].lenth)+1,temp.lenth);
            temp.data+=(s[1]).data;
            fa[temp.id][i]=(s[1]).id;
            get();
        }
        if(!r)break;
        fa[temp.id][k]=1;
        insert(temp);
    }
}
void dfs(ll x){
    if(x<=n)return;

    for(int i=0;i<k;i++){
        fa[fa[x][i]][k]+=fa[x][k];
        dfs(fa[x][i]);
    }
}
int main(){
    init();
    dfs(all-1);
    for(int i=1;i<=n;i++){
        asize+=fa[i][k]*w[i];
        size=max(size,fa[i][k]);
    }
    cout<<asize<<endl<<size;
}