結果

問題 No.36 素数が嫌い！
ユーザー tonegawatonegawa
提出日時 2020-08-16 06:27:18
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 3 ms / 5,000 ms
コード長 2,811 bytes
コンパイル時間 1,370 ms
コンパイル使用メモリ 124,084 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-04-19 18:15:20
合計ジャッジ時間 2,409 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 1 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 1 ms
5,376 KB
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 1 ms
5,376 KB
testcase_10 AC 1 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 3 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 1 ms
5,376 KB
testcase_17 AC 1 ms
5,376 KB
testcase_18 AC 1 ms
5,376 KB
testcase_19 AC 1 ms
5,376 KB
testcase_20 AC 1 ms
5,376 KB
testcase_21 AC 1 ms
5,376 KB
testcase_22 AC 2 ms
5,376 KB
testcase_23 AC 1 ms
5,376 KB
testcase_24 AC 1 ms
5,376 KB
testcase_25 AC 1 ms
5,376 KB
testcase_26 AC 1 ms
5,376 KB
testcase_27 AC 1 ms
5,376 KB
testcase_28 AC 2 ms
5,376 KB
testcase_29 AC 2 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string>
#include <vector>
#include <queue>
#include <deque>
#include <algorithm>
#include <set>
#include <map>
#include <bitset>
#include <cmath>
#include <functional>
#include <iomanip>
#define vll vector<ll>
#define vvvl vector<vvl>
#define vvl vector<vector<ll>>
#define VV(a, b, c, d) vector<vector<d>>(a, vector<d>(b, c))
#define VVV(a, b, c, d) vector<vvl>(a, vvl(b, vll (c, d)));
#define re(c, b) for(ll c=0;c<b;c++)
#define rep(a,b,c) for(ll a=b;a<c;a++)
#define all(obj) (obj).begin(), (obj).end()
typedef long long int ll;
typedef long double ld;
using namespace std;

void get(vll &a){re(i,a.size()) scanf("%lld",&a[i]);}
void get(vvl &a){re(i,a.size()) re(j,a[0].size()) scanf("%lld",&a[i][j]);}
void print(vll &a){re(i,a.size()) cout<<a[i]<<(i==a.size()-1?"\n":" ");}
void print(vvl &a){re(i,a.size()) re(j,a[0].size())cout<<a[i][j]<<(j==a[0].size()-1?"\n":" ");}

namespace fact{
  ll gcd(ll _a, ll _b) {
    unsigned long long a = abs(_a), b = abs(_b);
    if(a == 0) return b;
    if(b == 0) return a;
    int shift = __builtin_ctzll(a | b);
    a >>= __builtin_ctzll(a);
    do {
      b >>= __builtin_ctzll(b);
      if(a > b) std::swap(a, b);
      b -= a;
    }while(b);
    return (a << shift);
  }
  template <class T, class U>
  T modpow(T x, U n, T md) {
      T r = 1 % md;
      x %= md;
      while(n) {
          if(n & 1) r = (r * x) % md;
          x = (x * x) % md;
          n >>= 1;
      }
      return r;
  }
  bool is_prime(ll n) {
    if(n<=1) return false;
    if(n==2) return true;
    if(n%2==0) return false;
    ll d = n - 1;
    while(d % 2 == 0) d /= 2;
    for(ll a : {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37}) {
      if(n <= a) break;
      ll t = d, y = modpow<__int128_t>(a, t, n);
      while(t != n - 1 && y != 1 && y != n - 1) {
        y = __int128_t(y) * y % n, t <<= 1;
      }
      if(y != n - 1 && t % 2 == 0) {
        return false;
      }
    }
    return true;
  }
  ll rho(ll n){
    if(is_prime(n)) return n;
    if(n % 2 == 0) return 2;
    ll c = 0;
    while(true){
      c = (3 * c + 1) % 24;
      auto f = [&](ll x){return(__int128_t(x)*x+c)%n;};
      ll st = 0;
      while(true){
        st++;
        ll x = st%n, y = f(x);
        while(true){
          ll g = gcd(abs(x-y), n);
          if(g==1) x = f(x), y = f(f(y));
          else if(g==n) break;
          else return g;
        }
      }
    }
  }
  vll factorize(ll n) {
    if(n == 1) return {};
    ll x = rho(n);
    if(x == n) return {x};
    vll le = factorize(x);
    vll ri = factorize(n / x);
    le.insert(le.end(), ri.begin(), ri.end());
    return le;
  }
}

int main(int argc, char const *argv[]) {
  ll n;std::cin >> n;
  vll d = fact::factorize(n);
  if(d.size()>=3) std::cout << "YES" << '\n';
  else std::cout << "NO" << '\n';
}
0