#include #include #include using namespace std; using namespace atcoder; using mint = modint; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf32 1000000005 #define Inf64 4000000000000000001LL int main(){ int N,P; cin>>N>>P; mint::set_mod(P); vector dp(N+1); mint sum = 0; for(int i=3;i<=N;i++){ mint ts = 0; int tc = 0; vector ys; for(int j=1;j*j<=i;j++){ if(i%j==0){ ys.push_back(j); if(j != i/j) ys.push_back(i/j); } } ts =sum; rep(j,ys.size()){ ts -= dp[ys[j]]; } mint tdp = ts *mint(i - ys.size()).inv(); mint tp = i-ys.size(); tp /= i; tdp += tp.inv(); rep(j,ys.size()){ sum -= dp[ys[j]]; dp[ys[j]] = tdp; sum += tdp; } } cout<