기본 콘텐츠로 건너뛰기

11051 : 이항계수2 (Binomial Coefficient) (Dynamic Programming) [C++]

저번 이항 계수의 문제와 비슷하지만

입력 값의 범위가 1000까지 늘어서 DP를 이용하기에 딱 좋은 문제이다.


왜 DP를 이용하는가?


이항계수의 원리를 보면 파스칼의 삼각형 형태!

1
1(1C0) 1(1C1)
1(2C0) 2(2C1) 1(2C2)
1(3C0) 3(3C1) 3(3C2) 1(3C3)
1(4C0) 4(4C1) 6(4C2) 4(4C3) 1(4C4)
1(5C0) 5(5C1) 10(5C2) 10(5C3) 5(5C4) 1(5C5)
1(6C0) 6(6C1) 15(6C2) 20(6C3) 15(6C4) 6(6C5) 1(6C6)


<파스칼의 삼각형(이항원리)>


DP식을 써보면

///////////////////////////////////////////////////////////
if (k == 0 || k == n)
       DP[i][j] = 1;
else
       DP[i][j] = DP[i-1][j-1] + DP[i-1][j];
///////////////////////////////////////////////////////////



BinomialCoefficient::BinomialCoefficient(int N)
{
    for(int i = 0; i < N; i++)
    {
        vector<int> Temp;
        Temp.resize(i+2);
        DP.push_back(Temp);
    }
    
    for(int i = 0; i < N; i++)
    {
        for(int j = 0; j < i + 2; j++)
        {
            if(j == 0 || j == i + 1)
                DP[i][j] = 1;
            else
                DP[i][j] = DP[i-1][j-1] % 10007 + DP[i-1][j] % 10007;
        }
    }
}

void BinomialCoefficient::Print(int N, int K)
{
    cout << DP[N-1][K] % 10007;
}
<소스 코드>

*Source of the problem = https://www.acmicpc.net/problem/11051
*문제 출처 : BAEKJOON ONLINE JUDGE

댓글

이 블로그의 인기 게시물

11401 : 이항계수 3 (Dynamic Programming, Divide and Conquer) [C++]

1번, 2번 문제들과 확연히 차이나는 입력의 범위. 400만 ! DP를 사용해서 풀 수 없는 문제이다 . 하지만 DP가 쓰이긴 한다! 수학은 너무 어렵다. 곱셈의 역원을 공부해보다가 모르겠어서 도움을 구했다. 왜 곱셈의 역원을 구해야하는가? N! / (K! * (N-K)!) 에서 K! * (N-K)! 의 역원을 구해야 하기 때문! 곱셈의 역원을 구하는 정리인 페르마의 소정리를 이용하면 p가 1000000007 이지만 분할 정복을 이용한 제곱 수 계산 덕분에 logP 시간 소요. 분할 정복을 이용한 제곱 수 계산은 계속 써먹을 것 같아서 따로 올려놓았다. DP가 쓰이는 부분은 구해준 400만의 역원을 바탕으로 모든 역원을 구하는 부분이다. 그러므로 총 시간 소요는 O(N+LogP) long long BinomialCoefficient :: GetNum ( int N , int K ) { this - > Factorial [ 1 ] = 1 ; for ( int i = 2 ; i < = 4000000 ; i + + ) this - > Factorial [ i ] = ( this - > Factorial [ i - 1 ] * i ) % P ; this - > Invert [ 4000000 ] = this - > Pow_DC ( this - > Factorial [ 4000000 ] , P - 2 ) ; for ( int i = 4000000 - 1 ; i > 0 ; i - - ) this - > Invert [ i ] = ( this - > Invert [ i + 1 ] * ( i + 1 ) ) % P ; if ( N = = K | | K = ...

1978 : 소수 찾기 [C++]

# include < iostream > # include < vector > using namespace std ; int main ( ) { cin . tie ( NULL ) ; vector < int > Primes ; Primes . push_back ( 2 ) ; Primes . push_back ( 3 ) ; for ( int i = 4 ; i < 1000 ; i + + ) { bool IsPrime = true ; if ( i % 2 = = 0 | | i % 3 = = 0 ) continue ; for ( int j = 4 ; j < i ; j + + ) { if ( i % j = = 0 ) { IsPrime = false ; break ; } } if ( IsPrime ) Primes . push_back ( i ) ; } int N , Count = 0 ; cin > > N ; for ( int i = 0 ; i < N ; i + + ) { int Input ; cin > > Input ; for ( int j = 0 ; j < Primes . size ( ) ; j + + ) if ( Input = = Primes [ j ] ) Count + + ; } cout < < Count < < " \n " ; return 0 ; }