結果

問題 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
コンパイル時間 2,915 ms
コンパイル使用メモリ 232,964 KB
実行使用メモリ 37,240 KB
最終ジャッジ日時 2024-11-15 14:56:10
合計ジャッジ時間 49,634 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 842 ms
24,368 KB
testcase_01 AC 753 ms
24,360 KB
testcase_02 WA -
testcase_03 AC 1,283 ms
27,440 KB
testcase_04 WA -
testcase_05 AC 433 ms
37,208 KB
testcase_06 AC 365 ms
37,240 KB
testcase_07 WA -
testcase_08 AC 1,563 ms
20,796 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 AC 668 ms
14,488 KB
testcase_14 AC 1,338 ms
27,732 KB
testcase_15 WA -
testcase_16 WA -
testcase_17 AC 1,132 ms
18,108 KB
testcase_18 WA -
testcase_19 AC 1,160 ms
18,728 KB
testcase_20 WA -
testcase_21 WA -
testcase_22 AC 1,106 ms
18,092 KB
testcase_23 WA -
testcase_24 AC 654 ms
15,316 KB
testcase_25 AC 1,378 ms
26,980 KB
testcase_26 WA -
testcase_27 AC 1,553 ms
26,492 KB
testcase_28 WA -
testcase_29 AC 1,462 ms
24,772 KB
testcase_30 WA -
testcase_31 AC 425 ms
12,952 KB
testcase_32 AC 482 ms
12,920 KB
testcase_33 AC 553 ms
13,276 KB
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 WA -
testcase_39 AC 550 ms
13,824 KB
testcase_40 AC 115 ms
10,288 KB
testcase_41 AC 1,055 ms
17,272 KB
testcase_42 AC 2 ms
6,820 KB
testcase_43 AC 2 ms
6,816 KB
testcase_44 AC 1 ms
6,820 KB
testcase_45 WA -
testcase_46 WA -
testcase_47 WA -
testcase_48 AC 1,157 ms
27,688 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