結果

問題 No.1611 Minimum Multiple with Double Divisors
コンテスト
ユーザー 蜜蜂
提出日時 2021-07-21 21:27:21
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,853 bytes
コンパイル時間 2,210 ms
コンパイル使用メモリ 180,752 KB
実行使用メモリ 14,208 KB
最終ジャッジ日時 2024-07-17 16:01:17
合計ジャッジ時間 8,908 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 2
other TLE * 1 -- * 36
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'std::vector<long long int> divisor(ll)':
main.cpp:63:10: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
   63 |     auto [a,b]=fac[i];
      |          ^

ソースコード

diff # 

//g++ CF3.cpp -std=c++14 -O2 -I .
#include <bits/stdc++.h>
using namespace std;

using ll = long long;
using ld = long double;

using vi = vector<int>;
using vvi = vector<vi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using vld = vector<ld>;
using vvld = vector<vld>;

#define fi first
#define se second
#define pb push_back
#define pq_big(T) priority_queue<T,vector<T>,less<T>>
#define pq_small(T) priority_queue<T,vector<T>,greater<T>>
#define all(a) a.begin(),a.end()
#define rep(i,start,end) for(ll i=start;i<(ll)(end);i++)
#define per(i,start,end) for(ll i=start;i>=(ll)(end);i--)

//nを素因数分解 {素因数,個数} で管理 (1は空配列)
vector<pair<ll,ll>> factor(ll n){
  vector<pair<ll,ll>> res;

  int a=0;
  while(n%2==0){
    n/=2;
    a++;
  }

  if(a>0){
    res.pb({2,a});
  }

  for(ll x=3;x*x<=n;x+=2){
    if(n%x==0){
      a=0;
      while(n%x==0){
        n/=x;
        a++;
      }
      res.pb({x,a});
    }
  }

  if(n>1){
    res.pb({n,1});
  }

  return res;
}

//nの約数配列(昇順)
vector<ll> divisor(ll n){
  auto fac=move(factor(n));
  int sz=fac.size();
  
  vector<vector<ll>> x(sz);
  for(int i=0;i<sz;i++){
    auto [a,b]=fac[i];
    ll now=1;
    for(int j=0;j<=b;j++){
      x[i].pb(now);
      now*=a;
    }
  }
  
  vector<vector<ll>> res(sz);
  res[0]=x[0];
  
  for(int i=1;i<sz;i++){
    for(int j=0;j<res[i-1].size();j++){
      for(int k=0;k<x[i].size();k++){
        res[i].pb(res[i-1][j]*x[i][k]);
      }
    }
  }

  sort(all(res[sz-1]));
  
  return res[sz-1];
}

void solve(){
  ll n;
  cin>>n;

  auto x=divisor(n);
  int z=x.size();

  rep(i,2,100){
    auto check=divisor(i*n);
    if(check.size()==2*z){
      cout<<i*n<<"\n";
      return;
    }
  }
}

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  int t;
  cin>>t;
  //t=1;
  while(t--){
    solve();
  }
}
0