/* #include namespace my{entry void main(){ LL(N); LinearSievels(N); pp(ls.mobius_enumerate().sum()); }}*/ #include #if __has_include() #endif using namespace std; #define eb emplace_back #define RD(T,...) T __VA_ARGS__;lin(__VA_ARGS__) #define LL(...) RD(ll,__VA_ARGS__) #define fo(i,...) for(auto[i,i##stop,i##step]=for_range(0,__VA_ARGS__);ivoid pp(const auto&...a){[[maybe_unused]]const char*c="";((o<(a...);} #define entry defpp void main();void main2();}int main(){my::io();my::main();}namespace my{ namespace my{ void io(){cin.tie(nullptr)->sync_with_stdio(0);cout<constexpr auto for_range(T s,T a,T b,T c=1){return array{a-s,b,(1-s*2)*c};} void lin(auto&...a){(cin>>...>>a);} templatestruct pair{ A a;B b; pair()=default; pair(A a,B b):a(a),b(b){} auto operator<=>(const pair&)const=default; }; templateusing pack_back_t=tuple_element_t>; } namespace my{ templateconcept vectorial=is_base_of_v::value_type>,remove_cvref_t>; templatestruct core_t_helper{using type=T;}; templateusing core_t=core_t_helper::type; templatestruct vec; templatestruct hvec_helper{using type=vec::type>;}; templatestruct hvec_helper<0,T>{using type=T;}; templateusing hvec=hvec_helper::type; templatestruct vec:vector{ using C=core_t; using vector::vector; auto&emplace_back(auto&&...a){vector::emplace_back(std::forward(a)...);return*this;} templateauto fold(auto&&f)const{ pairres{}; fe(*this,e){ if constexpr(!vectorial){ if(res.b)f(res.a,e); else res={e,1}; }else { } } return res; } templateauto sum()const{return fold([](auto&a,const auto&b){a+=b;}).a;} }; templaterequires(sizeof...(A)>=2)vec(const A&...a)->vec>>; } namespace my{ ll mobius_prime_pow(ll,i8 k,ll){return-(k==1);} } namespace my{ templateclass LinearSieve{ public: T n; veclpf; veclpf_ord; veclpf_pow; veclpf_pow_except; vecprimes; LinearSieve(T n):n(n),lpf(n+1,-1),lpf_ord(n+1),lpf_pow(n+1),lpf_pow_except(n+1){ lpf[1]=lpf_ord[1]=lpf_pow[1]=lpf_pow_except[1]=1; fo(i,2,n+1){ if(lpf[i]==-1)primes.eb(lpf[i]=i); fe(primes,p){ if(p*i>n||p>lpf[i])break; lpf[p*i]=p; } int j=i/lpf[i]; lpf_ord[i]=lpf_ord[j]*(lpf[i]==lpf[j])+1; lpf_pow[i]=((lpf_pow[j]-1)*(lpf[i]==lpf[j])+1)*lpf[i]; lpf_pow_except[i]=i/lpf_pow[i]; } } auto multiplicative_function_enumerate(const auto&f)const{ vecr(n+1); r[1]=1; fo(i,2,n+1)r[i]=f(lpf[i],lpf_ord[i],lpf_pow[i])*r[lpf_pow_except[i]]; return r; } auto mobius_enumerate()const{return multiplicative_function_enumerate(mobius_prime_pow);} }; } namespace my{entry void main(){ LL(N); LinearSievels(N); pp(ls.mobius_enumerate().sum()); }}