1653: [Usaco2006 Feb]Backward Digit Sums

1653: [Usaco2006 Feb]Backward Digit Sums

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 285  Solved: 215
[Submit][Status]

Description

FJ and his cows enjoy playing a mental game. They write down the numbers from 1 to N (1 <= N <= 10) in a certain order and then sum adjacent numbers to produce a new list with one fewer number. They repeat this until only a single number is left. For example, one instance of the game (when N=4) might go like this: 3 1 2 4 4 3 6 7 9 16 Behind FJ's back, the cows have started playing a more difficult game, in which they try to determine the starting sequence from only the final total and the number N. Unfortunately, the game is a bit above FJ's mental arithmetic capabilities. Write a program to help FJ play the game and keep up with the cows.

Input

* Line 1: Two space-separated integers: N and the final sum.

Output

* Line 1: An ordering of the integers 1..N that leads to the given sum. If there are multiple solutions, choose the one that is lexicographically least, i.e., that puts smaller numbers first.

Sample Input

4 16

Sample Output

3 1 2 4

OUTPUT DETAILS:

There are other possible sequences, such as 3 2 1 4, but 3 1 2 4
is the lexicographically smallest.

HINT

 

Source

Silver

 题解：这个嘛，本来还想什么高端洋气的算法的，可是再想想果断决定——弃疗——10！不过才3628800而已嘛，（具体说算法嘛，很显然对于原数列，每个数依次的最终累计次数即是杨辉三角形第N行的对应数字，别的没了），可是再一想，有点不对——你不是每次还要判断此组解是否合法么？这样子复杂度可还要再×10哦（36288000，3kW多了，这下子可危险啊，虽然事实上只要有解的话，由于杨辉三角形的对称性，那么最多理论上一半的时间即可找到解）。。。可是结果是——228kb 60ms Accept我也是醉了。。。

var
   i,j,k,l,m,n:longint;
   a:array[0..20] of longint;
   b:array[0..20,0..20] of longint;
procedure swap(var x,Y:longint);
          var z:longint;
          begin
               z:=x;x:=y;y:=z;
          end;
procedure sort(l,r:longint);
          var i,j,x,y:longint;
          begin
               i:=l;j:=r;x:=a[(l+r) div 2];
               repeat
                     while a[i]<x do inc(i);
                     while a[j]>x do dec(j);
                     if i<=j then
                        begin
                             swap(a[i],a[j]);
                             inc(i);dec(j);
                        end;
               until i>j;
               if l<j then sort(l,j);
               if i<r then sort(i,r);
          end;
begin
     readln(n,m);
     for i:=0 to n do a[i]:=i;
     fillchar(b,sizeof(b),0);
     b[1,1]:=1;
     for i:=2 to n do
         for j:=1 to i do b[i,j]:=b[i-1,j-1]+b[i-1,j];
     while a[0]=0 do //萌萌哒生成法全排列，小学时学的现在居然还记得*_*
           begin
                l:=0;
                for i:=1 to n do l:=l+a[i]*b[n,i];
                if l=m then
                   begin
                        for i:=1 to n-1 do
                            write(a[i],' ');
                        writeln(a[n]);
                        halt;
                   end;
                j:=n;
                while a[j-1]>a[j] do dec(j);
                k:=n;
                while a[j-1]>a[k] do dec(k);
                swap(a[j-1],a[k]);
                sort(j,n);
           end;
end.