Problem:
D：歪みの国 - Country In Distortion -

物語

アリスはすっかり退屈していました。いつも一緒に遊んでいる白ウサギが、トランプ城へ出かけているからです。（あーあ、こんなことなら一緒に出かければよかったわ。）とアリスは思いました。けれども、ここは歪みの国。そうやすやすと出かけては、ひどく疲れてしまいます。どういうことかというと—

白ウサギは困っていました。トランプ城への道中、ちょっとだけ道を間違えてしまったのです。（困ったなあ。一度戻って、道を良く知っているアリスに案内してもらうことにしよう。）と白ウサギは考えました。それにしても、白ウサギはとっても疲れているみたい。白ウサギは懐中時計を見ながら、

「なんだってこの国の道は、通る度に長さが変わったように感じるんだ！この道、さっきは100メートルくらいに感じたけど、今度は1キロはあるんじゃないか？それも、通り抜けるのにかかる時間はほとんど変わってないって言うのに！」

と叫んでいます。そう、ここは歪みの国。まさに、時々刻々と道が歪み続ける不思議な国なのです。
やっとの思いでアリスの家に着いた白ウサギが開口一番、

「ごめんよアリス。トランプ城まで、道案内をお願いできるかい。」

と言うと、退屈のあまりオフトゥンに突っ伏してしていたアリスは元気に飛び起きました。

「もちろんよ！一緒に行きましょう！」

嬉しい一方で、疲れた様子の白ウサギを見て（でも、あまりに疲れてしまうのは悪いわね。）と思ったアリスは、家の奥から歪みの国の地図を取り出してきました。

「さあ、トランプ城までの疲れない道のりを探すわ！」

問題

　n頂点とm辺からなる無向グラフが与えられる。各辺iの移動にかかる時間はt_iで表される。また、初めの体感移動速度としてv_0が与えられる。その後、辺を1本通る度に体感移動速度が変化する。辺をj本通った後の体感移動速度は、整数a, b, cをパラメータとして以下の式で与えられる。

v_j = (a \times v_{j-1} + b) mod c

辺iを速度v_jで通ると体感移動距離がt_i \times v_jとなる。なお、歪みの国では体感移動速度が0でも移動可能であり、そのときの体感移動距離は0である。

　今、頂点1から出発し頂点nを経由して頂点1に戻ってくるような経路を考える。考えうる経路の中で体感移動距離の総和を最小化せよ。

入力形式

n m
x_1 y_1 t_1
...
x_m y_m t_m
v_0
a b c

入力は全て整数からなる。
1行目には、グラフの頂点数と辺数を表すn,mがこの順で空白区切りで与えられる。
続くm行のうちi行目には、i番目の辺が接続する2頂点x_i,y_iと、通過するのにかかる時間t_iがこの順で空白区切りで与えられる。
m+2行目には、初速度を表すv_0が与えられる。
m+3行目には、速度の計算式のパラメータa,b,cがこの順に空白区切りで与えられる。

制約

各値は整数である。

2 \leq n \leq 500

1 \leq m \leq min\{5000, n(n - 1) / 2\}

各i=1,...,mについて1 \leq x_i < y_i \leq nである。また、1 \leq i < j \leq mを満たすi, jについて、x_i \neq x_jまたはy_i \neq y_jである。つまり自己ループおよび多重辺はない。

各i=1,...,mについて1 \leq t_i \leq 10^8

0 \leq v_0 \leq 49

0 < a < c, 0 \leq b < c, v_0 < c \leq 50

グラフは連結である。すなわち任意の2頂点間に経路があることが保証される。

出力形式

体感移動距離の総和の最小値を1行に出力せよ。

入力例1

4 4
1 2 3
2 3 1
2 4 5
3 4 2
5
4 5 6

出力例1

34

入力例2

6 5
1 2 1
2 3 100
3 4 2
4 5 200
5 6 1
3
4 5 9

出力例2

341

通過するのに時間がかかる辺 (2,3) と (4,5) を通る前に、短い時間で通れる道を往復して、なるべく速度を落とした方が体感移動距離は短い。

入力例3

