#if __INCLUDE_LEVEL__ == 0 #include __BASE_FILE__ void Solve() { int n; IN(n); vector p(n); IN(OneBased(p)); basic_string done(n, false); vector cyc; cyc.reserve(n); vector ans(n); for (int d : Rep1(1, n)) { for (int i = d; i <= n; i += d) { CsrArray::Add(i, d); } } auto ds = CsrArray::Build(n + 1); for (int s : Rep(0, n)) { if (done[s]) { continue; } cyc.clear(); { int i = s; do { cyc.push_back(i); i = p[i]; } while (i != s); } int g = 0; ranges::sort(cyc); for (int ci : Rep(1, Len(cyc))) { g = gcd(g, cyc[ci] - cyc[0]); } for (int d : ds[g]) { ans[d] += Len(cyc) - 1; } for (int i : cyc) { done[i] = true; } } ranges::for_each(ans | views::drop(1), LIFT(OUT)); } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); Solve(); } #elif __INCLUDE_LEVEL__ == 1 #include template class CsrArray { public: static void Reserve(int m) { buf_.reserve(m); } static void Add(int i, T x) { buf_.emplace_back(i, std::move(x)); } static CsrArray Build(int n) { CsrArray ret; ret.pos_.resize(n + 1); for (int i : buf_ | std::views::keys) { ++ret.pos_[i]; } std::partial_sum(ret.pos_.begin(), ret.pos_.end(), ret.pos_.begin()); ret.data_.resize(ret.pos_[n]); for (auto& [i, x] : buf_ | std::views::reverse) { ret.data_[--ret.pos_[i]] = std::move(x); } buf_.clear(); return ret; } int size() const { return int(pos_.size()) - 1; } auto operator[](int i) { return std::span(data_.data() + pos_[i], data_.data() + pos_[i + 1]); } auto operator[](int i) const { return std::span(data_.data() + pos_[i], data_.data() + pos_[i + 1]); } private: static thread_local inline std::vector> buf_; std::vector data_; std::vector pos_; }; template concept MyRange = std::ranges::range && !std::convertible_to && !std::convertible_to; template concept MyTuple = std::__is_tuple_like::value && !MyRange; namespace std { istream& operator>>(istream& is, MyRange auto&& r) { for (auto&& e : r) { is >> e; } return is; } istream& operator>>(istream& is, MyTuple auto&& t) { apply([&](auto&... xs) { (is >> ... >> xs); }, t); return is; } ostream& operator<<(ostream& os, MyRange auto&& r) { auto sep = ""; for (auto&& e : r) { os << exchange(sep, " ") << forward(e); } return os; } ostream& operator<<(ostream& os, MyTuple auto&& t) { auto sep = ""; apply([&](auto&... xs) { ((os << exchange(sep, " ") << xs), ...); }, t); return os; } } // namespace std template class OneBased { public: explicit OneBased(T&& x) : ref_(std::forward(x)) {} template requires(sizeof...(Ts) > 1) OneBased(Ts&&... xs) : ref_(std::forward_as_tuple(std::forward(xs)...)) {} friend std::istream& operator>>(std::istream& is, OneBased x) { if constexpr (MyRange) { for (auto&& e : x.ref_) { is >> ::OneBased(e); } } else if constexpr (MyTuple) { std::apply([&](auto&... xs) { (is >> ... >> ::OneBased(xs)); }, x.ref_); } else { is >> x.ref_; --x.ref_; } return is; } friend std::ostream& operator<<(std::ostream& os, OneBased x) { if constexpr (MyRange) { auto f = [](auto&& e) { return ::OneBased(std::forward(e)); }; os << (x.ref_ | std::views::transform(f)); } else if constexpr (MyTuple) { std::apply([&](auto&... xs) { os << std::tuple(::OneBased(xs)...); }, x.ref_); } else { os << ++x.ref_; --x.ref_; } return os; } private: T ref_; }; template OneBased(T&&) -> OneBased; template OneBased(Ts&&...) -> OneBased>; using namespace std; #define LIFT(f) ([&](auto&&... xs) -> decltype(auto) { return f(forward(xs)...); }) #define Rep(...) [](int l, int r) { return views::iota(min(l, r), r); }(__VA_ARGS__) #define Rep1(...) [](int l, int r) { return Rep(l, r + 1); }(__VA_ARGS__) #define Len(r) int(size(r)) #define IN(...) (cin >> forward_as_tuple(__VA_ARGS__)) #define OUT(...) (cout << forward_as_tuple(__VA_ARGS__) << '\n') #endif // __INCLUDE_LEVEL__ == 1