#include using namespace atcoder; #include #define all(a) (a.begin()),(a.end()) #define rall(a) (a.rbegin()),(a.rend()) using namespace std; using ll=long long; using ld=long double; using pll=pair; templatebool chmax(T &a, const T &b) {if(abool chmin(T &a, const T &b) {if(b> (istream& is, modint1000000007& x) { unsigned int t; is >> t; x=t; return is; } istream& operator >> (istream& is, modint998244353& x) { unsigned int t; is >> t; x=t; return is; } istream& operator >> (istream& is, modint& x) { unsigned int t; is >> t; x=t; return is; } ostream& operator << (ostream& os, const modint1000000007& x) { os << x.val(); return os; } ostream& operator << (ostream& os, const modint998244353& x) { os << x.val(); return os; } ostream& operator << (ostream& os, const modint& x) { os << x.val(); return os; } templateistream& operator >> (istream& is, pair& p) { is >> p.first >> p.second; return is; } templateostream& operator << (ostream& os, const pair& p) { os << p.first << " " << p.second; return os;} templateistream& operator >> (istream& is, vector& v) { for(T& x:v) is >> x; return is; } templateostream& operator << (ostream& os, const vector& v) {for(int i=0;i<(int)v.size();i++) {os << v[i] << (i+1 == v.size() ? "":" ");} return os;} template void pt() { std::cout << "\n"; } template void pt_rest() { std::cout << "\n"; } template void pt_rest(const T& first, const A&... rest) { std::cout << " " << first; pt_rest(rest...); } template void pt(const T& first, const A&... rest) { std::cout << first; pt_rest(rest...); } template void resize(vector& vec, const H head){ vec.resize(head); } template void resize(vector& vec, const H& head, const T ... tail){ vec.resize(head); for(auto& v: vec) resize(v, tail...); } template void fill(V& x, const T& val){ x = val; } template void fill(vector& vec, const T& val){ for(auto& v: vec) fill(v, val); } template void vin(istream& is, const int idx, vector& head){ is >> head[idx]; } template void vin(istream& is, const int idx, vector& head, T& ... tail){ vin(is >> head[idx], idx, tail...); } template void vin(istream& is, vector& head, T& ... tail){ for(int i=0; i<(int)head.size(); i++) vin(is, i, head, tail...); } template void vin(vector& head, T& ... tail){ vin(cin, head, tail...); } mapdivisor(ll n){mapres;for(ll i=1;i*i<=n;i++)if(n%i==0){res[i]++;res[n/i]++;}return res;} mapfactrization(ll n){mapres;for(ll i=2;i*i<=n;i++)while(n%i==0){res[i]++;n/=i;}if(n>1)res[n]++;return res;} static const ll inf=1e18+7; static const ld pi=acos(-1); static const ld eps=1e-10; using v1=vector; using v2=vector; using v3=vector; using v4=vector; using S=pll; S op(S a,S b){return max(a,b);} S e(){return {-inf,-inf};} class Edge { public: ll t,w; Edge(){} Edge(ll t,ll w): t(t),w(w){} }; v1 dijkstra(ll s,vector>g) { ll n=g.size(); v1 d(n,inf); priority_queue,greater<>>q; d[s]=0; q.push({0,s}); while(q.size()){ ll u=q.top().second; ll du=q.top().first; q.pop(); if(d[u]>N>>M>>K; vector>g(N); sets; for(ll i=0;i>r;r--;s.insert(r);} v2 es; for(ll i=0;i>a>>b>>c;a--;b--; if(a>b)swap(a,b); if(s.find(i)!=s.end())es.push_back({a,b,c}); g[a].push_back({b,c}); g[b].push_back({a,c}); } v3 d;resize(d,K,2,N); for(ll i=0;i>i&1))continue; for(ll j=0;j>j&1)continue; for(ll k=0;k<2;k++){ for(ll l=0;l<2;l++){ ll s=es[i][k]; ll w=es[j][2]; ll cost=min(d[j][l][s]+w*2,d[j][!l][s]+w); chmin(dp[S|1<