[题解]CLYZ2018省选训(bao)练(zha)模拟赛 Day 8

最新推荐文章于 2019-03-26 21:47:17 发布

xyz32768

最新推荐文章于 2019-03-26 21:47:17 发布

阅读量262

点赞数

分类专栏：题解文章标签： orz

本文链接：https://blog.csdn.net/xyz32768/article/details/79367801

版权

题解专栏收录该内容

29 篇文章 0 订阅

订阅专栏

竟然卡SPFA，出题人好毒啊。

题目

T1：BZOJ 2330 / SCOI 2011 糖果 (candy)（差分约束）
T2：BZOJ 2743 / HEOI 2012 采花 (flower)（离线+树状数组）
T3：BZOJ 3294 / CQOI 2011 放棋子 (chess)（DP+组合数学+容斥）

T1

分析

差分约束裸题。对于一个条件：
$1\>A\>B$ ：建边<script type="math/tex" id="MathJax-Element-2"> </script>。
$2\>A\>B$ ：建边 <script type="math/tex" id="MathJax-Element-4"> </script>。
$3\>A\>B$ ：建边 <script type="math/tex" id="MathJax-Element-6"> </script>。
$4\>A\>B$ ：建边 <script type="math/tex" id="MathJax-Element-8"> </script>。
$5\>A\>B$ ：建边 <script type="math/tex" id="MathJax-Element-10"> </script>。
然后新建一个虚拟点 $N+1$ ，向每个点 $i$ 建边 $<N+1,i,1>$ 。
然后从 $N+1$ 开始SPFA求解单源 最长路， $N+1$ 到其他每个点的最长路之和就是答案。
无解的情况：存在正权环（一个点入队超过 $N$ 次）。
但是SPFA被卡掉了一个点，把建边的顺序稍微改了一下之后过了……坑啊！

Source

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
inline int read() {
    int res = 0; bool bo = 0; char c;
    while (((c = getchar()) < '0' || c > '9') && c != '-');
    if (c == '-') bo = 1; else res = c - 48;
    while ((c = getchar()) >= '0' && c <= '9')
        res = (res << 3) + (res << 1) + (c - 48);
    return bo ? ~res + 1 : res;
}
const int N = 1e5 + 5, M = 8e5 + 5, P = 7e6 + 5;
int n, m, ecnt, nxt[M], adj[N], go[M], dis[N], val[M], len, que[P], cnt[N];
bool vis[N]; ll ans;
void add_edge(int u, int v, int w) {
    nxt[++ecnt] = adj[u]; adj[u] = ecnt; go[ecnt] = v; val[ecnt] = w;
}
void orzcyxdalao(int S) {
    int i; cnt[que[len = 1] = S] = 1;
    memset(dis, -1, sizeof(dis)); dis[S] = 0;
    for (i = 1; i <= len; i++) {
        int u = que[i]; vis[u] = 0;
        for (int e = adj[u], v; e; e = nxt[e])
            if (dis[u] + val[e] > dis[v = go[e]]) {
                dis[v] = dis[u] + val[e];
                if (!vis[v]) {
                    vis[que[++len] = v] = 1;
                    if ((++cnt[v]) == n) puts("-1"), exit(0);
                }
            }
    }
}
int main() {
    int i, op, x, y; n = read(); m = read(); while (m--) {
        op = read(); x = read(); y = read();
        if (op == 1) add_edge(x, y, 0), add_edge(y, x, 0);
        else if (op == 2) add_edge(x, y, 1);
        else if (op == 3) add_edge(y, x, 0);
        else if (op == 4) add_edge(y, x, 1);
        else if (op == 5) add_edge(x, y, 0);
        if (!(op & 1) && x == y) return puts("-1"), 0;
    }
    for (i = n; i; i--) add_edge(n + 1, i, 1); orzcyxdalao(n + 1);
    for (i = 1; i <= n; i++) ans += dis[i]; cout << ans << endl; return 0;
}

T2

分析

