某岛

… : "…アッカリ~ン . .. . " .. .
June 23, 2013

TCO 2013 Algorithm Round 3B

弱菜只能围观了。。)。。

250. SimilarSequences

Brief description:

… 两个序列被认为是相似的,当前仅当从两个序列中各删除一个元素后相等。(这里相等指相同位置上的元素相同。。
。。给定长度为 n 的初始序列 A 和上界 up。。。。问存在多少序列与之相似。。(每个元素取自 [1, up]…
( n <= 40 .. )

Analysis:

… 暴力美学。。。。
。枚举 A 序列中删除的元素。。再枚举要补上的元素位置。(用 -1 标记之。。)。。把这个 token 放进 set 里。。记之 h。。
那么答案就是 SZ(h) * up ?。。。

。。。。明显这样会重复计数。。。

。。。因此我们用 SI B(ALL(A)) 记录 A 数组中出现过的元素。。。
。。再建一个 seth2。(之前的现在记作 h1。。
。。用以保存那些所有可能会被重复计数的 token 。。(既在 h1 中原先用 -1标记的那个位置。。现在依次用 B 中的元素填上。。。

那么总的答案就是 SZ(h1) * (up - SZ(B)) + SZ(h2)。。。

//}/* .................................................................................................................................. */

class SimilarSequences {
public:
	int count(vector <int> A, int up) {

	    SI B(ALL(A)); set<VI> h1, h2; int n = SZ(A); REP(i, n){
	        VI AA = A; ERS(AA, i); REP(j, n){
	            INS(AA, j, -1); h1.insert(AA);
	            ECH(it, B) AA[j] = *it, h2.insert(AA);
	            ERS(AA, j);
	        }
	    }

	    return sum(pdt(SZ(h1), (up - SZ(B))), SZ(h2));
	}
};

450. AntlerSwapping

Brief description:

。。一对数被称为是平衡的。。如果它们的差 < 给定的 cap。
… 给定 n 对数。。问至少多少次交换可以使它们全部平衡。。(交换在某对数的任意一项和另一对数的任意一项之间进行。。。
..( n <= 16 .. )

Analysis:

。。很容易想偏。。。但是这个数据范围的话似乎应该是 SCDP?。。。

dp[s] 表示只考虑 s 集合中的数对的话。。要让他们平衡至少需要几次交换。
那么显然有…

..
	        REP_SS(i, s) checkMin(dp[s], dp[i] + dp[s^i]);
..

这里 REP_SS 是枚举子集。。(参见 [技巧]枚举子集的飘逸写法 by 柯神小天使

。。。考虑赋初值。。

。。对于单元素集合。。那么如果这两个元素的差 < cap 那么显然代价是 0 否则是无穷大。。。 。。对于双元素集合。。如果可以凑成 2 对的话。。。那么初值设为 1 次。。(。。当然也有可能是 0 次。。0 次的话将会由上面的转移得到。。。。 。。一般的。。对于 n 个元素的集合。。如果可以凑成 n 对的话。。那么初值设为 n-1。。(至少可以通过 n-1 次交换使之平衡。。。 [cpp collapse="false" firstline="1" highlight=""] //}/* .................................................................................................................................. */ int dp[1<<16]; int n; class AntlerSwapping { public: int getmin(vector <int> L, vector <int> R, int C) { n = SZ(L); FLC(dp, 0x3f); REP_1(s, _U(n)){ VI A; REP(i, n) if (_1(s, i)) A.PB(L[i]), A.PB(R[i]); SRT(A); int j; REP_N(j, SZ(A)){ if (A[j+1] - A[j] > C) break; ++j; } if (j == SZ(A)) dp[s] = count_bits(s) - 1; REP_SS(i, s) checkMin(dp[s], dp[i] + dp[s^i]); } return dp[_U(n)] == INF ? -1 : dp[_U(n)]; } }; [/cpp]

1000. ToastJumping

Brief description:

。。问从 (0, 0) 跳到 (x, y) 至少需要跳几步。。。
。。。跳跃只能在整点之间进行。。且每次跳跃的 欧几里得距离^2 不得超过 d。。

Analysis:

.. 一图流。。。

。。。k 步后的可达区域是凸的。。而且就是 1 步可达区域的那个凸包每个点乘以 k。。。因此算法就是构造出最里层的那个凸包。。然后枚举每条边。。checkMax 。。
。。实际根据对称性。。。构造出第一象限的一半就够了。。

const int N = 50;
VP P;

int f(int x0, int y0, int dd){
    if (!x0 && !y0) return 0;

    LL t = sqr((LL)x0)+sqr((LL)y0);

    CLR(P); int d = sqrt(dd), y = d, x = 0;
#define xx sqr(x)
#define yy sqr(y)
    while(1){
        while (xx+yy>dd) --y; Po p(x, y);
        while (SZ(P)>=2 && dett(P[SZ(P)-2], P.back(), p) >= 0) P.pop_back();
        P.PB(p); if (x>=y) break; ++x;
    }

    x0 = abs(x0), y0 = abs(y0); if (x0 > y0) swap(x0, y0);

    int k = 0; Po p0(x0, y0); FOR(i, 1, SZ(P)){
        Po p1(P[i]-P[i-1]); LL p = det(p1, p0), q = det(p1, P[i-1]);
        //cout << p << " " << q << endl;
        checkMax(k, int(ceil(p, q)));
    }

    return k;
}

class ToastJumping {
public:
    vector <int> minJumps(vector <int> x, vector <int> y, vector <int> d){
        VI res; int n = SZ(x); REP(i, n) res.PB(f(x[i], y[i], d[i]));
        return res;
    }
};

External link:

https://apps.topcoder.com/wiki/display/tc/TCO+2013+Round+3B
https://gist.github.com/lychees/5843231