1

我正在尝试为帕斯卡三角形制作一个程序,计算行中rth值的公式nthn!/r!(n-r)!

我一直在尝试像这样实现它:

          #include <iostream>     // std::cout, std::endl
          #include <iomanip>      // std::setw
          int Pascal(int ,int );
          int Factorial( int);
          int main () 
          {
           int n=0;//Number of rows in the triangle
           for(int i=12;i>0;i--)
            {
              std::cout << std::setw(i)<<std::endl;
              for (int j=1;j<12-i;j++)
               {
                int r=0; //rth element initialized for each row
                int P= Pascal(r,n);
                std::cout << P ;
                std::cout <<std::setw(2);
                r=r+1;
               }

              n=n+1;

           }
               std::cout<<std::endl;
                   return 0;
           }
           int Pascal(int r,int n)
           {
               int d = n-r; ///Difference of n with r
               int f1; ///factorial of n
               int f2; ///factorial of r
               int f3; ///factorial of (n-r)

               f1=Factorial(n);
               f2=Factorial(r);
               f3=Factorial(d);
               return f1/(f2*f3);

           }
           int Factorial( int begin )
           {
                   int F;
                   if ( begin == 0 )
                   {
                    return 1; 
                   }
                   else
                   {
                     F= begin*Factorial(begin-1);
                   }
              return F;         
             }

但不知何故,我得到了这个输出:

                   1
                  1 1
                 1 1 1
                1 1 1 1
               1 1 1 1 1
              1 1 1 1 1 1
             1 1 1 1 1 1 1
            1 1 1 1 1 1 1 1
           1 1 1 1 1 1 1 1 1
          1 1 1 1 1 1 1 1 1 112

有人可以指导我哪里出错了吗?

4

2 回答 2

15

您的第一个明显问题当然是您应该打印 Pascal(i,j)。你的第二个更微妙:

递归

帕斯卡三角形的全部意义在于它提供了一种计算二项式系数的方法,而无需计算阶乘。

问题是阶乘增长得非常快,而且像 Pascal(1,120) = 120!/(1!*119!} 这样的系数只等于 120,但它的分母和分母都在 10^198 的数量级上,这不可能以任何机器整数类型存储。

查看维基百科上的帕斯卡三角。重点是递归关系:

Pascal( r, n ) = Pascal( r-1, n-1 ) + Pascal( r, n-1 )

这是一个使用它的简单解决方案(这只打印 r,n 帕斯卡数):

#include <iostream>
#include <sstream>

int pascal( int r, int n ) {
    if( n == 0 )
        return 1;

    if( r == 0 || r == n )
        return 1;

    return pascal( r - 1, n - 1 ) + pascal( r, n - 1 );
}

int main( int argc, char *argv[] ) {
    if( argc != 3 ) {
        std::cout << "Expected exactly 3 arguments." << std::endl;
        return -1;
    }

    int r, n;

    std::stringstream ss;
    ss << argv[1] << ' ' << argv[2];
    ss >> r >> n;

    if( ss.bad() || r < 0 || n < 0 || r > n ) {
        std::cout << "Invalid argument values." << std::endl;
        return -2;
    }

    std::cout << pascal( r, n ) << std::endl;
    return 0;
}

记忆

这种方法还有第三个甚至更微妙的问题,它的复杂性比它需要的要糟糕得多。考虑计算 Pascal(3,6):

                       Pascal(3,6)                      =
=        Pascal(2,5)        +        Pascal(3,5)        =
= (Pascal(1,4)+Pascal(2,4)) + (Pascal(2,4)+Pascal(3,4)) = ...

在最后一行中,您会注意到 Pascal(2,4) 出现两次,这意味着您的代码将计算两次。此外,Pascal(3, 5) 实际上等于 Pascal(2,5)。所以你可以计算 Pascal(2,5) 两次,这意味着计算 Pascal(2,4) 四次。这意味着随着 r 和 n 变大,程序将非常缓慢。我们想计算每个 Pascal(i,j) 一次,然后保存它的值以供其他调用使用。为此,一种简单的方法是使用将 (r,n) 对映射到 Pascal(r,n) 值的映射:std::map< std::pair, int >。这种方法称为记忆。此外,将 int 更改为 long long 以获得大数,您将获得以下算法:

