結果
問題 | No.1207 グラフX |
ユーザー | kyoprouno |
提出日時 | 2020-08-30 14:10:24 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 446 ms / 2,000 ms |
コード長 | 4,869 bytes |
コンパイル時間 | 2,441 ms |
コンパイル使用メモリ | 195,732 KB |
実行使用メモリ | 53,948 KB |
最終ジャッジ日時 | 2024-04-27 07:23:03 |
合計ジャッジ時間 | 19,213 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 402 ms
45,592 KB |
testcase_01 | AC | 396 ms
45,464 KB |
testcase_02 | AC | 384 ms
45,620 KB |
testcase_03 | AC | 389 ms
45,596 KB |
testcase_04 | AC | 389 ms
45,588 KB |
testcase_05 | AC | 446 ms
53,904 KB |
testcase_06 | AC | 439 ms
53,736 KB |
testcase_07 | AC | 438 ms
53,948 KB |
testcase_08 | AC | 311 ms
37,304 KB |
testcase_09 | AC | 335 ms
41,236 KB |
testcase_10 | AC | 443 ms
50,588 KB |
testcase_11 | AC | 444 ms
53,908 KB |
testcase_12 | AC | 307 ms
37,948 KB |
testcase_13 | AC | 193 ms
27,412 KB |
testcase_14 | AC | 374 ms
44,416 KB |
testcase_15 | AC | 333 ms
39,872 KB |
testcase_16 | AC | 184 ms
27,996 KB |
testcase_17 | AC | 257 ms
34,768 KB |
testcase_18 | AC | 216 ms
34,008 KB |
testcase_19 | AC | 268 ms
32,432 KB |
testcase_20 | AC | 385 ms
45,076 KB |
testcase_21 | AC | 49 ms
20,112 KB |
testcase_22 | AC | 255 ms
35,104 KB |
testcase_23 | AC | 279 ms
37,156 KB |
testcase_24 | AC | 193 ms
32,256 KB |
testcase_25 | AC | 379 ms
44,828 KB |
testcase_26 | AC | 299 ms
38,592 KB |
testcase_27 | AC | 357 ms
42,316 KB |
testcase_28 | AC | 336 ms
41,312 KB |
testcase_29 | AC | 336 ms
41,576 KB |
testcase_30 | AC | 184 ms
29,428 KB |
testcase_31 | AC | 170 ms
25,804 KB |
testcase_32 | AC | 157 ms
29,004 KB |
testcase_33 | AC | 165 ms
28,984 KB |
testcase_34 | AC | 320 ms
39,444 KB |
testcase_35 | AC | 49 ms
20,176 KB |
testcase_36 | AC | 316 ms
40,488 KB |
testcase_37 | AC | 282 ms
36,140 KB |
testcase_38 | AC | 85 ms
23,380 KB |
testcase_39 | AC | 170 ms
30,256 KB |
testcase_40 | AC | 145 ms
23,816 KB |
testcase_41 | AC | 247 ms
31,880 KB |
testcase_42 | AC | 19 ms
19,028 KB |
testcase_43 | AC | 18 ms
18,964 KB |
testcase_44 | AC | 18 ms
18,876 KB |
testcase_45 | AC | 368 ms
45,496 KB |
testcase_46 | AC | 357 ms
46,088 KB |
testcase_47 | AC | 353 ms
45,472 KB |
testcase_48 | AC | 354 ms
45,500 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; }