Language: C++

#include<bits/stdc++.h>
//using namespace std;
#define rep(i,j,n) for(ll i=(ll)(j);i<(ll)(n);i++)
#define REP(i,j,n) for(ll i=(ll)(j);i<=(ll)(n);i++)
#define per(i,j,n) for(ll i=(ll)(j);(ll)(n)<=i;i--)
#define ll long long
#define ALL(a) (a).begin(),(a).end()
#define disup(A,key) distance(A.begin(),upper_bound(ALL(A),(ll)(key)))
#define dislow(A,key) distance(A.begin(),lower_bound(ALL(A),(ll)(key)))
#define pb push_back
#define mp std::make_pair
#define endl "\n"
using std::cin;
using std::cout;
using std::vector;
using std::string;
using std::upper_bound;
using std::lower_bound;
using vi=vector<ll>;
using vii=vector<vi>;
using pii=std::pair<ll,ll>;
const ll MOD=1e9+7;
//const ll MOD=998244353;
//const ll MOD=100000;
const ll MAX=1e7;
const ll INF=(1ll<<60);
template<class T>
class prique :public std::priority_queue<T, std::vector<T>, std::greater<T>> {};
struct Binary_indexed_tree{
    int N;
    vi bit;
    Binary_indexed_tree(int n):N(n){
        bit.resize(N+1,0);
    }
    void add(int x,int a){
        for(x;x<=N;x+=(x&-x)) bit[x]+=a;
    }
    ll sum(int x){
        ll ret=0;
        for(x;x>0;x-=(x&-x)) ret+=bit[x];
        return ret;
    }
};
struct Union_Find{
    ll N;
    vi par;
    vi siz;
    Union_Find(int n):N(n){
        par.resize(N);
        siz.resize(N,1);
        rep(i,0,N) par[i]=i;
    }
    ll root(ll X){
        if(par[X]==X) return X;
        return par[X]=root(par[X]);
    }
    bool same(ll X,ll Y){
        return root(X)==root(Y);
    }
    void unite(ll X,ll Y){
        X=root(X);
        Y=root(Y);
        if(X==Y) return;
        par[X]=Y;
        siz[Y]+=siz[X];
        siz[X]=0;
    }
    ll size(ll X){
        return siz[root(X)];
    }
};
ll modpow(ll X,ll Y){
    ll sum=X,cnt=1;
    vi A;
    while(cnt<=Y){
        A.pb(sum);
        sum*=sum;
        sum%=MOD;
        cnt*=2;
    }
    int M=A.size();
    ll ret=1;
    REP(i,1,M){
        if(Y>=(1ll<<M-i)){
            Y-=(1ll<<M-i);
            ret*=A[M-i];
            ret%=MOD;
        }
    }
    return ret;
}
ll fac[MAX],finv[MAX],inv[MAX];
void COMinit() {
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;
    for (int i = 2; i < MAX; i++){
        fac[i] = fac[i - 1] * i % MOD;
        inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD;
        finv[i] = finv[i - 1] * inv[i] % MOD;
    }
}
ll COM(ll n,ll r){
    if(n<r||n<0||r<0) return 0;
    return fac[n]*finv[r]%MOD*finv[n-r];
}
signed main(){
    ll N,M; cin>>N>>M;
    vector<vector<pii>> edge(N);
    vi p(2);
    rep(i,0,2) cin>>p[i];
    rep(i,0,M){
        ll X,Y,Z; cin>>X>>Y>>Z;
        edge[X-1].pb(mp(Y-1,Z));
        edge[Y-1].pb(mp(X-1,Z));
    }
    vii dist(2,vi(N,INF));
    rep(i,0,2){
        prique<pii> que;
        que.push(mp(0,p[i]-1));
        dist[i][p[i]-1]=0;
        while(!que.empty()){
            ll X=que.top().second,Y=que.top().first;
            que.pop();
            rep(j,0,edge[X].size()){
                ll P=edge[X][j].first,Q=edge[X][j].second+Y;
                if(dist[i][P]>Q){
                    dist[i][P]=Q;
                    que.push(mp(Q,P));
                }
            }
        }
    }
    vii dp(2,vi(N));
    rep(i,0,2){
        dp[i][p[i]-1]=1;
        vector<pii> C;
        rep(j,0,N){
            if(dist[0][j]+dist[1][j]==dist[0][p[1]-1] C.pb(mp(dist[i][j],j));
        }
        sort(ALL(C));
        rep(j,0,C.size()){
            ll X=C[j].second,Y=C[j].first;
            rep(k,0,edge[X].size()){
                ll Z=edge[X][k].first;
                if(dist[0][Z]+dist[1][Z]==dist[0][p[1]-1]&&dist[i][Z]-dist[i][X]==edge[X][k].second){
                    dp[i][Z]+=dp[i][X];
                    dp[i][Z]%=MOD;
                }
            }
        }
    }
    ll ans=dp[0][p[1]-1]*dp[1][p[0]-1]%MOD;
    rep(i,0,N){
        if(dist[0][i]+dist[1][i]==dist[0][p[1]-1]&&dist[0][i]==dist[1][i]){
            ans=(ans+MOD-dp[0][i]*dp[1][i]%MOD)%MOD;
        }
    }
    rep(i,0,N){
        if(dist[0][i]+dist[1][i]!=dist[0][p[1]-1]) continue;
        if(dist[0][i]>=dist[1][i]) continue;
        rep(j,0,edge[i].size()){
            ll X=edge[i][j].first;
            if(dist[0][X]+dist[1][X]==dist[0][p[1]-1]){
                if(dist[1][X]*2<dist[0][p[1]-1]&&dist[0][X]-dist[0][i]==edge[i][j].second){
                    ans=(ans+MOD-dp[0][i]*dp[1][X]%MOD)%MOD;
                }
            }
        }
    }
    cout<<ans<<endl;
}


Compiler output:
a.cc: In function 'int main()':
a.cc:139:54: error: expected ';' before 'C'
  139 |             if(dist[0][j]+dist[1][j]==dist[0][p[1]-1] C.pb(mp(dist[i][j],j));
      |                                                      ^~
      |                                                      ;
a.cc:140:9: error: expected primary-expression before '}' token
  140 |         }
      |         ^
a.cc:139:78: error: expected ')' before '}' token
  139 |             if(dist[0][j]+dist[1][j]==dist[0][p[1]-1] C.pb(mp(dist[i][j],j));
      |               ~                                                              ^
      |                                                                              )
  140 |         }
      |         ~                                                                     
a.cc:140:9: error: expected primary-expression before '}' token
  140 |         }
      |         ^