#include using namespace std; using ll =long long; #define all(v) v.begin(),v.end() #define rep(i,a,b) for(int i=a;i=b;i--) ll INF=2e18; struct UnionFind { vector> par; vector r; void init(ll n) { par=vector> (n); r=vector (n); for(ll i=0;i= mod) x -= mod; return *this; } mint& operator-=(const mint& a) { if ((x += mod-a.x) >= mod) x -= mod; return *this; } mint& operator*=(const mint& a) { (x *= a.x) %= mod; return *this; } mint operator+(const mint& a) const { mint res(*this); return res+=a; } mint operator-(const mint& a) const { mint res(*this); return res-=a; } mint operator*(const mint& a) const { mint res(*this); return res*=a; } mint pow(ll t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime mod mint inv() const { return pow(mod-2); } mint& operator/=(const mint& a) { return (*this) *= a.inv(); } mint operator/(const mint& a) const { mint res(*this); return res/=a; } friend ostream& operator<<(ostream& os, const mint& m){ os << m.x; return os; } }; int main() { ios::sync_with_stdio(false); cin.tie(0); ll N,M;cin>>N>>M; vector P(N+1); for(ll i=1;i<=N;i++) P[i]=i; for(ll i=0;i>T; vector S(T); for(ll j=0;j>S[j]; ll memo=P[S[T-1]]; for(ll j=T-1;j>=1;j--) { P[S[j]]=P[S[j-1]]; } P[S[0]]=memo; } UnionFind uf; uf.init(N+1); for(ll i=1;i<=N;i++) { uf.unite(P[i],i); } vector note(N+1); vector vec(0); for(ll i=1;i<=N;i++) { ll k=uf.find(i); if(!note[k]) { note[k]=true; vec.push_back(uf.par[k].second); } } vector prime(0); vector memo(1000000+100); for(ll i=2;i<=1000000;i++) { if(!memo[i]) { prime.push_back(i); for(ll j=1;j*i<=1000000;j++) memo[j*i]=true; } } vector>> d(1000000+100,vector> (0)); for(auto x:prime) { for(ll j=1;j*x<=1000000;j++) { ll count=0; ll now=j*x; while(now%x==0) { now/=x; count++; } d[j*x].push_back(make_pair(x,count)); } } mint ans=1; vector ma(1000000+10); for(ll i=0;i