对于每一个 $i$ ，预处理出两个值：
$pre_i$ 表示在 $i$ 的左边，颜色与第 $i$ 朵花的颜色相同并且与 $i$ 距离最近的位置。
$suf_i$ 表示在 $i$ 的右边，颜色与第 $i$ 朵花的颜色相同并且与 $i$ 距离最近的位置。
如样例中 $pre_3=2,suf_1=5$ 。
询问要求颜色数，因此可以对于询问区间 $[l,r]$ ，只考虑在 $[l,r]$ 里 第一次出现对应颜色的位置，也就是满足 $pre_i<l$ 的 $i$ 。又由于要求出现次数至少 $2$ 次，这个条件也就是，对于一个位置 $i$ （ $i$ 的颜色在 $[l,r]$ 内为第一次出现）， $i$ 的后面（到 $r$ 为止）还有与 $i$ 的颜色相同的花，即 $suf_i\leq r$ 。
因此，问题转化为：
每次询问 $[l,r]$ 区间内，有多少个 $i$ 同时满足 $pre_i\leq l-1$ 和 $suf_i\leq r$ 。
但这样还是不好做，所以把一个询问拆成两个： $[1,l-1]$ 和 $[1,r]$ 。
分析后可以发现，对于任意的 $suf_i\leq r$ 都一定有 $i\in[1,r]$ ，对于任意的 $i\in[1,l-1]$ 都一定有 $pre_i\leq l-1$ ，所以现在询问有两种：
1、询问有多少个 $i$ 同时满足 $pre_i\leq l-1$ 和 $suf_i\leq r$ 。
2、询问有多少个 $i$ 同时满足 $i\leq l-1$ 和 $suf_i\leq r$ 。
就第一个询问，可以把一朵花看成一个点 $(pre_i,suf_i)$ ，每次询问 $(l-1,r)$ 的左下角有多少个点。可以将所有的花按照 $pre_i$ 排序，所有的 $[1,r]$ 式的询问也按照对应的 $l$ 排序，然后用两个指针从左往右分别扫描花和询问，按顺序插入每个点，确保当前已经插入的点的 $x$ 坐标一定小于当前询问的 $l-1$ ，然后就可以用树状数组查询了。第二个询问基本相似。
注意树状数组修改时下标为 $0$ 会死循环，因此要将坐标值加 $1$ 。

Source

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
inline int read() {
    int res = 0; bool bo = 0; char c;
    while (((c = getchar()) < '0' || c > '9') && c != '-');
    if (c == '-') bo = 1; else res = c - 48;
    while ((c = getchar()) >= '0' && c <= '9')
        res = (res << 3) + (res << 1) + (c - 48);
    return bo ? ~res + 1 : res;
}
const int N = 1e6 + 5;
int n, c, m, a[N], nxt[N], pre[N], suf[N], T[N << 1], L[N],
R[N], ans[N], w[N];
struct cyx {
    int x, y, id;
    cyx() {}
    cyx(int _x, int _y, int _id) : x(_x), y(_y), id(_id) {}
} fl[N], qn[N], tmp[N];
void change(int x, int v) {
    for (; x <= 2000000; x += x & -x) T[x] += v;
}
int ask(int x) {
    int res = 0; for (; x; x -= x & -x) res += T[x]; return res;
}
int main() {
    int i, j; n = read(); c = read(); m = read();
    for (i = 1; i <= n; i++) a[i] = read();
    for (i = 1; i <= c; i++) nxt[i] = 0;
    for (i = 1; i <= n; i++) pre[i] = nxt[a[i]], nxt[a[i]] = i;
    for (i = 1; i <= c; i++) nxt[i] = n + 1;
    for (i = n; i; i--) suf[i] = nxt[a[i]], nxt[a[i]] = i;
    for (i = 1; i <= m; i++) L[i] = read(), R[i] = read();
    for (i = 1; i <= n; i++) fl[i] = cyx(pre[i] + 1, suf[i] + 1, i);
    for (i = 1; i <= m; i++) qn[i] = cyx(L[i], R[i] + 1, i);
    for (i = 1; i <= n; i++) w[fl[i].x]++;
    for (i = 2; i <= n; i++) w[i] += w[i - 1];
    for (i = 1; i <= n; i++) tmp[w[fl[i].x]--] = fl[i];
    for (i = 1; i <= n; i++) fl[i] = tmp[i];
    memset(w, 0, sizeof(w));
    for (i = 1; i <= m; i++) w[qn[i].x]++;
    for (i = 2; i <= n; i++) w[i] += w[i - 1];
    for (i = 1; i <= m; i++) tmp[w[qn[i].x]--] = qn[i];
    for (i = 1; i <= m; i++) qn[i] = tmp[i];
    fl[n + 1] = cyx(0x3f3f3f3f, 0x3f3f3f3f, 0x3f3f3f3f);
    for (i = j = 1; i <= m; i++) {
        while (fl[j].x <= qn[i].x) change(fl[j].y, 1), j++;
        ans[qn[i].id] = ask(qn[i].y);
    }
    memset(T, 0, sizeof(T));
    for (i = j = 1; i <= m; i++) {
        while (j + 1 <= qn[i].x) change(suf[j] + 1, 1), j++;
        ans[qn[i].id] -= ask(qn[i].y);
    }
    for (i = 1; i <= m; i++) printf("%d\n", ans[i]);
    return 0;
}

T3

分析

