【2015.12NOIPtj模拟】神奇的项链(fett)

题目描述

从前有一条神奇的项链，为什么说它神奇呢？因为它有两个性质：
1. 神奇的项链可以拉成一条线，线上依次是N 个珠子，每个珠子有一个能量值Ei；
2. 除了第一个和最后一个珠子，其他珠子都满足E i =(E i-1 +E i+1 )/2+D i 。

由于这条项链很长，我们只能知道其两端珠子的能量值。并且我们知道每个珠子的Di是多少。请聪明的你求出这N 个珠子的能量值分别是多少。

输入

第一行三个整数N、E 1、E N，表示珠子个数N，第一个珠子和第N 个珠子的能量值。
第二行N-2 个整数，表示第2 个珠子到第N-1 个珠子的D i。

输出

输出仅一行，N 个整数，表示1 到N 个这N 个珠子各自的能量值E i。
请放心，数据保证对于任意珠子满足(E i-1+E i+1)Mod 2=0

样例输入

Sample Input 1:
4 1 4
0 0


Sample Input 2:
10 1 22
1 2 -3 5 1 4 2 -1

样例输出

Sample Output 1:
1 2 3 4


Sample Output 2:
1 14 25 32 45 48 49 42 31 22

数据范围限制

40%的数据 1<N<=100。
70%的数据 1<N<=5,000，所有数据（包括计算中的）不超过10^9。
100%的数据 1<N<=500,000，|E1|、|EN|<=10^14，|Di|<=10^4。

题解

这道题一看数据范围，蒙了！

比赛时冥思苦想想不出什么好方法，只能暴力枚举第二个点，逆推求出其他的点，再判断是否合法。这种算法时间复杂度是O（n*max（max为最大点）），显然很慢，最后神奇地拿到了30分。

var
        n,i,j:longint;
        e,d:array[1..500000]of longint;
begin
        assign(input,'fett.in');reset(input);
        assign(output,'fett.out');rewrite(output);
        readln(n,e[1],e[n]);
        for i:=2 to n-1 do read(d[i]);
        repeat
                inc(e[2]);
                for i:=3 to n-1 do e[i]:=(e[i-1]-d[i-1])*2-e[i-2];
                if (e[n]+e[n-2])div 2+d[n-1]=e[n-1] then
                begin
                        for i:=1 to n do write(e[i],' ');
                        writeln;
                        break;
                end;
        until i=0;
        close(input);close(output);
end.

正解

比赛之后终于想通了。这道题其实是很经典的二分。l和r一开始分别赋为-maxlongint和maxlongint，然后取mid，即第二个点的能量值。最后看看啊a[n]是否等于输入的a[n]。如果小了l变成mid，大了r变成mid。轻松拿到100分。

var
        a,d:array[1..500000]of int64;
        i:longint;
        n,l,r,mid,x:int64;
begin
        assign(input,'fett.in');reset(input);
        assign(output,'fett.out');rewrite(output);
        readln(n,a[1],x);
        for i:=2 to n-1 do read(d[i]);
        l:=-maxlongint;
        r:=maxlongint;
        repeat
                mid:=(l+r)div 2;
                a[2]:=mid;
                for i:=2 to n-1 do a[i+1]:=(a[i]-d[i])*2-a[i-1];
                if a[n]=x then
                begin
                        for i:=1 to n do write(a[i],' ');
                        writeln;
                        break;
                end
                else
                if a[n]<x then l:=mid
                else r:=mid;
        until (l>=r)or(r-l=1);
        close(input);close(output);
end.

友情提示：注意循环退出条件，谨慎陷入死循环泥潭。