#include #include #include #include #include #include #include "math.h" #include #include #include #define ifor(i,a,b) for (int i=(a);i<(b);i++) #define rfor(i,a,b) for (int i=(b)-1;i>=(a);i--) #define rep(i,n) for (int i=0;i<(n);i++) #define rrep(i,n) for (int i=(n)-1;i>=0;i--) using namespace std; typedef long double ld; typedef long long int lli; typedef complex P; const double eps = 1e-11; int vex[4]={1,0,-1,0}; int vey[4]={0,1,0,-1}; typedef vector Vec; typedef vector vec; typedef vector MAT; typedef vector mat; lli MOD=1000000007; //Ax=bをとくAは正方行列 //rankA<=min(m,n)ならば配列0のvecが帰る Vec gauss_jordan(const MAT& A,const Vec&b){ int n = A.size(); MAT B (n,Vec(n+1)); rep(i,n)rep(j,n)B[i][j]=A[i][j]; rep(i,n)B[i][n]=b[i]; rep(i,n){ int pivot = i; for(int j= i;jabs(B[pivot][i]))pivot = j; } swap(B[i],B[pivot]); if(abs(B[i][i])abs(B[pivot][i]))pivot = j; } if(i!=pivot)ans *= -1; swap(B[i],B[pivot]); if(abs(B[i][i])abs(B[pivot][i]))pivot = j; } swap(B[i],B[pivot]); if(abs(B[i][i]) euler(lli m){ vector p; lli M = m; for(int i = 2;i<=m;i++){ if(m%i==0){ p.push_back(i); while(m%i==0)m/=i; } if(M0){ if(p&1)ans=(ans*a) % mod; a=(a*a)%mod; p >>=1; } return ans%mod; } lli inv(lli a){ return powm(a,MOD-2,MOD); } lli gcd(lli A,lli B){ if(A%B==0)return B; else return gcd(B,A%B); } mat mul_mat_mod(mat A,mat B,lli m){ int n = A.size(); mat C (n,vec(n)); rep(i,n)rep(j,n)rep(k,n){ C[i][j]+=A[i][k]*B[k][j] %m; C[i][j]%=m; } return C; } mat pow_mat(mat A,lli p,lli mod){ int n= A.size(); mat B = mat(n,vec(n)); while(p>0){ if(p&1){ B=mul_mat_mod(A,B,mod); } A = mul_mat_mod(A,A,mod); p>>=1; } return B; } lli comb(lli a,lli b ){ lli ans =1; rep(i,b){ ans = ans %MOD* (a-i)%MOD*inv(b-i)%MOD; } return ans ; } int main(){ int N; cin >>N; vector p = euler(N); vector d; rep(i,p.size()){ int b =0; while(N%p[i]==0){ b++; N/=p[i]; } d.push_back(b); } int s = d[0]; for(int i=1;i