結果

問題 No.1207 グラフX
ユーザー suzuken_wsuzuken_w
提出日時 2020-08-30 18:09:46
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 3,891 bytes
コンパイル時間 3,411 ms
コンパイル使用メモリ 231,800 KB
実行使用メモリ 37,056 KB
最終ジャッジ日時 2024-04-27 12:51:35
合計ジャッジ時間 58,109 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 961 ms
24,364 KB
testcase_01 AC 931 ms
24,488 KB
testcase_02 WA -
testcase_03 AC 1,484 ms
27,432 KB
testcase_04 WA -
testcase_05 AC 493 ms
37,040 KB
testcase_06 AC 401 ms
37,056 KB
testcase_07 WA -
testcase_08 AC 1,898 ms
20,664 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 AC 829 ms
14,360 KB
testcase_14 AC 1,662 ms
26,840 KB
testcase_15 WA -
testcase_16 WA -
testcase_17 AC 1,344 ms
17,984 KB
testcase_18 WA -
testcase_19 AC 1,430 ms
18,684 KB
testcase_20 WA -
testcase_21 WA -
testcase_22 AC 1,399 ms
18,088 KB
testcase_23 WA -
testcase_24 AC 757 ms
15,444 KB
testcase_25 AC 1,584 ms
27,240 KB
testcase_26 WA -
testcase_27 AC 1,778 ms
25,976 KB
testcase_28 WA -
testcase_29 AC 1,434 ms
24,772 KB
testcase_30 WA -
testcase_31 AC 516 ms
12,832 KB
testcase_32 AC 632 ms
13,044 KB
testcase_33 AC 744 ms
13,400 KB
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 WA -
testcase_39 AC 818 ms
13,696 KB
testcase_40 AC 133 ms
10,160 KB
testcase_41 AC 1,322 ms
17,272 KB
testcase_42 AC 2 ms
6,948 KB
testcase_43 AC 2 ms
6,944 KB
testcase_44 AC 2 ms
6,944 KB
testcase_45 WA -
testcase_46 WA -
testcase_47 WA -
testcase_48 AC 1,461 ms
27,684 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#pragma GCC optimize("O3")
#include <bits/stdc++.h> 
using namespace std;
using ll=long long;
using P=pair<ll,ll>;
template<class T> using V=vector<T>; 
#define fi first
#define se second
#define all(v) (v).begin(),(v).end()
const ll inf=(1e18);
//const ll mod=998244353;
const ll mod=1000000007;
ll GCD(ll a,ll b) {return b ? GCD(b,a%b):a;}
ll LCM(ll c,ll d){return c/GCD(c,d)*d;}
struct __INIT{__INIT(){cin.tie(0);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);}} __init;
template<class T> bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T> bool chmin(T &a, const T &b) { if (a>b) { a=b; return 1; } return 0; }
struct mint{
    using ull=unsigned long long int;
    ull v;
    mint(ll vv=0){s(vv%mod+mod);}
    mint& s(ull vv){
        v=vv<mod?vv:vv-mod;
        return *this;
    }
    //オーバーロード
    mint operator-()const{return mint()-*this;}//符号反転
    mint&operator+=(const mint&val){return s(v+val.v);}
    mint&operator-=(const mint&val){return s(v+mod-val.v);}
    mint&operator*=(const mint&val){
        v=ull(v)*val.v%mod;
        return *this;
    }
    mint&operator/=(const mint&val){return *this*=val.inv();}
    mint operator+(const mint&val){return mint(*this)+=val;}
    mint operator-(const mint&val){return mint(*this)-=val;}
    mint operator*(const mint&val){return mint(*this)*=val;}
    mint operator/(const mint&val){return mint(*this)/=val;}
    mint pow(ull n)const{
        mint res(1),x(*this);
        while(n){
            if(n&1)res*=x;
            x*=x;
            n>>=1;
        }
        return res;
    }
     mint inv()const{return pow(mod-2);}
    //拡張ユークリッドの互除法
    /* mint inv()const{
        int x,y;
        int g=extgcd(v,mod,x,y);
        assert(g==1);
        if(x<0)x+=mod;
        return mint(x);
    }*/
    friend ostream& operator<<(ostream&os,const mint&val){
        return os<<val.v;
    }//出力
    bool operator<(const mint&val)const{return v<val.v;}
    bool operator==(const mint&val)const{return v==val.v;}
    bool operator>(const mint&val)const{return v>val.v;}
};
struct edge{
    ll to,cost;
    edge(ll to,ll cost):to(to),cost(cost) {}
};
vector<vector<edge>> g;
vector<ll> d;
ll n,m;
ll mid;
V<bool> used;
bool dijkstra(ll s){
    for(int i=0;i<n;i++)d[i]=inf;
    used.assign(n,false);
    d[s]=0;
    priority_queue<P,vector<P>,greater<P>> q;
    q.push(P(0,s));//始点を0で置く
    while(!q.empty()){
        ll cur=q.top().se,val=q.top().fi;
        q.pop();
        if(used[cur]||d[cur]<val)continue;
        used[cur]=true;
        for(edge &e:g[cur]){
            if(e.cost>mid)continue;
            if(used[e.to])continue;
            if(d[e.to]>e.cost){
                d[e.to]=e.cost;
                q.push(P(d[e.to],e.to));
            }
        }
    }
    return (*max_element(all(d))<inf);
}
V<ll> get_path(ll t){//頂点tへの最短経路
    V<ll> path;
    for(;t!=-1;t=d[t]){
        path.emplace_back(t);
    }
    reverse(all(path));
    return path;
}
mint ans=0;
V<ll> cnt;
ll D;
void dfs(int cur,int par){
    for(edge &e:g[cur]){
        if(e.to==par)continue;
        if(d[e.to]!=e.cost)continue;
        dfs(e.to,cur);
    }
    for(edge &e:g[cur]){
        if(e.to==par)continue;
        if(d[e.to]!=e.cost)continue;
        ans+=mint(D).pow(e.cost)*mint((cnt[e.to])*(n-cnt[e.to]));
        cnt[cur]+=cnt[e.to];
    }
    cnt[cur]++;
}
int main(){
    cin>>n>>m>>D;
    cnt.assign(n,0);
    d.resize(n);
    g.resize(n);
    ll x,y,z;
    for(int i=0;i<m;i++){
        cin>>x>>y>>z;
        x--;y--;
        g[x].emplace_back(y,z);
        g[y].emplace_back(x,z);
    }
    ll  l=-1,r=(1e9)+7;
    while(r-l>1){
        mid=(l+r)/2;
        if(dijkstra(0))r=mid;
        else l=mid;
    }
    // for(int i=0;i<n;i++)cout<<d[i]<<" ";
    // cout<<"\n";
     dfs(0,-1);
    cout<<ans<<"\n";
}
0