use proconio::input; fn main() { input! { n: usize, } let primes = primes((n as f64).sqrt() as usize + 3); let mut ans = 0; for p in primes { let mut temp = p*p; while temp <= n { ans += temp; temp *= p; } } println!("{}", ans); } fn primes(max: usize) -> Vec { let mut is_prime = vec![true; max]; is_prime[0..=1].fill(false); for i in 2.. { if i*i >= max { break; } if is_prime[i] { let mut temp = i*2; while temp < max { is_prime[temp] = false; temp += i; } } } is_prime.into_iter().enumerate().filter(|&(_, is_prime)| is_prime).map(|(n, _)| n).collect() }