結果

問題 No.36 素数が嫌い！
ユーザー tonegawa
提出日時 2020-08-16 06:27:18
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 4 ms / 5,000 ms
コード長 2,811 bytes
コンパイル時間 1,478 ms
コンパイル使用メモリ 118,664 KB
最終ジャッジ日時 2025-01-13 01:54:27
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘void get(std::vector<long long int>&)’:
main.cpp:25:38: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   25 | void get(vll &a){re(i,a.size()) scanf("%lld",&a[i]);}
      |                                 ~~~~~^~~~~~~~~~~~~~
main.cpp: In function ‘void get(std::vector<std::vector<long long int> >&)’:
main.cpp:26:56: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   26 | void get(vvl &a){re(i,a.size()) re(j,a[0].size()) scanf("%lld",&a[i][j]);}
      |                                                   ~~~~~^~~~~~~~~~~~~~~~~

ソースコード

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