10 9
1 2 74150933
2 3 13214145
3 4 45636411
4 5 90778512
5 6 85676138
6 7 60913535
7 8 43917556
8 9 47519690
9 10 56593754
33
37 41 47

出力例3

23049391759

体感移動距離が非常に長くなる場合もある。

Solution in C++14 (judge verdict: Accepted)

#include <bits/stdc++.h>
using namespace std;

using VI = vector<int>;
using VVI = vector<VI>;
using PII = pair<int, int>;
using LL = long long;
using VL = vector<LL>;
using VVL = vector<VL>;
using PLL = pair<LL, LL>;
using VS = vector<string>;

#define ALL(a)  begin((a)),end((a))
#define RALL(a) (a).rbegin(), (a).rend()
#define PB push_back
#define EB emplace_back
#define MP make_pair
#define SZ(a) int((a).size())
#define SORT(c) sort(ALL((c)))
#define RSORT(c) sort(RALL((c)))

#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define REP(i,n)  FOR(i,0,n)

#define FF first
#define SS second
template<class S, class T>
istream& operator>>(istream& is, pair<S,T>& p){
  return is >> p.FF >> p.SS;
}
template<class S, class T>
ostream& operator<<(ostream& os, const pair<S,T>& p){
  return os << p.FF << " " << p.SS;
}
template<class T>
void maxi(T& x, T y){
  if(x < y) x = y;
}
template<class T>
void mini(T& x, T y){
  if(x > y) x = y;
}


const double EPS = 1e-10;
const double PI  = acos(-1.0);
const LL MOD = 1e9+7;
const LL INF = 1e15;

struct Edge{
  LL to, cost;

  Edge(LL t, LL c): to(t), cost(c)
  {}
  bool operator>(const Edge& rhs) const{
	return cost > rhs.cost;
  }
  bool operator<(const Edge& rhs) const{
	return cost < rhs.cost;
  }

};
typedef vector< vector<Edge> > Graph;

void Dijkstra(const Graph& edges, vector<LL>& d, int s){
  using Type = PLL;
  priority_queue<Type, vector<Type>, greater<Type> > pq;
  fill(d.begin(), d.end(), INF);
  d[s] = 0;
  pq.push(MP(0,s));

  while(!pq.empty()){
	Type pii = pq.top(); pq.pop();
	int v = pii.second;
	if(d[v] < pii.first) continue;
	for(unsigned int i=0;i<edges[v].size();++i){
	  const Edge& e = edges[v][i];
	  if(d[e.to] > d[v] + e.cost){
		d[e.to] = d[v] + e.cost;
		pq.push(MP(d[e.to], e.to));
	  }
	}
  }
}

int idx(int u, int v, int flag){
  return (u * 50 + v) * 2 + flag;
}

int main(){
  cin.tie(0);
  ios_base::sync_with_stdio(false);

  int N, M;
  cin >> N >> M;

  VL xs(M), ys(M), ts(M);
  REP(i,M){
	cin >> xs[i] >> ys[i] >> ts[i];
	--xs[i];
	--ys[i];
  }

  LL V, A, B, C;
  cin >> V >> A >> B >> C;
	
  Graph G(500*50*2);
  REP(i,M){
	REP(k,C){
	  G[idx(xs[i], k, 0)].EB(idx(ys[i], (A*k+B)%C, 0), k*ts[i]);
	  G[idx(ys[i], k, 0)].EB(idx(xs[i], (A*k+B)%C, 0), k*ts[i]);
	  G[idx(xs[i], k, 1)].EB(idx(ys[i], (A*k+B)%C, 1), k*ts[i]);
	  G[idx(ys[i], k, 1)].EB(idx(xs[i], (A*k+B)%C, 1), k*ts[i]);
	}
  }
  REP(k,C)
	G[idx(N-1, k, 0)].EB(idx(N-1, k, 1), 0);

  VL d(500*50*2);
  Dijkstra(G, d, idx(0, V, 0));

  LL ans = INF;
  REP(k,C)
	mini(ans, d[idx(0, k, 1)]);
  cout << ans << endl;

  return 0;
}