結果
問題 | No.811 約数の個数の最大化 |
ユーザー |
![]() |
提出日時 | 2019-04-13 00:25:44 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 76 ms / 2,000 ms |
コード長 | 1,490 bytes |
コンパイル時間 | 1,747 ms |
コンパイル使用メモリ | 169,612 KB |
実行使用メモリ | 14,080 KB |
最終ジャッジ日時 | 2024-09-15 06:29:34 |
合計ジャッジ時間 | 2,762 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 12 |
ソースコード
#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;}