fn get_primes(n: usize) -> Vec { let mut flgs = vec![true; n+1]; flgs[0] = false; flgs[1] = false; for i in 2..=(n as f64).sqrt().ceil() as usize { if !flgs[i] { continue; } for j in i..=n/i { flgs[i*j] = false; } } flgs } fn main() { let mut n = String::new(); std::io::stdin().read_line(&mut n).ok(); let n: usize = n.trim().parse().unwrap(); let primes = get_primes((n as f64).sqrt().ceil() as usize); let mut factors = (0..=n).collect::>(); for i in 2..primes.len() { if !primes[i] { continue; } let size = i*i; for j in (size..=n).step_by(size) { while factors[j] % size == 0 { factors[j] /= size; } } } let mut result = 0usize; for i in 1..=n { let target = n / factors[i]; result += (target as f64).sqrt().floor() as usize; } println!("{}", result); }