NOIopenjudge 407:Heritage

总时间限制: 

5000ms

内存限制: 

65536kB

描述

Your rich uncle died recently, and the heritage needs to be divided among your relatives and the church (your uncle insisted in his will that the church must get something). There are N relatives (N <= 18) that were mentioned in the will. They are sorted in descending order according to their importance (the first one is the most important). Since you are the computer scientist in the family, your relatives asked you to help them. They need help, because there are some blanks in the will left to be filled. Here is how the will looks:

Relative #1 will get 1 / ... of the whole heritage, 
Relative #2 will get 1 / ... of the whole heritage,
---------------------- ...
Relative #n will get 1 / ... of the whole heritage.

The logical desire of the relatives is to fill the blanks in such way that the uncle's will is preserved (i.e the fractions are non-ascending and the church gets something) and the amount of heritage left for the church is minimized.

输入

The only line of input contains the single integer N (1 <= N <= 18).

输出

Output the numbers that the blanks need to be filled (on separate lines), so that the heritage left for the church is minimized.

样例输入

2

样例输出

2
3

来源

ural 1108

贪心的来想如果只有一个人那就是1/2 ， 两个个人就是1/2，1/3 ，还剩下1/6，如果还要分给下一个人，就会1/7 ，还剩1/42，同理下一个人得1/43，然后看出规律了吧，

然后我们就涉及到了高精乘，然而普通高精乘是过不了的，所以可以用fft，或者ntt这种黑科技（奇技淫巧）。

以下是我ntt

 #include <iostream>
 #include <cstring>
 #include <cstdio>
 using namespace std ;
 typedef long long LL ;
 const int  P = (  <<  ) +  , N =  <<  , G =  , NUM =  ;
 int m , len , l1 , l2 ;
 LL a[N] , b[N] , wn[NUM] ; 

 void Init( )
 {
     scanf( "%d" , &m ) ;
     if( m >  ) l1 =  , l2 =  , a[] =  , b[] =  , len =  ;
 } 

 LL quick_mod( LL a , LL b , LL m )
 {
     LL ans = 1ll ;
     a %= m ;
     while( b )
     {
         if( b &  ) ans = ans * a % m ;
         b >>=  ;
         a = a * a % m ;
     }
     return ans ;
 } 

 void Rader( LL a[] , int n )
 {
     int j = n >>  ;
     for( int i =  ; i < n -  ; ++i )
     {
         if( i < j ) swap( a[i] , a[j] ) ;
         int k = n >>  ;
         while( j >=k )
         {
             j -= k ;
             k >>=  ;
         }
         if( j < k ) j += k ;
     }    

 }

 void NTT( LL a[] , int n , int on )
 {
     Rader( a , n ) ;
     int id =  ;
     for( int h =  ; h <= n ; h <<=  )
     {
         ++id ;
         for( int j =  ; j < n ; j += h )
         {
             LL w =  ;
             for( int k = j ; k < j + ( h >>  ) ; ++k )
             {
                 LL u = a[k] % P ;
                 LL t = w * a[k+(h>>)] % P ;
                 a[k] = ( u + t ) % P ;
                 a[k+(h>>)] = ( u - t + P ) % P ;
                 w = w * wn[id] % P ;
             }
         }
     }
     if( on == - )
     {
         for( int i =  ; i <  ( n >>  ) ; ++i ) swap( a[i] , a[n-i] ) ;
         LL inv = quick_mod( n , P -  , P ) ;
         for( int i =  ; i < n ; ++i  ) a[i] = a[i] * inv % P ;
     }

 }

 void Conv( LL a[] , LL b[] , int n )
 {
     NTT( a , n ,  ) ;
     NTT( b , n ,  ) ;
     for( int i =  ; i < n ; ++i )
         a[i] = a[i] * b[i] % P ;
     NTT( a , n , - ) ;
 }

 void Transfer( LL a[] , int n )
 {
     for( int i =  ; i < n ; ++i )
     {
         a[i+] += a[i]/ ;
         a[i] %=  ;
         if( a[i] ) l1 = i ;
     } ++l1 ;
 }

 void Update( LL b[] , int n )
 {
     int t =  ; l2 = n ;
     for( int i =  ; i <= n ; ++i )
     {
         b[i] = a[i] + t ;
         t = b[i]/ ; b[i] %=  ;
     }
     if( b[l1] ) ++l2 ; 

 }

 void print( LL b[] , int n )
 {
     for( int i = l2 -  ; i >=  ; --i ) putchar( b[i] + '' ) ;
     puts( "" ) ;
 }

 void Length( int l1 , int l2 )
 {
     while( len <= (l1<<) || len <= (l2<<) ) len <<=  ;
 }

 void Getwn( )
 {
     for( int i =  ; i < NUM ; ++i ) wn[i] = quick_mod( G , (P-)/(<<i) , P ) ;
 }

 void Solve( )
 {
     if( m >=  ) puts( "" ) ;
     if( m >=  ) puts( "" ) ;
     if( m <=  ) return ;
     Getwn( ) ;
     for( int i =  ; i <= m ; ++i )
     {
         for( int i = l2 ; i <= len ; ++i ) b[i] =  ;
         Conv( a , b , len ) ;
         Transfer( a , len ) ;
         Update( b , l1 ) ;
         print( b , l2 ) ;
         Length( l1 , l2 ) ;
     }

 }

 int main( )
 {
     Init( ) ;
     Solve( ) ;
     return  ;
 }