結果

問題 No.1207 グラフX
ユーザー Ricky_ponRicky_pon
提出日時 2020-08-30 13:49:16
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 208 ms / 2,000 ms
コード長 4,518 bytes
コンパイル時間 2,535 ms
コンパイル使用メモリ 210,260 KB
実行使用メモリ 30,544 KB
最終ジャッジ日時 2023-08-09 12:02:15
合計ジャッジ時間 11,495 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 155 ms
16,488 KB
testcase_01 AC 157 ms
16,468 KB
testcase_02 AC 154 ms
16,544 KB
testcase_03 AC 153 ms
16,560 KB
testcase_04 AC 156 ms
16,436 KB
testcase_05 AC 208 ms
30,392 KB
testcase_06 AC 204 ms
30,428 KB
testcase_07 AC 201 ms
30,364 KB
testcase_08 AC 123 ms
12,784 KB
testcase_09 AC 132 ms
14,836 KB
testcase_10 AC 198 ms
23,384 KB
testcase_11 AC 200 ms
30,544 KB
testcase_12 AC 125 ms
13,196 KB
testcase_13 AC 66 ms
9,676 KB
testcase_14 AC 143 ms
16,124 KB
testcase_15 AC 130 ms
14,056 KB
testcase_16 AC 68 ms
9,804 KB
testcase_17 AC 99 ms
12,396 KB
testcase_18 AC 77 ms
12,608 KB
testcase_19 AC 98 ms
10,908 KB
testcase_20 AC 153 ms
16,160 KB
testcase_21 AC 14 ms
8,408 KB
testcase_22 AC 101 ms
12,832 KB
testcase_23 AC 108 ms
13,288 KB
testcase_24 AC 72 ms
11,996 KB
testcase_25 AC 151 ms
16,256 KB
testcase_26 AC 117 ms
13,668 KB
testcase_27 AC 140 ms
15,072 KB
testcase_28 AC 134 ms
14,760 KB
testcase_29 AC 127 ms
15,108 KB
testcase_30 AC 66 ms
10,564 KB
testcase_31 AC 57 ms
9,048 KB
testcase_32 AC 58 ms
10,764 KB
testcase_33 AC 63 ms
10,564 KB
testcase_34 AC 128 ms
13,880 KB
testcase_35 AC 15 ms
8,604 KB
testcase_36 AC 127 ms
14,536 KB
testcase_37 AC 111 ms
13,028 KB
testcase_38 AC 30 ms
9,392 KB
testcase_39 AC 63 ms
11,480 KB
testcase_40 AC 45 ms
8,596 KB
testcase_41 AC 91 ms
11,000 KB
testcase_42 AC 5 ms
8,264 KB
testcase_43 AC 4 ms
8,236 KB
testcase_44 AC 4 ms
8,240 KB
testcase_45 AC 143 ms
16,324 KB
testcase_46 AC 141 ms
16,512 KB
testcase_47 AC 143 ms
16,408 KB
testcase_48 AC 137 ms
16,484 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define For(i, a, b) for (int(i) = (int)(a); (i) < (int)(b); ++(i))
#define rFor(i, a, b) for (int(i) = (int)(a)-1; (i) >= (int)(b); --(i))
#define rep(i, n) For((i), 0, (n))
#define rrep(i, n) rFor((i), (n), 0)
#define fi first
#define se second
using namespace std;
typedef long long lint;
typedef unsigned long long ulint;
typedef pair<int, int> pii;
typedef pair<lint, lint> pll;
template <class T>
bool chmax(T &a, const T &b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}
template <class T>
bool chmin(T &a, const T &b) {
    if (a > b) {
        a = b;
        return true;
    }
    return false;
}
template <class T>
T div_floor(T a, T b) {
    if (b < 0) a *= -1, b *= -1;
    return a >= 0 ? a / b : (a + 1) / b - 1;
}
template <class T>
T div_ceil(T a, T b) {
    if (b < 0) a *= -1, b *= -1;
    return a > 0 ? (a - 1) / b + 1 : a / b;
}

