ZOJ Monthly, May 2011

浙大5月1号月赛…… 我只A了4道……不爆成绩了。。。

 

 下面是某神队的解题报告：

 

 

 

ZOJ Monthly, May 2011
A	ZOJ3500	Electron Cloud	27.20% (139/511)
B	ZOJ3501	Roman Order	30.79% (400/1299)
C	ZOJ3502	Contest	15.11% (78/516)
D	ZOJ3503	Quadratic Surface	2.28% (4/175)
E	ZOJ3504	P-norm	49.43% (262/530)
F	ZOJ3505	Yet Another Set of Numbers	28.00% (42/150)
G	ZOJ3506	Cut the Tree	17.94% (21/117)
H	ZOJ3507	Fractal	13.20% (66/500)
I	ZOJ3508	The War	13.53% (226/1670)
J	ZOJ3509	Island Communication	4.24% (7/165)

 

 

一是准备给包括在去年参加过集训在内的xpies一次新手上路前热身的机会。二是因为神奇的原因，存有summer2009和summer2010数据的ZOJ服务器在省赛后被拔网线了。所以这次的题目来自校赛省赛未用题(ABCD)、summer2007(J)和summer2008(EFGHI)。就算题目总体看来水题很多啦，不过两个小时被圆满什么的＞＿＜，何方神圣啊orz(台大一队？)……不管怎样把五月份混过去了，下次Monthly大概在六月底了吧……

 

ZOJ3500. Electron Cloud

[math]

 求两个球的并的体积。

判断两个球的位置关系和圆没有区别，在相离、相切和包含的时候很简单，所以主要就是要考虑相交的情况。这时候的体积是两个球冠的体积，球冠的体积可以很简单的积分得到：

而要求h₁和h₂，也只要解：

 

vls特意把精度和数据都出得很低，所以数值积分也是行得通的。

 

#include <cmath> #include <cstdio> #include <algorithm> using namespace std; double sqr(double x) { return x * x; } double vol(double r, double h) { return M_PI * (r * r - h * h / 3) * h; } int main() { int re; double x1, y1, z1, r1, x2, y2, z2, r2, d, t, ans; scanf("%d", &re); for (int ri = 1; ri <= re; ++ri) { scanf("%lf%lf%lf%lf%lf%lf%lf%lf", &x1, &y1, &z1, &r1, &x2, &y2, &z2, &r2); d = sqrt(sqr(x1 - x2) + sqr(y1 - y2) + sqr(z1 - z2)); if (r1 < r2) { swap(r1, r2); } if (d >= r1 + r2) { ans = vol(r1, r1) - vol(r1, -r1) + vol(r2, r2) - vol(r2, -r2); } else if (d <= r1 - r2) { ans = vol(r1, r1) - vol(r1, -r1); } else { t = (r1 * r1 - r2 * r2) / d; z1 = (d + t) / 2; z2 = (d - t) / 2; ans = vol(r1, z1) - vol(r1, -r1) + vol(r2, z2) - vol(r2, -r2); } printf("%.2lf/n", ans); } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //744 2011-04-27 03:40:42 Accepted A C++ 0 180 watashi@ArcOfDream Source  

 

ZOJ3501. Roman Order

 [sorting]

将数字按它们的罗马数字表示排序。

将数字转成pair<roman, arabic>，然后直接sort就好了。最近才知道，原来这种常见的技术有个名字叫“Schwartzian Transform”。

 

#include <cstdio> #include <string> #include <vector> #include <cassert> #include <utility> #include <algorithm> using namespace std; const string ones [] = {"", "I", "II", "III", "IV", "V", "VI", "VII", "VIII", "IX"}; const string tens [] = {"", "X", "XX", "XXX", "XL", "L", "LX", "LXX", "LXXX", "XC"}; const string hundreds [] = {"", "C", "CC", "CCC", "CD", "D", "DC", "DCC", "DCCC", "CM"}; const string thousands[] = {"", "M", "MM", "MMM"}; string a2r(int x) { return thousands[x / 1000 % 10] + hundreds[x / 100 % 10] + tens[x / 10 % 10] + ones[x % 10]; } int main() { int re, n; vector<pair<string, int> > v; scanf("%d", &re); for (int ri = 1; ri <= re; ++ri) { assert(scanf("%d", &n) == 1); assert(1 <= n && n <= 10000); v.resize(n); for (int i = 0; i < n; ++i) { assert(scanf("%d", &v[i].second) == 1); v[i].first = a2r(v[i].second); } sort(v.begin(), v.end()); for (int i = 0; i < n; ++i) { printf("%d%c", v[i].second, i == n - 1 ? '/n' : ' '); } } assert(scanf("%d", &re) == EOF); return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //729 2011-04-22 17:03:25 Accepted B C++ 180 576 watashi@ArcOfDream Source

 

