#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #if __has_include() #include #else #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #endif #define rep(i,n) for(int i=0;i=0;i--) #define repv(i,v) for(auto i:v) #define all(x) x.begin(),x.end() templateinline bool chmax(T &a,T b){return((ainline bool chmin(T &a,T b){return((a>b)?(a=b,true):(false));} using namespace std; using ll=long long; using vi=vector;using vvi=vector;using vvvi=vector; using vl=vector;using vvl=vector;using vvvl=vector; using pii=pair;using vpii=vector;using pll=pair;using vpll=vector;using pli=pair;using pil=pair; using mpi=map;using mps=map;using mpc=map;using mpl=map; using sti=set;using stc=set;using sts=set;using stl=set; #if __has_include() #include using namespace atcoder; using mint=modint998244353; using mint1=modint1000000007; #endif using bll = __int128_t; // O(√n) verified bool is_prime(ll n){if(n<2)return false;if(n==2)return true;if(n%2==0)return false;for(ll i=3;i*i<=n;i+=2)if(n%i==0)return false;return true;} // O(√n) verified vectordivisors(ll n){vectorr;int i;for(i=1;i*i>prime_factorize(ll n){vector>r;int c=0;while(n%2==0){n/=2;c++;}r.push_back({2,c});for(ll p=3;p*p<=n;p+=2){if(n%p!=0)continue;int e=0;while(n%p==0){e++;n/=p;}r.push_back({p,e});}if(n!=1)r.push_back({n,1});return r;} // エラトステネスの篩 // O(n*log(log(n))) verified vectorEraSieve(ll n){vectorp(n+1,true);p[0]=p[1]=false;for(ll i=2;i*i<=n;i++){if(!p[i])continue;for(int q=i*2;q<=n;q+=i){p[q]=false;}}return p;} // エラトステネスの区間篩 verified // a以上b未満のcに対してp[c - a] vectorEraSecSieve(ll a,ll b){ll sq=(ll)ceil(sqrt(b));vectorps(sq,true);ps[0]=ps[1]=false;for(ll i=2;ip(b-a,true);for(ll i=2;iT modpow(T a,T n,T m){T r=1;while(n){if(n&1)r=r*a%m;a=a*a%m;n>>=1;}return r;} // 失敗率 1/(4^30) bool MillerRabin(ll n) {if(n==2)return true;if(n==1||(n>2&&n%2==0))return false;ll s=0,t=n-1;while(t%2==0){s++;t>>=1;}random_device gen;mt19937_64 mtrand(gen());for(int i=0;i<30;i++){ll a=(mtrand())%(n-1)+1;if(a%n==0)continue;bool ok=false;if(modpow(a,t,n)==1)ok=true;else for(int i=0;i(a,modpow(2,i,1ull<<63)*t,n)==n-1){ok=true;break;}if(!ok)return false;}return true;} ll _findprime(ll n){if(n%2==0)return 2;for(int c=1;cll {return((bll)a*a+c)%n;};ll x=0,y=0,g=1;while(g==1){x=f(x);y=f(f(y));g=gcd(abs(x-y),n);}if(g==n)continue;if(MillerRabin(g))return g;else if(MillerRabin(n/g))return n/g;else return _findprime(g);}} vector>PollardRho(ll n){vector>r;while(!MillerRabin(n)&&n>1){ll p=_findprime(n),s=0;while(n%p==0){n/=p;s++;}r.push_back({p,s});}if(n>1)r.push_back({n,1});return r;} // memo:PollardRhoからvectorを作った後にsortするのを忘れない int main() { ll n; cin >> n; if (n == 1) cout << "NO\n"; else cout << (MillerRabin(n) ? "NO\n" : "YES\n"); }