注意到条件「不同颜色的棋子不能放在同一行」，可以得出，一种颜色的棋子放完之后，这些棋子会占领一定的行和列，这些行和列不能继续放棋子。
所以定义状态 $f[k][i][j]$ 表示到第 $k$ 种颜色，还剩下 $i$ 行 $j$ 列没被占领。
边界为 $f[0][n][m]=1$ 。
转移为：

f [k] [i] [j] = \sum i < h \leq n \sum j < l \leq m f [k - 1] [h] [l] \times C i h \times C j l \times g [k] [h - i] [l - j]

$f[k][i][j]=\sum_{i<h\leq n}\sum_{j<l\leq m}f[k-1][h][l]\times C_h^i\times C_l^j\times g[k][h-i][l-j]$

g[k][i][j] g [ k ] [ i ] [ j ] $g[k][i][j]$ 表示第

k k $k$ 种颜色恰好占领

i

$i$ 行

j j $j$ 列的方案数。
关于

g

$g$ ， ~~一开始我想出了一个 $250\times n^2\times m^2$ 的求法然后TLE，~~可以用组合数+容斥原理求得。
先不考虑「恰好」，则为

Ccntki×j C i × j c n t k $C_{i\times j}^{cnt_k}$ （

cntk c n t k $cnt_k$ 为第

k k $k$ 种颜色的棋子个数）
然后除去不合法的方案数，即

g [k] [i] [j] = C_{i \times j}^{c n t_{k}} - \sum_{h = 1}^{i} \sum_{l = 1}^{j} g [k] [h] [l] \times C_{i}^{h} \times C_{j}^{l} [h \neq i or l \neq j]

$g[k][i][j]=C_{i\times j}^{cnt_k}-\sum_{h=1}^i\sum_{l=1}^jg[k][h][l]\times C_i^h\times C_j^l[h\neq i\mbox{ or }l\neq j]$
最后结果就是

∑ni=0∑mj=0f[c][i][j] ∑ i = 0 n ∑ j = 0 m f [ c ] [ i ] [ j ] $\sum_{i=0}^n\sum_{j=0}^mf[c][i][j]$ （不一定要占领所有的格子）。

Source

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
inline int read() {
    int res = 0; bool bo = 0; char c;
    while (((c = getchar()) < '0' || c > '9') && c != '-');
    if (c == '-') bo = 1; else res = c - 48;
    while ((c = getchar()) >= '0' && c <= '9')
        res = (res << 3) + (res << 1) + (c - 48);
    return bo ? ~res + 1 : res;
}
const int R = 2005, N = 33, E = 13, MX = 1e9 + 9;
int n, m, c, orz, h, cnt[E], C[R][R], g[E][N][N], f[E][N][N];
int main() {
    int i, j, k, h, l; n = read(); m = read(); c = read();
    for (i = 1; i <= c; i++) orz += (cnt[i] = read());
    for (i = 0; i <= 2000; i++) C[i][0] = 1;
    for (i = 1; i <= 2000; i++) for (j = 1; j <= i; j++)
        C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]) % MX;
    for (k = 1; k <= c; k++) for (i = 1; i <= n; i++)
    for (j = 1; j <= m; j++) {
        g[k][i][j] = C[i * j][cnt[k]];
        for (h = 1; h <= i; h++) for (l = 1; l <= j; l++) if (h != i || l != j)
            g[k][i][j] = (g[k][i][j] - 1ll * g[k][h][l]
            * C[i][h] % MX * C[j][l] % MX + MX) % MX;
    }
    f[0][n][m] = 1; for (k = 1; k <= c; k++) for (i = 0; i <= n; i++)
    for (j = 0; j <= m; j++) for (h = i + 1; h <= n; h++)
    for (l = j + 1; l <= m; l++)
        f[k][i][j] = (f[k][i][j] + 1ll * f[k - 1][h][l]
        * g[k][h - i][l - j] % MX * C[h][i] % MX * C[l][j] % MX) % MX;
    int ans = 0; for (i = 0; i <= n; i++) for (j = 0; j <= m; j++)
        ans = (ans + f[c][i][j]) % MX; cout << ans << endl;
    return 0;
}

xyz32768

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
[题解]CLYZ2018省选训(bao)练(zha)模拟赛 Day 8

竟然卡SPFA，出题人好毒啊。题目T1：BZOJ 2330 / SCOI 2011 糖果 (candy)（差分约束） T2：BZOJ 2743 / HEOI 2012 采花 (flower)（离线+树状数组） T3：BZOJ 3294 / CQOI 2011 放棋子 (chess)（DP+组合数学+容斥）T1分析差分约束裸题。对于一个条件： 1AB1AB1\>...
复制链接

扫一扫