ZOJ3502. Contest

 [DP]

按顺序尝试10道题，每道题的AC概率与前面尝试过的题目有关，由给定的概率表决定。要求字典序最小的最优尝试顺序。

很显然的动态规划，2^10表示尝试过的题目，记录字典序最小的最优方案和对应的期望。注意比较浮点数时要使用eps。实现得好搜索应该也没问题。

 

#include <cstdio> #include <string> #include <utility> #include <algorithm> using namespace std; const double EPS = 1e-8; const int MAXN = 12; double p[MAXN][MAXN]; double dp[1 << MAXN]; string pre[1 << MAXN]; int main() { int re, n; double q; scanf("%d", &re); for (int ri = 1; ri <= re; ++ri) { scanf("%d", &n); for (int i = 0; i < n; ++i) { for (int j = 0; j < n; ++j) { scanf("%lf", &p[i][j]); } } dp[0] = 0; pre[0] = ""; for (int k = 1; k < (1 << n); ++k) { dp[k] = -1; for (int i = 0; i < n; ++i) { if (k & (1 << i)) { q = 0; for (int j = 0; j < n; ++j) { if (k & (1 << j)) { q = max(q, p[j][i]); } } q = dp[k ^ (1 << i)] + q / 100; if ((q > dp[k] + EPS) || (q > dp[k] - EPS && pre[k ^ (1 << i)] < pre[k])) { dp[k] = q; pre[k] = pre[k ^ (1 << i)]; pre[k] += ('A' + i); } } } } printf("%.2f/n%s/n", dp[(1 << n) - 1], pre[(1 << n) - 1].c_str()); } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //740 2011-04-27 03:24:50 Accepted C C++ 40 228 watashi@ArcOfDream Source

ZOJ3503. Quadratic Surface

 [math, 秩, 行列式, 一元三次方程]

是Conic Section的姊妹题，原计划是一起出在省赛的。作为“圆锥曲线”的加强版，“二次曲面”没有b=0这样的简化条件，要求一个一般的三元二次方程对应的二次曲面类型。不过题目中给出了判别的算法，所以只要实现就好了。

根据题目中所给的算法，需要求矩阵的秩(rank)和行列式(det)，这个可以用高斯消元法求得。求k比较麻烦，可以通过求对称矩阵的特征根(eig)或解一元三次方程的根求得。求三次方程的根可以用卡丹公式，不过我用了不同的算法。首先，我们知道三次方程一定有一个实根（事实上对于这题三个都是实根），且对于首相系数为1的三次方程，f(-inf)=-inf，f(inf)=inf，所以我们可以用数值的二分法求出一个实根。然后方程降为我们熟悉的一元二次方程，轻松求解。

这题的数据是通过先随机出一个二次曲面类型，再随机出对应的标准方程，最后通过随机的镜面、旋转和平移后去掉可能的ill-formed case生成的。生成器是ruby写的ZOJ3503gen.rb。

 

