#include using namespace std; using ll=long long; using ull=unsigned long long; using P=pair; templateusing minque=priority_queue,greater>; templatebool chmax(T &a,const T &b){return (abool chmin(T &a,const T &b){return (a>b?(a=b,true):false);} templateistream &operator>>(istream &is,pair&p){is>>p.first>>p.second;return is;} templateistream &operator>>(istream &is,tuple&a){is>>std::get<0>(a)>>std::get<1>(a)>>std::get<2>(a);return is;} templateistream &operator>>(istream &is,array&a){for(auto&i:a)is>>i;return is;} templateistream &operator>>(istream &is,vector &a){for(auto &i:a)is>>i;return is;} templatevoid operator++(pair&a,int n){a.first++,a.second++;} templatevoid operator--(pair&a,int n){a.first--,a.second--;} templatevoid operator++(vector&a,int n){for(auto &i:a)i++;} templatevoid operator--(vector&a,int n){for(auto &i:a)i--;} #define overload3(_1,_2,_3,name,...) name #define rep1(i,n) for(int i=0;i<(int)(n);i++) #define rep2(i,l,r) for(int i=(int)(l);i<(int)(r);i++) #define rep(...) overload3(__VA_ARGS__,rep2,rep1)(__VA_ARGS__) #define reps(i,l,r) rep2(i,l,r) #define all(x) x.begin(),x.end() #define pcnt(x) __builtin_popcountll(x) #define fin(x) return cout<<(x)<<'\n',static_cast(0) #define yn(x) cout<<((x)?"Yes\n":"No\n") #define uniq(x) sort(all(x)),x.erase(unique(all(x)),x.end()) template inline int fkey(vector&z,T key){return lower_bound(z.begin(),z.end(),key)-z.begin();} ll myceil(ll a,ll b){return (a+b-1)/b;} template auto vec(const int (&d)[n],const T &init=T()){ if constexpr (id(d,init)); else return init; } #ifdef LOCAL #include #define SWITCH(a,b) (a) #else #define debug(...) static_cast(0) #define debugg(...) static_cast(0) #define SWITCH(a,b) (b) templateostream &operator<<(ostream &os,const pair&p){os<>testcase; for(int i=0;i struct Edge{ int from,to; T weight; int index; Edge(int from_,int to_,T weight_=T(),int index_=-1):from(from_),to(to_),weight(weight_),index(index_){} Edge():from(-1),to(-1),weight(),index(-1){} friend std::ostream &operator<<(std::ostream &os,const Edge&e){ os<<'['; os<<"from:"< struct Graph{ private: int n; std::vector>edge; std::vector>g; std::vectorptr; bool directed; struct graph_range{ using iterator=typename std::vector>::iterator; iterator l,r; iterator begin()const{return l;} iterator end()const{return r;} int size()const{return r-l;} Edge &operator[](int i)const{return l[i];} }; struct const_graph_range{ using iterator=typename std::vector>::const_iterator; iterator l,r; iterator begin()const{return l;} iterator end()const{return r;} int size()const{return r-l;} const Edge &operator[](int i)const{return l[i];} }; public: Graph(int n_,bool dir_):n(n_),directed(dir_){} Graph():n(0){} Graph(int n_,bool dir_,const std::vector>&e):n(n_),directed(dir_),edge(e){build();} template void read(int m){ edge.reserve(m); for(int i=0;i>u>>v; T w; if constexpr(index)u--,v--; if constexpr(weighted)std::cin>>w; else w=1; edge.emplace_back(u,v,w,i); } build(); } void add_edge(int u,int v){ int id=edge.size(); edge.emplace_back(u,v,1,id); } void add_edge(int u,int v,T w){ int id=edge.size(); edge.emplace_back(u,v,w,id); } void add_edge(int u,int v,T w,int index){ edge.emplace_back(u,v,w,index); } void build(){ std::vectorcnt(n+1,0); for(const Edge&e:edge){ cnt[e.from+1]++; if(!directed)cnt[e.to+1]++; } for(int i=1;i<=n;i++)cnt[i]+=cnt[i-1]; ptr=cnt; g.resize(cnt[n]); for(const Edge&e:edge){ g[cnt[e.from]++]=e; if(!directed)g[cnt[e.to]++]=Edge(e.to,e.from,e.weight,e.index); } } void reverse(){ if(directed){ for(Edge&e:edge)std::swap(e.from,e.to); build(); } } inline void to_directed(){ directed=true; build(); } inline void to_undirected(){ directed=false; build(); } void reserve(int m){edge.reserve(m);} graph_range operator[](int i){return graph_range{g.begin()+ptr[i],g.begin()+ptr[i+1]};} const_graph_range operator[](int i)const{return const_graph_range{g.begin()+ptr[i],g.begin()+ptr[i+1]};} const Edge& get_edge(int i)const{return edge[i];} std::vector>get_edges()const{return edge;} inline bool is_directed()const{return directed;} inline int size()const{return n;} inline int edge_size()const{return edge.size();} typename std::vector>::iterator begin(){return edge.begin();} typename std::vector>::iterator end(){return edge.end();} typename std::vector>::const_iterator begin()const{return edge.begin();} typename std::vector>::const_iterator end()const{return edge.end();} }; #include #include template constexpr std::enable_if_t::digits<=32,int>msb(T n){return n==0?-1:31-__builtin_clz(n);} template constexpr std::enable_if_t<(std::numeric_limits::digits>32),int>msb(T n){return n==0?-1:63-__builtin_clzll(n);} template constexpr std::enable_if_t::digits<=32,int>lsb(T n){return n==0?-1:__builtin_ctz(n);} template constexpr std::enable_if_t<(std::numeric_limits::digits>32),int>lsb(T n){return n==0?-1:__builtin_ctzll(n);} template constexpr std::enable_if_t,T>floor_pow2(T n){return n==0?0:T(1)< constexpr std::enable_if_t,T>ceil_pow2(T n){return n<=1?1:T(1)<<(msb(n-1)+1);} template constexpr T safe_div(T a,T b){return a/b-(a%b&&(a^b)<0);} template constexpr T safe_ceil(T a,T b){return a/b+(a%b&&(a^b)>0);} using namespace std; template struct radix_heap{ static_assert(is_integral_v); private: int sz; T last; vector>h[1+numeric_limits::digits]; vector>p; public: radix_heap(int n):p(n,make_pair(0,-1)),sz(0),last(0){} inline bool empty()const{return sz==0;} int pop(){ if(h[0].empty()){ int b=0; while(h[b].empty())b++; last=h[b].back().first; for(int j=0;j<(int)h[b].size();j++)if(last>h[b][j].first)last=h[b][j].first; for(int j=0;j<(int)h[b].size();j++){ int i=h[b][j].second; int nb=msb(h[b][j].first^last)+1; p[i]=make_pair(nb,(int)h[nb].size()); h[nb].emplace_back(move(h[b][j])); } h[b].clear(); } --sz; int i=h[0].back().second; p[i].second=-1; h[0].pop_back(); return i; } void decrease_key(int i,T v){ if(p[i].second<0){ int b=msb(v^last)+1; ++sz; p[i]=make_pair(b,(int)h[b].size()); h[b].emplace_back(v,i); } else if(h[p[i].first][p[i].second].first>v){ int preb=p[i].first,nxtb=msb(v^last)+1; if(nxtb,std::nullptr_t> =nullptr> std::pair,std::vector>dijkstra(const Graph&g,int s=0){ int n=g.size(); std::vectordst(n,std::numeric_limits::max()); std::vectorpre(n,-1); dst[s]=0; radix_heapheap(n); heap.decrease_key(s,0); while(!heap.empty()){ int x=heap.pop(); for(const Edge&e:g[x]){ if(dst[e.to]>dst[e.from]+e.weight){ pre[e.to]=e.from; dst[e.to]=dst[e.from]+e.weight; heap.decrease_key(e.to,dst[e.to]); } } } return std::make_pair(dst,pre); } template,std::nullptr_t> =nullptr> std::pair,std::vector>dijkstra(const Graph&g,int s=0){ int n=g.size(); std::vectordst(n,std::numeric_limits::max()); std::vectorpre(n,-1); std::priority_queue,std::vector>,std::greater>>que; dst[s]=0; que.push({std::make_pair(0,s)}); while(!que.empty()){ auto [d,x]=que.top(); que.pop(); if(dst[x]!=d)continue; for(const Edge&e:g[x]){ T nxt=d+e.weight; if(dst[e.to]>nxt){ pre[e.to]=x; dst[e.to]=nxt; que.push(std::make_pair(dst[e.to],e.to)); } } } return std::make_pair(dst,pre); } void SOLVE(){ int n,m,p; ll y; cin>>n>>m>>p>>y; Graphg(n,false); g.read<1>(m); vectorstore(n,-1); rep(i,p){ int d,e; cin>>d>>e; d--; store[d]=e; } auto dst=dijkstra(g).first; ll ans=0; rep(i,n)if(dst[i]<=y&&store[i]!=-1){ chmax(ans,(y-dst[i])/store[i]); } cout<