#include 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; /* */ struct D{ ll val,num,id; }; bool operator<(const D &l, const D &r){//全順序を定めないと,{1,0,0}=={1,1,0} if(l.num > r.num){ return true; }else if(l.num == r.num){ if(l.val > r.val){ return true; }else if(l.val == r.val){ return l.id > r.id; } } return false; } ll n,k,m; //2からprob_sizeまでの全ての数の素因数分解(ただし重複無視)を返す const ll prob_size=300000;//ここの100000は発生する数の範囲 vector fact[300010],pri; ll visited[prob_size]={}, a[100010]={},nu[100010],ans=1; //nまでの素数表を返す void prime(ll n, vector& p){ ll visited[300010]={}; reg(i,2,n){ if(visited[i]==0){ p.push_back(i); reg(j,2,n){ if(j*i>n)break; visited[i*j]=1; } } } } void init(){ cin>>n>>k; reg(i,2,prob_size){ if(visited[i]==0){ reg(j,1,prob_size){ if(j*i>prob_size)break; visited[i*j]=1; fact[i*j].push_back(i);//素因数リスト(重複は含まず) } } } prime(n,pri); rep(i,pri.size())nu[pri[i]]=i; m=n; cerr<<"初期化終了"<0){ rep(i,fact[n].size()){ if(m%fact[n][i]==0){ m/=fact[n][i]; ans*=fact[n][i]; a[nu[fact[n][i]]]++; // cerr<Q; rep(i,pri.size()){ Q.push({pri[i],a[i],i}); // cerr<