結果
問題 | No.1207 グラフX |
ユーザー | kyoprouno |
提出日時 | 2020-08-30 14:10:24 |
言語 | C++14 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 469 ms / 2,000 ms |
コード長 | 4,869 bytes |
コンパイル時間 | 2,466 ms |
コンパイル使用メモリ | 191,344 KB |
実行使用メモリ | 53,836 KB |
最終ジャッジ日時 | 2023-08-09 12:18:33 |
合計ジャッジ時間 | 20,351 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge14 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 393 ms
45,296 KB |
testcase_01 | AC | 397 ms
45,312 KB |
testcase_02 | AC | 391 ms
45,488 KB |
testcase_03 | AC | 391 ms
45,252 KB |
testcase_04 | AC | 393 ms
45,248 KB |
testcase_05 | AC | 455 ms
53,636 KB |
testcase_06 | AC | 459 ms
53,648 KB |
testcase_07 | AC | 462 ms
53,836 KB |
testcase_08 | AC | 310 ms
37,152 KB |
testcase_09 | AC | 333 ms
41,012 KB |
testcase_10 | AC | 457 ms
50,236 KB |
testcase_11 | AC | 469 ms
53,648 KB |
testcase_12 | AC | 302 ms
37,772 KB |
testcase_13 | AC | 184 ms
27,148 KB |
testcase_14 | AC | 375 ms
44,060 KB |
testcase_15 | AC | 335 ms
39,900 KB |
testcase_16 | AC | 180 ms
27,388 KB |
testcase_17 | AC | 259 ms
34,672 KB |
testcase_18 | AC | 218 ms
33,764 KB |
testcase_19 | AC | 262 ms
32,128 KB |
testcase_20 | AC | 381 ms
44,784 KB |
testcase_21 | AC | 48 ms
19,812 KB |
testcase_22 | AC | 255 ms
35,052 KB |
testcase_23 | AC | 278 ms
36,712 KB |
testcase_24 | AC | 195 ms
32,092 KB |
testcase_25 | AC | 388 ms
44,344 KB |
testcase_26 | AC | 300 ms
38,720 KB |
testcase_27 | AC | 360 ms
42,252 KB |
testcase_28 | AC | 341 ms
41,160 KB |
testcase_29 | AC | 337 ms
41,488 KB |
testcase_30 | AC | 179 ms
29,144 KB |
testcase_31 | AC | 159 ms
25,644 KB |
testcase_32 | AC | 149 ms
28,748 KB |
testcase_33 | AC | 163 ms
28,884 KB |
testcase_34 | AC | 326 ms
39,148 KB |
testcase_35 | AC | 50 ms
20,116 KB |
testcase_36 | AC | 325 ms
40,320 KB |
testcase_37 | AC | 282 ms
36,068 KB |
testcase_38 | AC | 84 ms
23,180 KB |
testcase_39 | AC | 170 ms
30,172 KB |
testcase_40 | AC | 139 ms
23,628 KB |
testcase_41 | AC | 233 ms
31,756 KB |
testcase_42 | AC | 19 ms
18,524 KB |
testcase_43 | AC | 18 ms
18,368 KB |
testcase_44 | AC | 20 ms
18,672 KB |
testcase_45 | AC | 370 ms
45,140 KB |
testcase_46 | AC | 360 ms
45,268 KB |
testcase_47 | AC | 366 ms
45,476 KB |
testcase_48 | AC | 372 ms
45,264 KB |
ソースコード
#pragma GCC optimize("O3") #include <bits/stdc++.h> #define ll long long #define rep(i,n) for(ll i=0;i<(n);i++) #define pll pair<ll,ll> #define pii pair<int,int> #define pq priority_queue #define pb push_back #define eb emplace_back #define fi first #define se second #define endl '\n' #define ios ios_base::sync_with_stdio(0),cin.tie(0),cout.tie(0); #define lb(c,x) distance(c.begin(),lower_bound(all(c),x)) #define ub(c,x) distance(c.begin(),upper_bound(all(c),x)) using namespace std; inline int topbit(unsigned long long x){ return x?63-__builtin_clzll(x):-1; } inline int popcount(unsigned long long x){ return __builtin_popcountll(x); } inline int parity(unsigned long long x){//popcount%2 return __builtin_parity(x); } template<class T> inline bool chmax(T& a,T b){if(a<b){a=b;return 1;}return 0;} template<class T> inline bool chmin(T& a,T b){if(a>b){a=b;return 1;}return 0;} const ll INF=1e15; const ll mod=1e9+7; struct mint { ll x; // typedef long long ll; mint(ll x=0):x((x%mod+mod)%mod){} mint operator-() const { return mint(-x);} mint& operator+=(const mint a) { if ((x += a.x) >= mod) x -= mod; return *this; } mint& operator-=(const mint a) { if ((x += mod-a.x) >= mod) x -= mod; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this;} mint operator+(const mint a) const { return mint(*this) += a;} mint operator-(const mint a) const { return mint(*this) -= a;} mint operator*(const mint a) const { return mint(*this) *= a;} mint pow(ll t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime mod mint inv() const { return pow(mod-2);} mint& operator/=(const mint a) { return *this *= a.inv();} mint operator/(const mint a) const { return mint(*this) /= a;} }; istream& operator>>(istream& is, mint& a) { return is >> a.x;} ostream& operator<<(ostream& os, const mint& a) { return os << a.x;} // combination mod prime // https://www.youtube.com/watch?v=8uowVvQ_-Mo&feature=youtu.be&t=1619 struct combination { vector<mint> fact, ifact; combination(ll n):fact(n+1),ifact(n+1) { assert(n < mod); fact[0] = 1; for (ll i = 1; i <= n; ++i) fact[i] = fact[i-1]*i; ifact[n] = fact[n].inv(); for (ll i = n; i >= 1; --i) ifact[i-1] = ifact[i]*i; } mint operator()(ll n, ll k) { if (k < 0 || k > n) return 0; return fact[n]*ifact[k]*ifact[n-k]; } mint p(ll n, ll k) { return fact[n]*ifact[n-k]; } } c(1000005); class DisjointSet{ public: vector<ll> rank,p,sz; DisjointSet(){} DisjointSet(ll size){ //作られうる木の頂点数の最大値を入れる。 rank.resize(size,0); p.resize(size,0); sz.resize(size,0); rep(i,size) makeSet(i); } void makeSet(ll x){ //xだけが属する木を作る。 p[x]=x; rank[x]=0; } bool same(ll x,ll y){ //xとyが同じ木に属するかどうか return findSet(x)==findSet(y); } void unite(ll x, ll y){ link(findSet(x),findSet(y)); } void link(ll x, ll y){ //rankが大きい方の根に小さい方の根をつける。 if(rank[x]>rank[y]){ p[y]=x; } else{ p[x]=y; if(rank[x]==rank[y]){ rank[y]++; //xとyの木のrankが同じであれば、統合するとrankが1増える。 } } sz[x]+=sz[y]; sz[y]=sz[x]; } ll findSet(ll x){ //xが属する木の根の番号を返す if(x!=p[x]){ p[x]=findSet(p[x]); } return p[x]; } }; vector<vector<pll>> to; vector<ll> memo; ll dfs(ll v,ll p=-1){ ll res=1; for(auto q:to[v]){ if(q.fi==p) continue; res+=dfs(q.fi,v); } return memo[v]=res; } int main(){ ll n,m; mint X; cin >> n >> m >> X; DisjointSet ds=DisjointSet(n+1); vector<pair<ll,pll>> edges(m); memo.resize(n); to=vector<vector<pll>>(n); vector<ll> not_need(m); rep(i,m){ cin >> edges[i].se.fi >> edges[i].se.se >> edges[i].fi; edges[i].se.fi--; edges[i].se.se--; if(ds.same(edges[i].se.fi,edges[i].se.se)){ not_need[i]=1; } else{ ds.unite(edges[i].se.fi,edges[i].se.se); } } rep(i,m){ if(not_need[i]) continue; to[edges[i].se.fi].push_back({edges[i].se.se,edges[i].fi}); to[edges[i].se.se].push_back({edges[i].se.fi,edges[i].fi}); } ll d=dfs(0LL); mint ans=0; rep(i,m){ if(not_need[i]) continue; mint val=X.pow(edges[i].fi); ll s=memo[edges[i].se.fi]; ll t=memo[edges[i].se.se]; if(s<t) swap(s,t); mint u=t; ans+=u*(n-u.x)*val; } cout << ans << endl; return 0; }