結果

問題 No.847 Divisors of Power
コンテスト
ユーザー tarattata1
提出日時 2019-07-05 23:45:44
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,983 bytes
コンパイル時間 833 ms
コンパイル使用メモリ 77,772 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-06 23:21:56
合計ジャッジ時間 1,416 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main(int, char**)’:
main.cpp:102:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
  102 |     scanf("%d%d%d", &n,&k,&m);
      |     ~~~~~^~~~~~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <stdio.h>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <list>
#include <iterator>
#include <assert.h>
#pragma warning(disable:4996) 

typedef long long ll;
#define MIN(a, b) ((a)>(b)? (b): (a))
#define MAX(a, b) ((a)<(b)? (b): (a))
#define LINF 9223300000000000000
#define INF 2140000000
const long long MOD = 1000000007;
using namespace std;

ll gcd(ll a, ll b) {
    if(b == 0) return a;
    return gcd(b,a%b);
}
 
ll lcm(ll a, ll b) {
    ll g = gcd(a,b);
    return a / g * b;   // Be careful not to overflow
}


vector<int> pp;
void preparePrime( int maxp )
{
    int i;
    vector<int> a(maxp+1);
    for(i=2; i<=maxp; i++) {
        int k=i+i;
        while(k<=maxp) {
            a[k]=1;
            k+=i;
        }
    }
    for(i=2; i<=maxp; i++) {
        if(a[i]==0) pp.push_back( i );
    }
    return;
}

void decompositPrime( ll a, vector<pair<ll,ll> >& vv )
{
    preparePrime( (int)sqrt((double)a)+1 );

    int i;
    for(i=0; i<(int)pp.size(); i++) {
        int count = 0;
        while(a % pp[i] == 0 ) {
            a /= pp[i];
            count ++;
        }
        if(count>0) vv.push_back( make_pair( pp[i], count ) );
        if (a==1) break;
    }
    if(a>1) {
        vv.push_back( make_pair( a, 1 ) );
    }

    return;
}

int n,k,m;
vector<pair<ll,ll> > vv;
int ans;

void dfs( int i, ll curr )
{
    if(i==vv.size()) {
        if(curr<=m) {
#ifdef _DEBUG
            printf("%d ", curr);
#endif
            ans++;
        }
        return;
    }
    int p;
    ll tmp=curr;
    for(p=0; p<=vv[i].second*k; p++) {
        dfs( i+1, tmp );
        tmp*=vv[i].first;
        if(tmp>m) break;
    }

    return;
}

int main(int argc, char* argv[])
{
    scanf("%d%d%d", &n,&k,&m);
    decompositPrime( n, vv);

    ans=0;
    dfs( 0, 1 );
#ifdef _DEBUG
    printf("\n");
#endif
    printf("%d\n", ans);
    return 0;
}
0