#include <cmath> #include <cstdio> #include <vector> #include <cassert> #include <numeric> #include <algorithm> using namespace std; const double EPS = 1e-6; struct Matrix { int n; double a[5][5]; double* operator[](int i) { return a[i]; } } e, E; int rank(Matrix a, double& d) { d = 1; for (int k = 0; k < a.n; ++k) { int ii = k, jj = k; for (int i = k; i < a.n; ++i) { for (int j = k; j < a.n; ++j) { if (fabs(a[i][j]) > fabs(a[ii][jj])) { ii = i; jj = j; } } } if (ii != k) { d = -d; for (int j = 0; j < a.n; ++j) { swap(a[k][j], a[ii][j]); } } if (jj != k) { d = -d; for (int i = 0; i < a.n; ++i) { swap(a[i][k], a[i][jj]); } } /* for (int i = 0; i < a.n; ++i) { for (int j = 0; j < a.n; ++j) { printf("%f ", a[i][j]); } puts(""); } puts(""); */ if (fabs(a[k][k]) < EPS) { d = 0; return k; } else { d *= a[k][k]; } for (int i = k + 1; i < a.n; ++i) { a[i][k] /= a[k][k]; for (int j = k + 1; j < a.n; ++j) { a[i][j] -= a[i][k] * a[k][j]; } } } return a.n; } inline double eval(int n, double c[], double x) { double y = c[0]; for (int i = 1; i < n; ++i) { y = y * x + c[i]; } return y; } bool gao(int n, double c[]) { double l = -1, r = 1, m; while (eval(n, c, l) > 0) { l *= 2; } while (eval(n, c, r) < 0) { r *= 2; } for (int loop = 0; loop < 64; ++loop) { // while (r - l > EPS) { m = (l + r) / 2; if (eval(n, c, m) < 0) { l = m; } else { r = m; } } m = (l + r) / 2; for (int i = 1; i < n; ++i) { c[i] += c[i - 1] * m; } assert(c[0] == 1); assert(fabs(c[n - 1]) < EPS); if (fabs(c[2]) < EPS) { return fabs(m) < EPS || fabs(c[1]) < EPS || m * c[1] < 0; } else { return c[2] > 0 && (fabs(m) < EPS || m * c[1] < 0); } } int main() { int re, rho3, rho4; double c[10], a[4], delta; Matrix m; scanf("%d", &re); for (int ri = 1; ri <= re; ++ri) { for (int i = 0; i < 10; ++i) { assert(scanf("%lf", &c[i]) == 1); assert(fabs(c[i]) <= 1000); } assert(fabs(c[0]) + fabs(c[1]) + fabs(c[2]) > 1e-4); m.n = 3; m[0][0] = c[0]; m[1][1] = c[1]; m[2][2] = c[2]; m[0][1] = m[1][0] = c[5]; m[0][2] = m[2][0] = c[4]; m[1][2] = m[2][1] = c[3]; rho3 = rank(m, delta); m.n = 4; m[3][3] = c[9]; m[0][3] = m[3][0] = c[6]; m[1][3] = m[3][1] = c[7]; m[2][3] = m[3][2] = c[8]; rho4 = rank(m, delta); a[0] = 1; a[1] = c[0] + c[1] + c[2]; a[2] = c[0] * c[1] + c[1] * c[2] + c[2] * c[0] - c[3] * c[3] - c[4] * c[4] - c[5] * c[5]; a[3] = c[0] * c[1] * c[2] + 2 * c[3] * c[4] * c[5] - c[0] * c[3] * c[3] - c[1] * c[4] * c[4] - c[2] * c[5] * c[5]; // printf("%d %d %lf/n", rho3, rho4, delta); switch (rho3) { case 1: switch (rho4) { case 1: puts("coincident planes"); break; case 2: puts("parallel planes (imaginary or real)"); break; case 3: puts("parabolic cylinder"); break; default: assert(false); } break; case 2: switch (rho4) { case 2: printf("intersecting planes (%s)/n", gao(4, a) ? "imaginary" : "real"); break; case 3: if (gao(4, a)) { puts("elliptic cylinder (imaginary or real)"); } else { puts("hyperbolic cylinder"); } break; case 4: puts(delta > 0 ? "hyperbolic paraboloid" : "elliptic paraboloid"); break; default: assert(false); } break; case 3: switch (rho4) { case 3: printf("elliptic cone (%s)/n", gao(4, a) ? "imaginary" : "real"); break; case 4: if (gao(4, a)) { printf("ellipsoid (%s)/n", delta > 0 ? "imaginary" : "real"); } else { puts(delta > 0 ? "hyperboloid of one sheet" : "hyperboloid of two sheets"); } break; default: assert(false); } break; default: assert(false); } } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //728 2011-04-22 17:02:29 Accepted D C++ 100 180 watashi@ArcOfDream Source  

 

 

ZOJ3504. P-norm

求两个N维复向量的差的P-范数。

dead do。比较麻烦的就是输入了。用std::complex和istream可以直接读入”(Re,Im)”格式的复数，所以如果用istringstream读入的话，是一件很轻松的事。

 