#include <iostream>
#include <sstream>
#include <map>

typedef long long                                          Integer;
typedef std::map< std::pair< Integer, Integer >, Integer > MemoMap;
typedef MemoMap::iterator                                  MemoIter;

MemoMap memo;

Integer pascal( Integer r, Integer n ) {
    // make sure r <= n/2 using the fact that pascal(r,n)==pascal(n-r,n)
    if( r > n / 2LL )
        r = n - r;

    // base cases
    if( n == 0LL || r == 0LL )
        return 1LL;

    // search our map for a precalculated value of pascal(r,n)
    MemoIter miter = memo.find( std::make_pair( r, n ) );

    // if it exists return the precalculated value
    if( miter != memo.end() )
        return miter->second;

    // otherwise run our function as before
    Integer result = pascal( r - 1LL, n - 1LL ) + pascal( r, n - 1LL );

    // save the value and return it
    memo.insert( std::make_pair( std::make_pair( r, n ), result ) );

    return result;
}

int main( int argc, char *argv[] ) {
    if( argc != 3 ) {
        std::cout << "Expected exactly 3 arguments." << std::endl;
        return -1;
    }

    Integer r, n;

    std::stringstream ss;
    ss << argv[ 1 ] << ' ' << argv[ 2 ];
    ss >> r >> n;

    if( ss.bad() || r < 0LL || n < 0LL || r > n ) {
        std::cout << "Invalid argument values." << std::endl;
        return -2;
    }

    std::cout << pascal( r, n ) << std::endl;

    return 0;
}

以前,Pascal(5,100) 永远不会终止。现在它在我的电脑上立即完成。将此代码集成到您的代码中,您将成为一只快乐的熊猫。

自下而上

记忆是自上而下的动态编程,即您首先想到最难的问题,然后将其划分为更简单的重叠问题,同时保存结果。自下而上的解决方案将从帕斯卡三角形的第一行开始并继续计算后续行 - 这在速度和内存方面都更有效(仅存储两个数组)。然而,自上而下的方法更容易实现(您不必考虑需要什么值,只需保存所有内容),它允许您在独立的 pascal() 调用之间保存中间结果。在您的情况下,由于您试图通过多次独立调用 pascal 来打印 Pascal 的三角形,因此这是合适的方法。如果您同时打印和计算,自下而上是迄今为止最有效的方法:

#include <iostream>
#include <sstream>
#include <vector>


typedef long long Integer;

void print_pascal( Integer maxn ) {
    std::vector< Integer > prevRow, currentRow;

    prevRow.resize( maxn + 1 );
    currentRow.resize( maxn + 1);

    prevRow[ 0 ] = 1;
    // print first row.
    std::cout << 1 << std::endl;
    for( Integer currentN = 1 ; currentN <= maxn ; ++ currentN ) {
        // compute & print current row
        currentRow[ 0 ] = currentRow[ currentN ] = 1;
        std::cout << 1;
        for( Integer r = 1 ; r < currentN ; ++ r ) {
            currentRow[ r ] = prevRow[ r - 1 ] + prevRow[ r ];
            std::cout << ' ' << currentRow[ r ];
        }
        std::cout << ' ' << 1 << std::endl;

        // constant time because swap() only swaps internal ptrs.
        currentRow.swap( prevRow );
    }
}

int main( int argc, char *argv[] ) {
    if( argc != 2 ) {
        std::cout << "Expected exactly 1 argument." << std::endl;
        return -1;
    }

    Integer maxn;

    std::stringstream ss;
    ss << argv[ 1 ]; ss >> maxn;

    if( ss.bad() || maxn < 0LL ) {
        std::cout << "Invalid argument values." << std::endl;
        return -2;
    }

    print_pascal( maxn );

    return 0;
}
于 2013-05-23T10:28:29.680 回答
0

您出错的第一件事是代码的格式。这太可怕了。(对不起,确实如此。)

第二件事是你总是在打印Pascal(r, n),而它们是不变的0, 0——你应该打印Pascal(i, j),因为ij是循环计数器。

顺便说一句,您最好迭代地计算阶乘并使用足够长的整数,您的代码在 IDEOne上抛出 SIGXCPU 并在我的计算机上出现段错误。

于 2013-05-23T09:09:05.850 回答