fn main() {
    proconio::input! {
        n: usize,
        p: usize,
    }

    let mut divisors = vec![vec![]; n + 1];
    for i in 1..=n {
        for j in (i..=n).step_by(i) {
            divisors[j].push(i);
        }
    }

    let mut invs = vec![0; n + 1];
    invs[1] = 1;
    for i in 2..=n {
        invs[i] = (p - p / i) * invs[p % i] % p;
    }

    let mut dp = vec![0; n + 1];
    let mut total = 0;
    for i in 3..=n {
        let d_sum = divisors[i].iter().map(|&d| dp[d]).sum::<usize>() % p;
        let dp_i = (total + i + p - d_sum) * invs[i - divisors[i].len()] % p;
        for &d in &divisors[i] {
            total = (total + dp_i + p - dp[d]) % p;
            dp[d] = dp_i;
        }
    }

    println!("{}", dp[n]);
}