#include <cmath> #include <cstdio> #include <vector> #include <string> #include <sstream> #include <iostream> #include <complex> using namespace std; typedef complex<long double> Complex; typedef vector<Complex> Vector; Vector toV(const string& str) { istringstream iss(str); Vector ret; Complex c; while (iss >> c) { ret.push_back(c); } return ret; } int main(void) { string buf; long double p; long double ans; while (getline(cin, buf)) { Vector a, b; a = toV(buf); getline(cin, buf); b = toV(buf); cin >> p; while (cin.get() != '/n'); if (a.size() != b.size()) { throw 1; } ans = 0.0; for (size_t i = 0; i < a.size(); i++) { ans += powl(abs(a[i] - b[i]), p); } printf("%.16Le/n", powl(ans, 1.0 / p)); } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //733 2011-04-24 02:01:21 Accepted E C++ 0 188 watashi@ArcOfDream Source

 

 

ZOJ3505. Yet Another Set of Numbers

 [DP, counting]

略

比较经典的DP加计数问题模型，可以算这类问题里比较简单的一个例子。只要实现str2id和id2str两个函数就好了。DP只是简单的sum{3^i}，而计数部分，也就是str2id和id2str这两个函数，也和最经典的那种求第k个C(n, m)组合的算法一样。PS: 题目所描述的集合其实是相当于是给一个深度为n的三叉树的每个节点进行遍历标号。

 

#include <cstdio> #include <algorithm> using namespace std; const int MAXN = 30; long long dp[30]; long long id(int n, char c, char* s) { if (*s == 0) { return 0; } else { long long m = 1; for (char i = '0'; i < *s; ++i) { if (i == c) { continue; } m += dp[n - 1]; } m += id(n - 1, *s, s + 1); return m; } } void gao(int n, long long m, char c, char* s) { if (m == 0) { *s = '/0'; } else { --m; for (*s = '0'; ; ++*s) { if (*s == c) { continue; } if (m < dp[n - 1]) { gao(n - 1, m, *s, s + 1); break; } else { m -= dp[n - 1]; } } } } int main() { int n; long long m; char s[80]; dp[0] = 1; for (int i = 1; i < MAXN; ++i) { dp[i] = dp[i - 1] * 3; } for (int i = 1; i < MAXN; ++i) { dp[i] += dp[i - 1]; } while (scanf("%d%lld%s", &n, &m, s) != EOF) { m = id(n, '0', s) - m; gao(n, m, '0', s); puts(s); } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //734 2011-04-26 22:25:13 Accepted F C++ 60 180 watashi@ArcOfDream Source  

 

 

ZOJ3506. Cut the Tree

、 [DP]

对一棵树做恰好k次cut操作，每次可以cut一条边，然后丢掉其中一部分而得到一棵新的树。要求最后得到的树的最小/最大权值和。

比较经典的树形DP模型，其中每个节点对所有儿子又是一个二维DP。DP记录以该节点为根，cut了0 ~ k次所能得到的最小/最大权值。每次有两种选择，一种是保留child，这样update(dp[parent][a + b], dp[parent][a] + dp[child][b])。另一种是直接舍去这棵子树，这样update(dp[parent][a + b], dp[parent][a])，其中b可以取1 ~ sizeOfSubtree(j)，而不仅仅是1。最后的答案包括两种情况，一种是根在最优解中，那么就是dp[root][k]。否则假设i是最优解中离根最近的点，那么最优解是dp[i][k - j]，其中j可以取1 ~ (n – sizeOfSubtree(i))，同样j不仅仅可以取1。

附上更强大的测试数据

show source

