結果
問題 | No.811 約数の個数の最大化 |
ユーザー | ahe100 |
提出日時 | 2019-04-13 00:25:44 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 70 ms / 2,000 ms |
コード長 | 1,490 bytes |
コンパイル時間 | 1,389 ms |
コンパイル使用メモリ | 167,688 KB |
実行使用メモリ | 14,012 KB |
最終ジャッジ日時 | 2023-10-13 09:22:18 |
合計ジャッジ時間 | 2,623 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 4 ms
7,988 KB |
testcase_01 | AC | 5 ms
8,040 KB |
testcase_02 | AC | 61 ms
13,504 KB |
testcase_03 | AC | 5 ms
7,996 KB |
testcase_04 | AC | 5 ms
8,076 KB |
testcase_05 | AC | 6 ms
8,200 KB |
testcase_06 | AC | 10 ms
8,608 KB |
testcase_07 | AC | 12 ms
8,784 KB |
testcase_08 | AC | 33 ms
10,920 KB |
testcase_09 | AC | 36 ms
11,072 KB |
testcase_10 | AC | 25 ms
10,156 KB |
testcase_11 | AC | 63 ms
13,396 KB |
testcase_12 | AC | 21 ms
9,704 KB |
testcase_13 | AC | 69 ms
14,012 KB |
testcase_14 | AC | 70 ms
13,968 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i,n) for(int i = 0;i<((int)(n));i++) #define reg(i,a,b) for(int i = ((int)(a));i<=((int)(b));i++) #define irep(i,n) for(int i = ((int)(n)-1);i>=0;i--) #define ireg(i,a,b) for(int i = ((int)(b));i>=((int)(a));i--) typedef long long ll; /* */ ll n,k,divisor_num[100010],prime_num[100010],ans=0,ans_num=0; vector<ll> fact[100010]; void factoring(ll n,vector<ll> v[]){ ll visit[300010]={}; reg(i,2,n){ if(visit[i]==0){ reg(j,1,n){ if(j*i>n)break; visit[i*j]=1; v[i*j].push_back(i);//素因数リスト(重複は含まず) } } } } void calc_divisor_num(ll n, vector<ll> fact[], ll divisor_num[]){ reg(i,2,n){ ll m=i; divisor_num[i]=1; rep(j,fact[i].size()){ ll t=0; while(m%fact[i][j]==0){ m/=fact[i][j]; t++; } divisor_num[i]*=(t+1); } } } void calc_prime_num(ll n, vector<ll> fact[], ll prime_num[]){ reg(i,1,n){ ll m=i; prime_num[i]=0; rep(j,fact[i].size()){ while(m%fact[i][j]==0){ m/=fact[i][j]; prime_num[i]++; } } } } //ユークリッドの互除法 long long gcd(long long a,long long b){ if(b==0)return a; return gcd(b,a%b); } void init(){ cin>>n>>k; factoring(n,fact); calc_divisor_num(n,fact,divisor_num); calc_prime_num(n,fact,prime_num); } int main(void){ init(); reg(i,1,n-1){ ll g = gcd(i,n); if(prime_num[g]>=k && ans_num<divisor_num[i]){ ans_num = divisor_num[i]; ans = i; } } cout<<ans<<endl; return 0; }