constexpr lint mod = 1000000007;
constexpr lint INF = mod * mod;
constexpr int MAX = 200010;

template <int_fast64_t MOD>
struct modint {
    using i64 = int_fast64_t;
    i64 a;
    modint(const i64 a_ = 0) : a(a_) {
        if (a > MOD)
            a %= MOD;
        else if (a < 0)
            (a %= MOD) += MOD;
    }
    modint inv() {
        i64 t = 1, n = MOD - 2, x = a;
        while (n) {
            if (n & 1) (t *= x) %= MOD;
            (x *= x) %= MOD;
            n >>= 1;
        }
        modint ret(t);
        return ret;
    }
    bool operator==(const modint x) const { return a == x.a; }
    bool operator!=(const modint x) const { return a != x.a; }
    modint operator+(const modint x) const { return modint(*this) += x; }
    modint operator-(const modint x) const { return modint(*this) -= x; }
    modint operator*(const modint x) const { return modint(*this) *= x; }
    modint operator/(const modint x) const { return modint(*this) /= x; }
    modint operator^(const lint x) const { return modint(*this) ^= x; }
    modint &operator+=(const modint &x) {
        a += x.a;
        if (a >= MOD) a -= MOD;
        return *this;
    }
    modint &operator-=(const modint &x) {
        a -= x.a;
        if (a < 0) a += MOD;
        return *this;
    }
    modint &operator*=(const modint &x) {
        (a *= x.a) %= MOD;
        return *this;
    }
    modint &operator/=(modint x) {
        (a *= x.inv().a) %= MOD;
        return *this;
    }
    modint &operator^=(lint n) {
        i64 ret = 1;
        while (n) {
            if (n & 1) (ret *= a) %= MOD;
            (a *= a) %= MOD;
            n >>= 1;
        }
        a = ret;
        return *this;
    }
    modint operator-() const { return modint(0) - *this; }
    modint &operator++() { return *this += 1; }
    modint &operator--() { return *this -= 1; }
    bool operator<(const modint x) const { return a < x.a; }
};

using mint = modint<1000000007>;

vector<mint> fact;
vector<mint> revfact;

void setfact(int n) {
    fact.resize(n + 1);
    revfact.resize(n + 1);
    fact[0] = 1;
    rep(i, n) fact[i + 1] = fact[i] * mint(i + 1);

    revfact[n] = fact[n].inv();
    for (int i = n - 1; i >= 0; i--) revfact[i] = revfact[i + 1] * mint(i + 1);
}

mint getC(int n, int r) {
    if (n < r) return 0;
    return fact[n] * revfact[r] * revfact[n - r];
}

typedef struct UnionFindTree {
    vector<int> par;

    UnionFindTree(int n) { par.resize(n, -1); }

    bool is_root(int x) { return par[x] < 0; }

    int find(int x) {
        if (is_root(x)) return x;
        return par[x] = find(par[x]);
    }

    int size(int x) { return -par[find(x)]; }

    bool unite(int x, int y) {
        x = find(x);
        y = find(y);
        if (x == y) return false;
        if (size(x) < size(y)) swap(x, y);
        par[x] += par[y];
        par[y] = x;
        return true;
    }

    bool same(int x, int y) { return find(x) == find(y); }
} UF;

int n, m;
mint X, ans;
vector<pii> G[MAX];

int dfs(int v, int pv) {
    int ret = 1;
    for (auto [nv, z] : G[v]) {
        if (nv != pv) {
            int tmp = dfs(nv, v);
            ans += (X ^ z) * mint(tmp) * mint(n - tmp);
            ret += tmp;
        }
    }
    return ret;
}

int main() {
    scanf("%d%d%lld", &n, &m, &X.a);
    UF uf(n);
    rep(_, m) {
        int x, y, z;
        scanf("%d%d%d", &x, &y, &z);
        --x;
        --y;
        int sx = uf.size(x), sy = uf.size(y);
        if (uf.unite(x, y)) {
            G[x].emplace_back(y, z);
            G[y].emplace_back(x, z);
        }
    }
    dfs(0, -1);
    printf("%lld\n", ans.a);
}
0