#include <cstdio> #include <vector> #include <algorithm> #include <functional> using namespace std; const int MAXN = 1024; const int MAXM = 24; const int INF = 1000000007; int n, m, v[MAXN]; vector<int> e[MAXN]; int d[MAXN], dp[MAXN][MAXM]; void gao(int s, int p) { e[s].erase(remove(e[s].begin(), e[s].end(), p), e[s].end()); d[s] = 0; dp[s][0] = 0; fill(dp[s] + 1, dp[s] + 1 + m, INF); for (vector<int>::const_iterator t = e[s].begin(); t != e[s].end(); ++t) { gao(*t, s); d[s] += d[*t] + 1; for (int j = m; j >= 0; --j) { dp[s][j] += dp[*t][0]; for (int k = j - 1; k >= 0; --k) { dp[s][j] = min(dp[s][j], dp[s][k] + dp[*t][j - k]); } for (int k = 1; k <= d[*t] + 1 && k <= j; ++k) { // !!! dp[s][j] = min(dp[s][j], dp[s][j - k]); } } } transform(dp[s], dp[s] + 1 + m, dp[s], bind2nd(plus<int>(), v[s])); } int gao() { int y = INF; gao(0, -1); y = min(y, dp[0][m]); if (m > 0) { for (int i = 1; i < n; ++i) { for (int j = 1; j <= n - d[i] && j <= m; ++j) { // !!! y = min(y, dp[i][m - j]); } } } return y; } int main() { int a, b; while (scanf("%d%d", &n, &m) != EOF) { for (int i = 0; i < n; ++i) { scanf("%d", &v[i]); e[i].clear(); } for (int i = 1; i < n; ++i) { scanf("%d%d", &a, &b); --a; --b; e[a].push_back(b); e[b].push_back(a); } printf("%d ", gao()); for (int i = 0; i < n; ++i) { v[i] = -v[i]; } printf("%d/n", -gao()); } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //748 2011-04-27 04:57:21 Accepted G C++ 40 296 watashi@ArcOfDream Source  

 

 

ZOJ3507. Fractal

 [recursion]

打印分形，要求rstrip(trailing spaces)。

对于分形问题，递归什么的也是很经典的算法了。不能输出trailing spaces是稍微麻烦一点的地方，最简单的办法是先把答案存到一个string/char[]里，去掉trailing spaces后再输出。事实上也可以通过相似的递归算法的判断一个字符是不是这一行的最后一个非空格字符，这样就不需要用那么多内存了。

 

#include <cstdio> #include <cassert> // mem : 1024k int h, w; bool blk[12], end[12][12]; char str[12][12]; int row[12][12], col[12][12]; bool isblk(int x, int n) { int i; while(n > 1) { i = x / row[n][1]; if(blk[i]) return true; x -= row[n][i]; --n; } return blk[x]; } bool show(int x, int y, int n) { int i, j; bool ret = true; while(n > 1) { i = x / row[n][1]; j = y / col[n][1]; if(str[i][j] == ' ') { putchar(' '); return false; } else { x -= row[n][i]; y -= col[n][j]; if(!end[i][j]) ret = false; } --n; } putchar(str[x][y]); if(!end[x][y]) ret = false; return ret; } int main(void) { int ri = 0; int n, hh, ww; while(scanf("%d%d", &h, &w) != EOF) { fgets(str[0], sizeof(str[0]), stdin); for (int i = 0; i < h; i++) { fgets(str[i], sizeof(str[i]), stdin); assert(str[i][w] == '/n'); blk[i] = false; for (int j = 0; j < w; j++) end[i][j] = false; for (int j = w - 1; j >= -1; j--) { if(j == -1) { blk[i] = true; // blank line break; } if(str[i][j] != ' ') { end[i][j] = true; break; } } } scanf("%d", &n); row[0][h] = col[0][w] = 1; for (int i = 1; i <= n; i++) { for (int j = 0; j <= h; j++) row[i][j] = j * row[i - 1][h]; for (int j = 0; j <= w; j++) col[i][j] = j * col[i - 1][w]; } if(ri > 0) putchar('/n'); printf("Fractal #%d:/n", ++ri); hh = row[n][h]; ww = col[n][w]; for (int i = 0; i < hh; i++) { if(isblk(i, n)) { putchar('/n'); continue; } for (int j = 0; j < ww; j++) if(show(i, j, n)) { putchar('/n'); break; } } } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //738 2011-04-26 22:37:45 Accepted H C++ 230 180 watashi@ArcOfDream Source

 

 

ZOJ3508. The War

source code (ZOJ3508.cpp) [greedy]

n个士兵，每个士兵有可以装备重量在[a_i, b_i]的装备。然后有m个重量在w_i的装备。问最多可以装备多少士兵。

很经典的排序加贪心模型。首先按b排序，然后扫描，每次选取尽量小的满足a<w<b的装备装备该士兵，如果不存在，则–ans。正确性挺容易证明的，不再赘述。

 

