結果
問題 | No.1424 Ultrapalindrome |
ユーザー |
|
提出日時 | 2021-03-12 22:15:39 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 10,551 bytes |
コンパイル時間 | 2,287 ms |
コンパイル使用メモリ | 209,604 KB |
最終ジャッジ日時 | 2025-01-19 14:58:00 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 WA * 2 |
ソースコード
#include <bits/stdc++.h>#define REP(i, n) for(int i = 0;i < n;i++)#define ll long longusing namespace std;typedef pair<int, int> P;typedef pair<ll,ll> LP;const int dx[8] = {1, 0, -1, 0, 1, -1, -1, 1};const int dy[8] = {0, 1, 0, -1, 1, 1, -1, -1};const int INF = 1000000000;const ll LINF = 1000000000000000000;//1e18const double PI = acos(-1.0);const double EPS = 1e-10;template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; }template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; }inline constexpr ll ceil_div(const ll a, const ll b) {return (a + b - 1) / b - !!((a + b - 1) % b < 0);}// return ceil(a/b)inline constexpr ll floor_div(const ll a, const ll b) {return a / b - !!(a % b < 0);}// return floor(a/b)/* ----------------------- DEBUG FUNCTION ---------------------------- */#define DUMPOUT cerrvoid dump_function() { DUMPOUT << ' '; }void dump_function(bool a) { DUMPOUT << a; }void dump_function(int a) { DUMPOUT << a; }void dump_function(long long a) { DUMPOUT << a; }void dump_function(char a) { DUMPOUT << a; }void dump_function(string &a) { DUMPOUT << a; }void dump_function(double a) { DUMPOUT << a; }template <class T> void dump_function(const vector<T> &);template <class T, size_t size> void dump_function(const array<T, size> &);template <class T, class L> void dump_function(const pair<T, L> &p);template <class T, size_t size> void dump_function(const T (&)[size]);template <class T> void dump_function(const vector<T> &a) {if(a.empty()) return;dump_function(a[0]);for(auto i = a.begin(); ++i != a.end();) {DUMPOUT << " ";dump_function(*i);}DUMPOUT << endl;}template <class T> void dump_function(const deque<T> &a) {if(a.empty()) return;dump_function(a[0]);for(auto i = a.begin(); ++i != a.end();) {DUMPOUT << " ";dump_function(*i);}}template <class T, size_t size> void dump_function(const array<T, size> &a) {dump_function(a[0]);for(auto i = a.begin(); ++i != a.end();) {DUMPOUT << " ";dump_function(*i);}}template <class T, class L> void dump_function(const pair<T, L> &p) {DUMPOUT << '(';dump_function(p.first);DUMPOUT << ",";dump_function(p.second);DUMPOUT << ')';}template <class T> void dump_function(set<T> &x) {for(auto e : x) dump_function(e), DUMPOUT << " ";DUMPOUT << endl;}template <class T> void dump_function(multiset<T> &x) {for(auto e : x) dump_function(e), DUMPOUT << " ";DUMPOUT << endl;}template <class T, size_t size> void dump_function(const T (&a)[size]) {dump_function(a[0]);for(auto i = a; ++i != end(a);) {DUMPOUT << " ";dump_function(*i);}}template <class T> void dump_function(const T &a) { DUMPOUT << a; }int dump_out() {DUMPOUT << '\n';return 0;}template <class T> int dump_out(const T &t) {dump_function(t);DUMPOUT << '\n';return 0;}template <class Head, class... Tail> int dump_out(const Head &head, const Tail &... tail) {dump_function(head);DUMPOUT << ' ';dump_out(tail...);return 0;}#ifdef DEBUG_#define dump(x)\DUMPOUT << #x << ": ";\dump_function(x);\DUMPOUT << endl;void dumps() {}template <class T> void dumps(const T &t) {dump_function(t);DUMPOUT << " ";}template <class Head, class... Tail> void dumps(const Head &head, const Tail &... tail) {dump_function(head);DUMPOUT << ' ';dump_out(tail...);}#else#define dump(x)template <class... T> void dumps(const T &...) {}#endif/* ----------------------- DEBUG FUNCTION ---------------------------- */template<int MOD> struct Fp {long long val;constexpr Fp(long long v = 0) noexcept : val(v % MOD) {if (val < 0) val += MOD;}constexpr int getmod() const { return MOD; }constexpr Fp operator - () const noexcept {return val ? MOD - val : 0;}constexpr Fp operator + (const Fp& r) const noexcept { return Fp(*this) += r; }constexpr Fp operator - (const Fp& r) const noexcept { return Fp(*this) -= r; }constexpr Fp operator * (const Fp& r) const noexcept { return Fp(*this) *= r; }constexpr Fp operator / (const Fp& r) const noexcept { return Fp(*this) /= r; }constexpr Fp& operator += (const Fp& r) noexcept {val += r.val;if (val >= MOD) val -= MOD;return *this;}constexpr Fp& operator -= (const Fp& r) noexcept {val -= r.val;if (val < 0) val += MOD;return *this;}constexpr Fp& operator *= (const Fp& r) noexcept {val = val * r.val % MOD;return *this;}constexpr Fp& operator /= (const Fp& r) noexcept {long long a = r.val, b = MOD, u = 1, v = 0;while (b) {long long t = a / b;a -= t * b, swap(a, b);u -= t * v, swap(u, v);}val = val * u % MOD;if (val < 0) val += MOD;return *this;}constexpr bool operator == (const Fp& r) const noexcept {return this->val == r.val;}constexpr bool operator != (const Fp& r) const noexcept {return this->val != r.val;}friend constexpr istream& operator >> (istream& is, Fp<MOD>& x) noexcept {is >> x.val;x.val %= MOD;if (x.val < 0) x.val += MOD;return is;}friend constexpr ostream& operator << (ostream& os, const Fp<MOD>& x) noexcept {return os << x.val;}friend constexpr Fp<MOD> modpow(const Fp<MOD>& r, long long n) noexcept {if (n == 0) return 1;auto t = modpow(r, n / 2);t = t * t;if (n & 1) t = t * r;return t;}friend constexpr Fp<MOD> modinv(const Fp<MOD>& r) noexcept {long long a = r.val, b = MOD, u = 1, v = 0;while (b) {long long t = a / b;a -= t * b, swap(a, b);u -= t * v, swap(u, v);}return Fp<MOD>(u);}};template<class T> struct BiCoef {vector<T> fact_, inv_, finv_;constexpr BiCoef() {}constexpr BiCoef(int n) noexcept : fact_(n, 1), inv_(n, 1), finv_(n, 1) {init(n);}constexpr void init(int n) noexcept {fact_.assign(n, 1), inv_.assign(n, 1), finv_.assign(n, 1);int MOD = fact_[0].getmod();for(int i = 2; i < n; i++){fact_[i] = fact_[i-1] * i;inv_[i] = -inv_[MOD%i] * (MOD/i);finv_[i] = finv_[i-1] * inv_[i];}}constexpr T com(int n, int k) const noexcept {if (n < k || n < 0 || k < 0) return 0;return fact_[n] * finv_[k] * finv_[n-k];}constexpr T fact(int n) const noexcept {if (n < 0) return 0;return fact_[n];}constexpr T inv(int n) const noexcept {if (n < 0) return 0;return inv_[n];}constexpr T finv(int n) const noexcept {if (n < 0) return 0;return finv_[n];}constexpr T perm(int n, int k) const noexcept {if (n < k || n < 0 || k < 0) return 0;return fact_[n] * finv_[n-k];}};/* ----------------------------- MOD ----------------------------------- */const int MOD = 1000000007;const int MOD2 = 998244353;using mint = Fp<MOD>;BiCoef<mint> bc;// using vec = vector<int>;// using mat = vector<vec>;/* ----------------------------- MOD ----------------------------------- *//* ----------------------- AtCoder Library ---------------------------- */// #include <atcoder/all>// using namespace atcoder;/* ----------------------- AtCoder Library ---------------------------- */struct treeDiameter{int N;vector<int> dist;vector<vector<int>> G;int v, w; // 直径のパスのv, wint diameter; // パスv, wの直径の長さtreeDiameter(const int& N) : N(N), dist(N, INF), G(N), v(0), diameter(0){}void add_edge(int u, int v){G[u].emplace_back(v);G[v].emplace_back(u);}pair<int, int> bfs(int root, vector<int>& dist){int v = root;int ma = 0;queue<tuple<int, int, int>> q;q.emplace(0, root, -1);while(!q.empty()){auto [d, cur, par] = q.front(); q.pop();dist[cur] = d;if(dist[cur] > ma){ma = dist[cur];v = cur;}for(int nv : G[cur]){if(nv == par) continue;q.emplace(d+1, nv, cur);}}return make_pair(ma, v);}void build(int root){dist.assign(N, INF);pair<int, int> p = bfs(root, dist);w = p.second;dist.assign(N, INF);pair<int, int> res = bfs(p.second, dist);diameter = res.first;v = res.second;}// 木の中心を求める(最遠距離を最小にする点、高々2個)vector<int> get_center(int root){vector<int> res;vector<int> dist_v(N, INF);vector<int> dist_w(N, INF);auto res_v = bfs(v, dist_v);auto res_w = bfs(w, dist_w);for(int i=0;i<N;i++){if(dist_v[i] == diameter / 2 and dist_w[i] == diameter - diameter / 2){res.emplace_back(i);}}return res;}};void solve(){int N;cin >> N;assert(N != 1);treeDiameter tree(N);vector<int> deg(N);REP(i,N-1){int a, b;cin >> a >> b;tree.add_edge(--a, --b);deg[a]++, deg[b]++;}tree.build(0);int v = tree.v;int w = tree.w;tree.build(v);auto d1 = tree.dist;tree.build(w);auto d2 = tree.dist;REP(i,N) if(i != v and i != w and deg[i] == 1){if(d1[i] != tree.diameter or d2[i] != tree.diameter){cout << "No" << endl;return;}}cout << "Yes" << endl;}int main(){cin.tie(0);ios::sync_with_stdio(false);solve();// int T; cin >> T; REP(t,T) solve();}