結果
問題 | No.811 約数の個数の最大化 |
ユーザー |
![]() |
提出日時 | 2019-04-13 16:22:01 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 12 ms / 2,000 ms |
コード長 | 8,669 bytes |
コンパイル時間 | 1,697 ms |
コンパイル使用メモリ | 164,992 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-15 11:11:06 |
合計ジャッジ時間 | 2,273 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 12 |
ソースコード
#include <bits/stdc++.h>// #include <boost/multiprecision/cpp_int.hpp>#define int long long#define inf 1000000007#define pa pair<int,int>#define ll long long#define pal pair<double,double>#define ppap pair<pa,int>#define PI 3.14159265358979323846#define paa pair<int,char>#define mp make_pair#define pb push_back#define EPS (1e-10)int dx[8]={0,1,0,-1,1,1,-1,-1};int dy[8]={1,0,-1,0,-1,1,1,-1};using namespace std;class pa3{public:int x;int y,z;pa3(int x=0,int y=0,int z=0):x(x),y(y),z(z) {}bool operator < (const pa3 &p) const{if(x!=p.x) return x<p.x;if(y!=p.y) return y<p.y;return z<p.z;//return x != p.x ? x<p.x: y<p.y;}bool operator > (const pa3 &p) const{if(x!=p.x) return x>p.x;if(y!=p.y) return y>p.y;return z>p.z;//return x != p.x ? x<p.x: y<p.y;}bool operator == (const pa3 &p) const{return x==p.x && y==p.y && z==p.z;}bool operator != (const pa3 &p) const{return !( x==p.x && y==p.y && z==p.z);}};class pa4{public:int x;int y,z,w;pa4(int x=0,int y=0,int z=0,int w=0):x(x),y(y),z(z),w(w) {}bool operator < (const pa4 &p) const{if(x!=p.x) return x<p.x;if(y!=p.y) return y<p.y;if(z!=p.z)return z<p.z;return w<p.w;//return x != p.x ? x<p.x: y<p.y;}bool operator > (const pa4 &p) const{if(x!=p.x) return x>p.x;if(y!=p.y) return y>p.y;if(z!=p.z)return z>p.z;return w>p.w;//return x != p.x ? x<p.x: y<p.y;}bool operator == (const pa4 &p) const{return x==p.x && y==p.y && z==p.z &&w==p.w;}};class pa2{public:int x,y;pa2(int x=0,int y=0):x(x),y(y) {}pa2 operator + (pa2 p) {return pa2(x+p.x,y+p.y);}pa2 operator - (pa2 p) {return pa2(x-p.x,y-p.y);}bool operator < (const pa2 &p) const{return y != p.y ? y<p.y: x<p.x;}bool operator > (const pa2 &p) const{return x != p.x ? x<p.x: y<p.y;}bool operator == (const pa2 &p) const{return abs(x-p.x)==0 && abs(y-p.y)==0;}bool operator != (const pa2 &p) const{return !(abs(x-p.x)==0 && abs(y-p.y)==0);}};string itos( int i ) {ostringstream s ;s << i ;return s.str() ;}int gcd(int v,int b){if(v>b) return gcd(b,v);if(v==b) return b;if(b%v==0) return v;return gcd(v,b%v);}int mod;int extgcd(int a, int b, int &x, int &y) {if (b == 0) {x = 1;y = 0;return a;}int d = extgcd(b, a%b, y, x);y -= a/b * x;return d;}pa operator+(const pa & l,const pa & r) {return {l.first+r.first,l.second+r.second};}pa operator-(const pa & l,const pa & r) {return {l.first-r.first,l.second-r.second};}int pr[1200010];int inv[1200010];int beki(int wa,int rr,int warukazu){if(rr==0) return 1%warukazu;if(rr==1) return wa%warukazu;wa%=warukazu;if(rr%2==1) return ((ll)beki(wa,rr-1,warukazu)*(ll)wa)%warukazu;ll zx=beki(wa,rr/2,warukazu);return (zx*zx)%warukazu;}int comb(int nn,int rr){int r=pr[nn]*inv[rr];r%=mod;r*=inv[nn-rr];r%=mod;return r;}void gya(int ert){pr[0]=1;for(int i=1;i<=ert;i++){pr[i]=(pr[i-1]*i)%mod;}inv[ert]=beki(pr[ert],mod-2,mod);for(int i=ert-1;i>=0;i--){inv[i]=inv[i+1]*(i+1)%mod;}}// cin.tie(0);// ios::sync_with_stdio(false);//priority_queue<pa3,vector<pa3>,greater<pa3>> pq;//sort(ve.begin(),ve.end(),greater<int>());//----------------kokomade tenpure------------vector<int>sosu,kosu;int a[100020]={};int yaku[100020]={};signed main(){cin.tie(0);ios::sync_with_stdio(false);int n;cin>>n;int nn=n;for(int i=2;i*i<=n;i++)if(nn%i==0){int c=0;while(nn%i==0){nn/=i;c++;}sosu.pb(i);kosu.pb(c);}if(nn>1){sosu.pb(nn);kosu.pb(1);}int l=sosu.size();for(int i=0;i<l;i++){// cout<<sosu[i]<<" "<<kosu[i]<<endl;for(int j=1;j<=n;j++){int jj=j;int c=0;while(jj%sosu[i]==0){jj/=sosu[i];c++;}a[j]+=min(c,kosu[i]);}}for(int i=1;i<=n;i++){for(int j=i;j<=n;j+=i)yaku[j]++;}int val=-1;int k;cin>>k;int ans=-1;for(int i=1;i<n;i++)if(a[i]>=k){// cout<<i<<endl;if(val<yaku[i]){val=yaku[i];ans=i;}}// cout<<val<<endl;cout<<ans<<endl;return 0;}