{题解}[jzoj3396] 【NOIP2013模拟】Rainbow的信号

jzoj3396

Description
   Freda发明了传呼机之后, rainbow 进一步改了传呼机发送信息所使用的号。 由于现在是数字、信息时代, rainbow 发明的信号用 N个自然数表示。 为了避免两个人的对话被 大坏蛋 VariantF偷听 T_T,rainbow 把对话分成 对话分成A、B、C三部分 ,分别用  a、b、c三个密码 加密 。现在 Freda接到了 rainbow的信息,她首要工作就是解密。 Freda了解到,这三部 分的密码计算方式如下:
    在 1~N这 N个数 中,等概率地选取两个数l、r,如果 l>r,则交换 l、r。把信号中的第 l个数到第 r个数 取出来,构成一数列 P。
    A部分对话的密码是数列 P的 xor 和的数学 期望值 。xor 和就是 数列 P中各个数异或之 后得到的数 ;xor和的期望 就是对于所有可能选取 的 l、r,所得到的 数列xor和的平均数 。
    A部分对话 占接收到的信息总量 的 40% ,因此 如果你计算出密码 a,将获得该测试点40% 的分数。
   B部分对话的密码是数列 P的 and 和的期望 ,定义类似于 xor和,占信息总 量的 30% 。
    C部分对话的密码是数列 P的 or 和的期望 ,定义类似于 xor 和,占信息总量的 30% 。

40分做法:
直接数位DP
枚举每个数列求值,
更新答案。

100分做法:
在40分基础上多想想,
我们可以再套一个DP
设f[i],g[i],d[i]分别为以第i个数结尾的区间的xor和,and和,or和。
(暂时不考虑期望)
想想。似乎不能直接考虑,需要分位考虑。
再想想。咦!好像和40分做法并没有什么关系呀!
的确没有什么关系。
方程来了!!!
设n1为最近的1(之前),n0为最近的0(之前)
于是就有
{题解}[jzoj3396] 【NOIP2013模拟】Rainbow的信号_第1张图片
多么美妙啊!!!
都很好推,多多思考吧!
Tips
小心精度问题

Code

const
        maxn=100000;
        maxnum=1000000000var
        n:longint;
        i:longint;
        ansa,ansb,ansc:real;
        a:array[0..maxn] of longint;
        f,g,d:array[-1..maxn] of int64;
        x:array[0..maxn,0..32] of byte;
procedure init;
var
        i,l:longint;
begin
        readln(n);
        for i:=1 to n do
        begin
                read(a[i]);
                l:=a[i];
                while(l>0)do
                begin
                        inc(x[i,0]);
                        x[i,x[i,0]]:=l and 1;
                        l:=l shr 1;
                end;
        end;
end;
procedure work(o:longint);
var
        i:longint;
        near0,near1:longint;
begin
        near0:=0;near1:=0;
        for i:=1 to n do
        begin
                if x[i,o]=0 then
                begin
                        f[i]:=f[i-1];
                        g[i]:=0;
                        d[i]:=near1;
                        near0:=i;
                end
                else
                if x[i,o]=1 then
                begin
                        f[i]:=f[near1-1]+(i-near1);
                        g[i]:=i-near0;
                        d[i]:=i;
                        near1:=i;
                end;
                ansa:=ansa+(f[i] shl o);
                ansb:=ansb+(g[i] shl o);
                ansc:=ansc+(d[i] shl o);
        end;
end;
begin
        init;
        for i:=1 to 32 do
                work(i);
        for i:=1 to n do
        begin
                ansa:=ansa-a[i];
                ansb:=ansb-a[i];
                ansc:=ansc-a[i];
        end;
        writeln(ansa/n/n:0:3,' ',ansb/n/n:0:3,' ',ansc/n/n:0:3);
end.

你可能感兴趣的:({题解}[jzoj3396] 【NOIP2013模拟】Rainbow的信号)