[BZOJ4576]262144 分治/区间DP

只有我是鬼畜的分治做法吗。。。

先形象的理解，假如某一个很小的数孤零零的在那里，它阻断了左右两边的合并，导致问题可以被分治解决

考虑从小的数字往大的数字合并，若有一2*k长度的相同数字块，我们可以直接将数字加一而长度变为k，若有一2*k+1长度的，则必定会阻断左右，则要分合并后的k个分给左边和右边两种情况讨论。

solve(l,r,b)表示把l到r的中的b元素全部合并所能的到的最大值。用双向链表维护上述数字块的前后关系，每合并一种数字b会导致一些块被阻断（奇数个时），分治即可。注意当奇数个且大于1时，要有两个节点（都表示k个）（就再新建一个），来分给两个子问题。

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=270000;
int n,pre[maxn],suc[maxn];
struct node
{
    int t,num;
}a[maxn];
void update(int l,int r)
{
    for(int i=l;i<=r;i=suc[i])
    {
        if(i!=l&&a[pre[i]].t==a[i].t) 
                    {
                        a[pre[i]].num+=a[i].num;
                        suc[pre[i]]=suc[i];
                        pre[suc[i]]=pre[i];
                    }
    }
}
int solve(int l,int r,int b)
{
    if(l>r||b>60) return 0;
    update(l,r);
    int last=l,re=0,minb=(maxn<<1);
    for(int i=l;i<=r;i=suc[i])
    {
        re=max(re,a[i].t);
        if(a[i].t==b) 
            {
                if(a[i].num==1)
                {
                    re=max(re,solve(last,pre[i],minb));minb=(maxn<<1);
                    last=suc[i];
                    continue;
                }
                int pd=(a[i].num&1);
                a[i].t++;a[i].num>>=1;minb=b+1;
                if(pd)
                {
                    int tmp=a[i].num;
                    suc[i+tmp]=suc[i];pre[i+tmp]=i;
                    suc[i]=i+tmp;pre[suc[i+tmp]]=i+tmp;
                    a[i+tmp]=a[i];
                    re=max(re,solve(last,i,minb));minb=(maxn<<1);
                    last=i+tmp;
                }
            }
            else minb=min(minb,a[i].t);
    }

    if(re<b) return re;
    re=max(re,solve(last,r,minb));
    return re;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i].t);
        a[i].num=1;
        suc[i]=i+1;
        pre[i]=i-1;
    }
    printf("%d",solve(1,n,1));
    return 0;
}

当然还有正解的区间DP做法（搬自洛谷题解）
f[i][j]表示从j开始合并到i这个数字序列的末尾的下表是什么。
那么因为合并的总是一段连续的区间，就有f[i][j]=f[i-1][f[i-1][j]];
代码：

#include <cstdio>
#define max(a,b) a>b?a:b
int x,n,ans=0;
int f[60][270000];
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;++i)
        scanf("%d",&x),f[x][i]=i+1;
    for (int i=2;i<=58;++i)
        for (int j=1;j<=n;++j)
        {
            if (!f[i][j]) f[i][j]=f[i-1][f[i-1][j]];
            if (f[i][j]) ans=i;
        }
    printf("%d\n",ans);
}