3641: 货车运输
Time Limit: 20 Sec Memory Limit: 256 MBSec Special Judge
Submit: 64 Solved: 27
[Submit][Status][Discuss]
Description
saffah所在的国家一共有N个城市，通过N条双向道路连接，使得城市之间两两可以到达。第i条道路连接了A[i]与B[i]，其长度为L[i] km。
这些道路分为M个等级。第i条道路的等级为X[i]。在第j级道路上行驶，你的限速是V[j] km/h，并且每行驶1小时会收取W[j]元的费用。根据道路分级的意义，V[j]与W[j]都随着j单调变化，即对于1 ≤ j < M，满足V[j] > V[j+1], W[j] > W[j+1]。
现在有Q辆货车需要运送货物。第k辆货车要从S[k]前往T[k]，并且其最大速度为U[k] km/h。你需要对每辆货车求出其最小运送花费。
Input
第一行三个整数N, M, Q。
接下来N行，每行四个整数A[i], B[i], L[i], X[i]，描述了一条道路。
接下来M行，每行两个整数V[j], W[j]，描述了一个道路等级的限速和费用。
接下来Q行，每行三个整数S[k], T[k], U[k]，描述了一辆货车。
Output
输出Q行，即每辆货车的最小运送花费。
你可以输出任意多位的实数，只要与标准答案的相对或绝对误差不超过0.0001就算正确。
Sample Input
4 2 2
1 2 50 1
2 3 50 1
1 3 50 2
3 4 50 1
100 20
10 10
1 4 100
1 4 10
Sample Output
30
150
HINT
对于所有的数据，1 ≤ N,M,Q ≤ 100,000, 1 ≤ A[i],B[i],S[k],T[k] ≤ N, 1 ≤ X[i] ≤ M, 1 ≤ L[i],V[j],W[j],U[k] ≤ 500,000。
样例解释：
第一辆货车应该沿1→2→3→4行驶，共在1级道路上行驶1.5小时，收费30元。
第二辆货车应该沿1→3→4行驶，在2级道路上行驶5小时，在1级道路上行驶5小时，收费150元。
出题人真是调皮，上来就给8个数组。。
傻逼代码题。。调了蛮久（主要是数据。。找了半天。。CH各种坑啊。。）。。最后发现
rep(i, 1, cCnt) {
dep[cir[i]] = 1;
}
写成：
rep(i, 1, cCnt) {
dep[i] = 1;
}
23333
rk3的代码。。
1 #include <cstdio>
2 #include <cstring>
3 #include <algorithm>
4 #include <iostream>
5 #include <cmath>
6 #include <string>
7 #include <vector>
8 #include <map>
9 #include <set>
10 #include <queue>
11 #include <climits>
12 #include <cstdlib>
13 #include <ctime>
14 using namespace std;
15 namespace my_useful_tools {
16 #define rep(_i, _k, _j) for(int _i = _k; _i <= _j; ++_i)
17 #define reu(_i, _k, _j) for(int _i = _k; _i < _j; ++_i)
18 #define red(_i, _k, _j) for(int _i = _k; _j <= _i; --_i)
19 #define foreach(_i, _s) for(typeof(_s.begin()) _i = _s.begin(); _i != _s.end(); ++_i)
20 #define pb push_back
21 #define mp make_pair
22 #define ipir pair<int, int>
23 #define ivec vector<int>
24 #define clr(t) memset(t, 0, sizeof t)
25 #define pse(t, v) memset(t, v, sizeof t)
26 #define brl puts("")
27 #define file(x) freopen(#x".in", "r", stdin), freopen(#x".out", "w", stdout)
28 #define file_hza freopen("input.txt", "r", stdin), freopen("output.txt", "w", stdout);
29 #define file_gen(x) freopen(#x".in", "w", stdout);
30 const int INF = 0x3f3f3f3f;
31 typedef long long LL;
32 typedef long double DB;
33 inline void pc(char c) { putchar(c); }
34 template<class T> inline T gcd(T a, T b) { return b == 0 ? a : gcd(b, a % b); }
35 inline char gchar() { char ret = getchar(); for(; ret == ‘\n‘ || ret == ‘\r‘ || ret == ‘ ‘; ret = getchar()); return ret; }
36 template<class T> inline void fr(T&ret) { char c = ‘ ‘; int flag = 1; for(c = getchar(); c != ‘-‘ && !(‘0‘ <= c && c <= ‘9‘); c = getchar());
37 if(c == ‘-‘) flag = -1, ret = 0; else ret = c - ‘0‘; for(c = getchar(); ‘0‘ <= c && c <= ‘9‘; c = getchar()) ret = ret * 10 + c - ‘0‘;
38 ret = ret * flag;
39 }
40 inline int fr() { int x; fr(x); return x; }
41 template<class T> inline void fr(T&a, T&b) { fr(a), fr(b); } template<class T> inline void fr(T&a, T&b, T&c) { fr(a), fr(b), fr(c); }
42 template<class T> inline T fast_pow(T base, T index, T mod = 2147483647, T ret = 1) {
43 for(; index; index >>= 1, base = base * base % mod) if(index & 1) ret = ret * base % mod;
44 return ret;
45 }
46 };
47 using namespace my_useful_tools;
48
49 const int maxn = 1e5 + 100;
50
51
52 int e, h[maxn], to[maxn*2], nxt[maxn*2], dis[maxn*2], rnk[maxn*2];
53 void addEdge(int u, int v, int c, int d) {
54 nxt[e] = h[u], to[e] = v, rnk[e] = c, dis[e] = d, h[u] = e++;
55 nxt[e] = h[v], to[e] = u, rnk[e] = c, dis[e] = d, h[v] = e++;
56 }
57
58 int V[maxn], W[maxn];
59 int n, m, q;
60
61 int fa[maxn], pa[maxn][21], cir[maxn], cCnt, onCir[maxn], dep[maxn], owner[maxn], lnk[maxn];
62
63 void get_cir(int u) {
64 //printf("%d\n", u);
65 for (int i=h[u]; i !=-1; i=nxt[i]) {
66 if (i != fa[u]) {
67 if (dep[to[i]] > 0) {
68 cir[++cCnt] = u;
69 lnk[cCnt] = fa[u];
70 while (u != to[i]) {
71 u = to[fa[u]];
72 cir[++cCnt] = u;
73 lnk[cCnt] = fa[u];
74 }
75 lnk[cCnt] = i;
76 break;
77 }
78 fa[to[i]] = i^1;
79 dep[to[i]] = dep[u]+1;
80 get_cir(to[i]);
81 if (cCnt != 0) break;
82 }
83 }
84 }
85
86 const int maxNode = ((maxn<<2)+maxn*30)*2;
87 #define cn cmt_node
88 #define lc ch[0]
89 #define rc ch[1]
90 struct cmt_node {
91 cn*ch[2];
92 long double a, b;
93 cn() {
94 ch[0] = ch[1] = NULL;
95 a = b = 0.0;
96 }
97 } pool[maxNode], *loc=pool, *root, *trt[maxn], *crt[maxn];
98
99 void build(int l, int r, cn*&rt) {
100 if (rt == NULL)
101 rt = loc++;
102 if (l == r) return ;
103 int mid=(l+r)>>1;
104 build(l, mid, rt->lc);
105 build(mid+1, r, rt->rc);
106 }
107
108 int rk;
109 long double aa, ab;
110 void insert(int l, int r, cn*&rt, cn*pre) {
111 rt = loc++;
112 rt->a = pre->a+aa, rt->b = pre->b+ab;
113 if (l == r) return ;
114 int mid = (l+r)>>1;
115 if (rk <= mid) rt->rc=pre->rc, insert(l, mid, rt->lc, pre->lc);
116 else rt->lc=pre->lc, insert(mid+1, r, rt->rc, pre->rc);
117 }
118 long double queryA(int ql, int qr, int l, int r, cn*ra, cn*rb, cn*f) {
119 if (ql <= l && r <= qr) {
120 return ra->a + rb->a - 2.0*f->a;
121 }
122 int mid = (l+r)>>1;
123 long double ret = 0.0;
124 if (ql <= mid) ret += queryA(ql, qr, l, mid, ra->lc, rb->lc, f->lc);
125 if (mid < qr) ret += queryA(ql, qr, mid+1, r, ra->rc, rb->rc, f->rc);
126 return ret;
127 }
128 long double queryB(int ql, int qr, int l, int r, cn*ra, cn*rb, cn*f) {
129 if (ql <= l && r <= qr) {
130 return ra->b + rb->b - 2.0*f->b;
131 }
132 int mid = (l+r)>>1;
133 long double ret = 0.0;
134 if (ql <= mid) ret += queryB(ql, qr, l, mid, ra->lc, rb->lc, f->lc);
135 if (mid < qr) ret += queryB(ql, qr, mid+1, r, ra->rc, rb->rc, f->rc);
136 return ret;
137 }
138 long double queryA(int ql, int qr, int l, int r, cn*ra, cn*rb) {
139 if (ql <= l && r <= qr) {
140 return rb->a - ra->a;
141 }
142 int mid = (l+r)>>1;
143 long double ret = 0.0;
144 if (ql <= mid) ret += queryA(ql, qr, l, mid, ra->lc, rb->lc);
145 if (mid < qr) ret += queryA(ql, qr, mid+1, r, ra->rc, rb->rc);
146 return ret;
147 }
148 long double queryB(int ql, int qr, int l, int r, cn*ra, cn*rb) {
149 if (ql <= l && r <= qr) {
150 return rb->b - ra->b;
151 }
152 int mid = (l+r)>>1;
153 long double ret = 0.0;
154 if (ql <= mid) ret += queryB(ql, qr, l, mid, ra->lc, rb->lc);
155 if (mid < qr) ret += queryB(ql, qr, mid+1, r, ra->rc, rb->rc);
156 return ret;
157 }
158
159 void dfs(int u, int f, int top) {
160 /*
161 if (u == 93078) {
162 printf("93078 %d\n", top);
163 }
164 if (u == 14257) {
165 printf("14257 %d\n", top);
166 }
167 */
168 owner[u] = top;
169 //pa[u][0] = f;
170 for (int i=h[u]; i!=-1; i=nxt[i])
171 if (to[i] != f && !onCir[to[i]]) {
172 pa[to[i]][0] = u;
173 dep[to[i]] = dep[u]+1;
174 rk = rnk[i];
175 aa = (DB)dis[i]*W[rk];
176 ab = (DB)dis[i]/V[rk]*W[rk];
177 insert(1, m, trt[to[i]], trt[u]);
178 dfs(to[i], u, top);
179 }
180 }
181
182 int lca(int u, int v) {
183 if (dep[u] < dep[v]) swap(u, v);
184 red(i, 20, 0)
185 if (dep[pa[u][i]] >= dep[v])
186 u = pa[u][i];
187 if (u == v) return u;
188 red(i, 20, 0)
189 if (pa[u][i] != pa[v][i])
190 u=pa[u][i], v=pa[v][i];
191 return pa[u][0];
192 }
193
194 vector<int> rV;
195 vector<int>::iterator it;
196 int getRank(int u) {
197 it = upper_bound(rV.begin(), rV.end(), u);
198 if (it == rV.end()) return 1;
199 int v = it-rV.begin();
200 return m-v+1;
201 }
202 long double qMinCost(int s, int t, int u, int r) {
203 if (s == t) return 0;
204 int f = lca(s, t);
205 long double ret = 0.0;
206 if (1 <= r-1) ret += queryA(1, r-1, 1, m, trt[s], trt[t], trt[f])/u;
207 if (r <= m) ret += queryB(r, m, 1, m, trt[s], trt[t], trt[f]);
208 return ret;
209 }
210 long double qMinCost2(int s, int t, int u, int r) {
211 if (t <= s) return 0;
212 long double ret = 0.0;
213 if (1 <= r-1) ret += queryA(1, r-1, 1, m, crt[s], crt[t])/(long double)u;
214 if (r <= m) ret += queryB(r, m, 1, m, crt[s], crt[t]);
215 return ret;
216 }
217
218 int main() {
219 fr(n, m, q);
220 root = NULL;
221 build(1, m, root);
222 trt[0] = crt[0] = root;
223 rep(i, 1, n) trt[i] = NULL;
224 e = 0, pse(h, -1);
225 int a, b, c, d;
226 rep(i, 1, n) {
227 fr(a, b), fr(d, c);
228 addEdge(a, b, c, d);
229 }
230 rep(i, 1, m)
231 fr(V[i], W[i]);
232 pse(fa, -1);
233 dep[1] = 1;
234 get_cir(1);
235 rep(i, 1, cCnt)
236 onCir[cir[i]] = i;
237 dep[0] = 0;
238 crt[0] = root;
239 rep(i, 1, cCnt) {
240 dep[cir[i]] = 1;
241 trt[cir[i]] = root;
242 dfs(cir[i], 0, cir[i]);
243 crt[i] = NULL;
244 int p = (i==1?lnk[cCnt]:lnk[i-1]);
245 rk = rnk[p];
246 aa = (DB)dis[p]*W[rk];
247 ab = (DB)dis[p]*W[rk]/V[rk];
248 insert(1, m, crt[i], crt[i-1]);
249 }
250 rep(j, 1, 20)
251 rep(i, 1, n)
252 pa[i][j] = pa[pa[i][j-1]][j-1];
253
254 rep(i, 1, m)
255 rV.pb(V[i]);
256 reverse(rV.begin(), rV.end());
257 int s, t, u, r;
258 while (q--) {
259 fr(s, t, u);
260 long double ans = 0.0;
261 r = getRank(u);
262 if (owner[s] == owner[t]) {
263 ans = qMinCost(s, t, u, r);
264 } else {
265 if (onCir[owner[t]] < onCir[owner[s]])
266 swap(s, t);
267 ans += qMinCost(s, owner[s], u, r);
268 ans += qMinCost(t, owner[t], u, r);
269 long double x = qMinCost2(onCir[owner[s]], onCir[owner[t]], u, r);
270 long double y = qMinCost2(1, onCir[owner[s]], u, r);
271 long double z = qMinCost2(onCir[owner[t]], cCnt, u, r);
272 long double w = qMinCost2(0, 1, u, r);
273 ans += min(x, y+z+w);
274 }
275 printf("%lf\n", (double)ans);
276 }
277
278 return 0;
279 }
时间： 2024-08-26 05:55:27