zoj 2711 Get Out of the Glass

/*
zoj_2771    dp+高精度
簡單dp，不過又是要結合上高精度。。

思路: 對玻璃四層（即四個交接面）分別標記爲0,1,2,3
      dp[i][j] 表示第j個反射點在第i層的可能數。
      則易知：
        j%2==0時有： dp[i][j]=sum( dp[k][j-1] ); i<k<4
              否則： dp[i][j]=sum( dp[k][j-1] ); 0=<k<i
*/
#include <iostream>
#include <cstdio>
#include <string.h>
using namespace std;
int dp[4][61][21];

void print( int n )
{
    int i,j,y;
    int sum[21]={0};
    for( i=0;i<4;i++ )    //sum+=dp[i][n]
    {
        y=0;
        for( j=20;j>=0;j-- )
        {
            sum[j]+=( dp[i][n][j]+y );
            y=sum[j]/10000;
            sum[j]%=10000;
        }
    }
    for( i=0;i<21;i++ )
        if( sum[i]!=0 )
        {
            cout<<sum[i];   i++;
            break;
        }
    for( ;i<21;i++ )
    {
        if( sum[i]<1000 ) cout<<'0';
        if( sum[i]<100 )  cout<<'0';
        if( sum[i]<10 )   cout<<'0';
        cout<<sum[i];
    }
    cout<<endl;
}

int main()
{
    int n,i,j,k,l,y;
    memset( dp,0,sizeof(dp) );
    for( i=1;i<4;i++ )  dp[i][1][20]=1;
    for( j=2;j<61;j++ )
    {
        for( i=0;i<4;i++ )
        {
            if( j%2==0 )
                for( k=i+1;k<4;k++ ) //dp[i][j]+=dp[k][j-1]
                {
                    y=0;
                    for( l=20;l>=0;l-- )
                    {
                        dp[i][j][l]+=( dp[k][j-1][l]+y );
                        y=dp[i][j][l]/10000;
                        dp[i][j][l]%=10000;
                    }
                }
            else
                for( k=0;k<i;k++ ) //dp[i][j]+=dp[k][j-1]
                {
                    y=0;
                    for( l=20;l>=0;l-- )
                    {
                        dp[i][j][l]+=( dp[k][j-1][l]+y );
                        y=dp[i][j][l]/10000;
                        dp[i][j][l]%=10000;
                    }
                }
        }
    }
    while( cin>>n )
    {
        if( n==0 )  cout<<"1\n";
        else print(n);
    }
    return 0;
}