#include <set> #include <cstdio> #include <utility> #include <algorithm> using namespace std; const int MAXN = 65536; pair<int, int> p[MAXN]; int main() { int n, m, t; multiset<int> s; multiset<int>::iterator it; while (scanf("%d%d", &n, &m) != EOF) { for (int i = 0; i < n; ++i) { scanf("%d%d", &p[i].second, &p[i].first); } s.clear(); for (int i = 0; i < m; ++i) { scanf("%d", &t); s.insert(t); } sort(p, p + n); t = 0; for (int i = 0; i < n; ++i) { it = s.lower_bound(p[i].second); if (it != s.end() && *it <= p[i].first) { ++t; s.erase(it); } } printf("%d/n", t); } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //745 2011-04-27 03:53:16 Accepted I C++ 30 1484 watashi@ArcOfDream Source  

 

 

ZOJ3509. Island Communication

 [graph]

有n个岛，初始没有桥，每次可以在i和j间建一座桥，拆一座桥，或查询i和j是否联通。点数n=500，操作数m=50000。i和j之间最多只会建一次桥。

由于点数只有500，其实这题不需要任何高级的数据结构就有一个O(n*m)的算法。注意到点数很少，而要求i和j是否联通，我们只要保留这幅图的生成森林就足够了，即控制边的个数为O(n)。对于查询操作，直接dfs一次即可。对于加边操作，如果i和j不联通，那么直接加上这条边就好了。否则，加上这条边后一定形成一个圈，如果我们把这个圈中最早被删去的边提前删掉，是不会影响后面操作的结果的，所以我们删掉这条边，以保证图中始终没有圈。对于删边操作，如果这个边已经被提前删掉了，那么什么也不做，否则简单删掉就好了。要知道哪条边先被删掉，可以先读入所有操作，预处理之。对于每次操作，复杂度均为O(n)，总的复杂度O(n*m)。

 

#include <cstdio> #include <vector> #include <algorithm> using namespace std; const int MAXN = 512; const int MAXM = 65536; vector<int> e[MAXN]; int d[MAXN][MAXN]; int xx, yy, zz; bool dfs(int p, int s, int t) { if (s == t) { return true; } else { for (vector<int>::const_iterator i = e[s].begin(); i != e[s].end(); ++i) { if (*i != p && dfs(s, *i, t)) { if (zz > d[s][*i]) { xx = s; yy = *i; zz = d[xx][yy]; } return true; } } return false; } } char c[MAXM]; int x[MAXM], y[MAXM]; int main() { int ri = 0; int n, m; while (scanf("%d%d", &n, &m) != EOF) { for (int i = 0; i < n; ++i) { e[i].clear(); fill(d[i], d[i] + n, MAXM); } for (int i = 0; i < m; ++i) { scanf(" %c%d%d", &c[i], &x[i], &y[i]); --x[i]; --y[i]; if (c[i] == 'D') { d[x[i]][y[i]] = d[y[i]][x[i]] = i; } } if (ri > 0) { puts(""); } printf("Case %d:/n", ++ri); for (int i = 0; i < m; ++i) { if (c[i] == 'I') { zz = MAXM; if (dfs(-1, x[i], y[i])) { if (zz < d[x[i]][y[i]]) { e[xx].erase(remove(e[xx].begin(), e[xx].end(), yy), e[xx].end()); e[yy].erase(remove(e[yy].begin(), e[yy].end(), xx), e[yy].end()); e[x[i]].push_back(y[i]); e[y[i]].push_back(x[i]); } } else { e[x[i]].push_back(y[i]); e[y[i]].push_back(x[i]); } } else if (c[i] == 'D') { e[x[i]].erase(remove(e[x[i]].begin(), e[x[i]].end(), y[i]), e[x[i]].end()); e[y[i]].erase(remove(e[y[i]].begin(), e[y[i]].end(), x[i]), e[y[i]].end()); } else { puts(dfs(-1, x[i], y[i]) ? "YES" : "NO"); } } } return 0; } //Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB) User Name Admin //747 2011-04-27 04:20:17 Accepted J C++ 590 1784 watashi@ArcOfDream Source  

 

 

神队啊……Orz。Orz。Orz。Orz。……