POJ 2992 Divisors 求一個數的因數的個數

來源:互聯網
上載者:User

題意:輸入C(n,k), 求該數的因數個數。

題解:對於任意質數p, n!中有(n/p+n/p^2+n/p^3+...)個質因子p。

#include<iostream>#include<cstdio>#include<cstring>using namespace std;const int MAX = 450;int a[MAX], p[MAX], pn;void prime (){    pn = 0;    memset(a,0,sizeof(a));    int i, j;    for ( i = 2; i < MAX; i++ )    {        if ( a[i] == 0 )            p[pn++] = i;        for ( j = 0; j < pn && i * p[j] < MAX && (p[j] <= a[i] || a[i] == 0); j++ )            a[i*p[j]] = p[j];    }}int count ( int x, int pri ){    int ret = 0, tmp = pri;    while ( x >= tmp )    {        ret += x / tmp;        tmp *= pri;    }    return ret;}int main(){    int n, k;    prime();    while ( scanf("%d%d",&n,&k) != EOF )    {        int a, b, c;        __int64 res = 1;        for ( int i = 0; i < pn && p[i] <= n; i++ )        {            a = count ( n, p[i] );            b = count ( k, p[i] );            c = count ( n-k, p[i] );            res *= ( a - b - c + 1 );        }        printf("%I64d\n", res );    }    return 0;}

 

聯繫我們

該頁面正文內容均來源於網絡整理,並不代表阿里雲官方的觀點,該頁面所提到的產品和服務也與阿里云無關,如果該頁面內容對您造成了困擾,歡迎寫郵件給我們,收到郵件我們將在5個工作日內處理。

如果您發現本社區中有涉嫌抄襲的內容,歡迎發送郵件至: info-contact@alibabacloud.com 進行舉報並提供相關證據,工作人員會在 5 個工作天內聯絡您,一經查實,本站將立刻刪除涉嫌侵權內容。

A Free Trial That Lets You Build Big!

Start building with 50+ products and up to 12 months usage for Elastic Compute Service

  • Sales Support

    1 on 1 presale consultation

  • After-Sales Support

    24/7 Technical Support 6 Free Tickets per Quarter Faster Response

  • Alibaba Cloud offers highly flexible support services tailored to meet your exact needs.