JZOJ 4937 【GDKOI2017模拟1.12】与运算

与运算

题目大意

给出一个长度为$n$的序列$a_1$、$a_2$、$......$$a_n$。$f_i$表示序列前$i$项依次进行与运算后的结果。
定义一个序列的价值为$\sum_{i=1}^{n}$$f_i$，现在让你重新排列这个序列，使得这个序列的价值最大，并求出序列最大价值。

数据范围

$a_i$<=$10^6$ , $n$<=$10^6$

题解

一个显然的结论，$f_i$会随着i的递增而不上升。
设$g_i$表示满足$a_j$ $\bigvee$ $i$=$a_j$的$a_j$的个数。
设$h_i$表示当前找到的与运算前缀和的末位为$i$的答案贡献，不难发现转移式为
$h_i$=$max$($h_i$,$h_j$+($g_i$-$g_j$)*$i$)其中$i$ $\bigvee$ $j$=$i$
也就是枚举上一次找到的不同的前缀和为j，则中间的那一段的前缀和都为i，中间相差的序列长度为$g_i$-$g_j$，因此中间那段序列的答案贡献为($g_i$-$g_j$)*$i$。
至于求$g_i$可以用分治实现，详见代码。

Code(Pascal)

var
    a,num,f:array[0..1600000] of int64;
    n,m,j,k,l,i,o:longint;
    u,ans:int64;
    m2:array[0..30] of int64;
function max(a,b:int64):int64;
begin  if a>b then exit(a) else exit(b); end;

procedure dg(l,r:longint);//所谓的g数组预处理
    var
        i,mid:longint;
    begin
        if l=r then exit;
        mid:=(l+r) div 2;
        for i:=l to mid do
        num[i]:=num[mid-l+i+1]+num[i];
        dg(l,mid);  dg(mid+1,r);
    end;
begin
    readln(n);
    for i:=1 to n do
    begin
        read(a[i]);
        inc(num[a[i]]);
        o:=max(o,a[i]);
    end;
    o:=trunc(ln(o)/ln(2))+1;
    m2[1]:=1; for i:=2 to o+1 do m2[i]:=m2[i-1]*2;
    u:=m2[o+1];
    dg(0,u-1);
    f[u-1]:=(u-1)*num[u-1]; ans:=f[u-1];
    for i:=u-1 downto 1 do
    begin
        ans:=max(ans,f[i]);
        for l:=1 to o do
        if i and m2[l]>0 then
        f[i-m2[l]]:=max(f[i-m2[l]],f[i]+(num[i-m2[l]]-num[i])*(i-m2[l]));
    end;
    ans:=max(ans,f[0]);
    writeln(ans);
end.