2020 Multi-University Training Contest 5

Tetrahedron

$\frac{1}{h^2}=\frac{1}{a^2}+\frac{1}{b^2}+\frac{1}{c^2}$

$E(x+y)=E(x)+E(y)$

// 百思不得其解 赛中3秒过了 赛后T cin改成scanf后4秒过了 hdu这是什么机制？？？
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod = 998244353;
const int N = 6e6 + 10;
int n;

ll qpow(ll a, ll b) {
    a %= mod;
    ll res = 1;
    while (b) {
        if (b & 1) res = res * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return res;
}

ll f[N], sum[N];

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);

    for (int i = 1, lim = 6e6; i <= lim; i++) {
        f[i] = qpow(1ll * i * i, mod - 2);
        sum[i] = (sum[i - 1] + f[i]) % mod;
    }

    int T;
    cin >> T;
    while (T--) {
        cin >> n;
        ll a = (qpow(n, mod - 2) * 3) % mod;
        cout << ((a * sum[n]) % mod) << endl;
    }
    return 0;
}

---

Funny String

Boring Game

官方题解

#include <bits/stdc++.h>
using namespace std;
const int N = 1e6 + 10;
vector<int> v[N];
int a[N];
int n, m, k;

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0), cout.tie(0);
    
    int T;
    cin >> T;
    for (int cs = 1; cs <= T; cs++) {
  
        cin >> n >> k;
        int lim = 2 * n * (1 << k);
        m = (1 << k);
        for (int i = 1; i <= m; i++) {
            v[i].clear();
        }
        for (int i = 1; i <= lim; i++) {
            cin >> a[i];
            v[m].push_back(a[i]);
        }
        reverse(v[m].begin(), v[m].end());
        for (int t = 0; t < k; t++) {
            for (int i = m - (1 << t) + 1, j = i - 1; i <= m; i++, j--) {
                int mid = v[i].size() / 2;
                for (int l = 0; l < mid; l++) {
                    v[j].push_back(v[i].back());
                    v[i].pop_back();
                }
            }
        }

        for (int i = 2 * n - 1; i >= 0; i--) {
            for (int j = 1; j <= m; j++) {
                cout << v[j][i] << (j == m && i == 0 ? "\n" : " ");
            }
        }
    }
    return 0;
}

Expression
Array Repairing
Alice and Bob
Tree
Set2

---

Paperfolding

经画图 _{和撕了N张纸} 后发现，左右对折、上下对折之间的顺序没有任何关系，

假设左右对折了 $x$ 次，上下对折了 $y$ 次，最后横竖裁开，纸片个数 $=(2^x+1)(2^y+1)$

$$E(x)=\sum x\times p(x)=\sum _{x=0}^n\frac{C_n^x(2^x+1)(2^{n-x}+1)}{2^n}$$

其中分子部分
$$\begin{gathered} \sum _{x=0}^n{C}_n^x(2^x+1)(2^{n-x}+1) \\ =\sum _{x=0}^n{C}_n^x(2^n+2^x+2^{n-x}+1) \\ =2^n\sum _{x=0}^n{C}_n^x+(\sum _{x=0}^n{C}_n^x{2}^x+\sum _{x=0}^n{C}_n^x{2}^{n-x})+\sum _{x=0}^n{C}_n^x \end{gathered}$$

观察发现
$$\sum _{x=0}^n{C}_n^x{2}^x=\sum _{x=0}^n{C}_n^x{2}^x\times 1^{n-x}=(2+1{)}^n=3^n$$

对，就是化成二项式，就是下面这个↓

$$(a+b{)}^n=\sum _{k=0}^n{C}_n^k{a}^k{b}^{n-k}$$
当 $a=1$ 且 $b=1$ 时，

$$(1+1{)}^n=\sum _{k=0}^n{C}_n^k{1}^k{1}^{n-k}=\sum _{k=0}^n{C}_n^k=2^n$$

然后回到本题，同理可得
$$\sum _{x=0}^n{C}_n^x{2}^{n-x}=\sum _{x=0}^n{C}_n^x{2}^{n-x}\times 1^x=(1+2{)}^n=3^n$$

因此化简原来的公式可以得出

$$E(x)=2^n+1+\frac{2\times 3^n}{2^n}$$

记个和这题没什么关系的笔记：

裁纸数个数的时候发现，如果只有左右对折 or 上下对折时，记 $a_x$ 为对折了 $x$ 次后纸片的个数

有如下数列：
$a_1=6$、$a_2=10$、$a_3=18$、$a_4=34$ …

$a_2-a_1=4$、$a_3-a_2=8$、$a_4-a_3=16$…

求像这种公差为等比数列的通项公式，一般分为如下步骤：

以上面的数列为例，记
$$a_{n+1}=a_n+2^{n+1}$$

然后让等式左右两边都减去同一个数 $x$，

$$a_{n+1}-x=a_n+2^{n+1}-x$$

这个数是自己构造的，构造出来的目的是使等式左右两边相似，这里令 $x=2^{n+2}$

$$a_{n+1}-2^{n+2}=a_n+2^{n+1}-2^{n+2}=a_n+2^{n+1}-2\times 2^{n+1}=a_n-2^{n+1}$$

发现 $a_n-2^{n+1}$ 是个公比 $q=1$ 的等比数列，

因此可以直接往前推

$$a_n-2^{n+1}=...=a_1-2^2=6-4=2$$

$$a_n=2^{n+1}+2$$

总结：若存在 $a_{n+1}=a_n+b^k$，则构造通项公式 $a_n-\frac{b^{k+1}}{b-1}$

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n;

namespace MOD { // 取模运算板子
    const ll mod = 998244353;

    ll qpow(ll a, ll b) {
        a %= mod;
        ll res = 1;
        while (b) {
            if (b & 1) res = res * a % mod;
            a = a * a % mod;
            b >>= 1;
        }
        return res;
    }

    ll Mul(ll a, ll b) {
        return (a * b) % mod;
    }

    ll Sub(ll a, ll b) {
        return (a - b + mod) % mod;
    }

    ll Add(ll a, ll b) {
        return (a + b) % mod;
    }
}
using namespace MOD;

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);

    int T;
    cin >> T;
    while (T--) {
        cin >> n;
        if (n == 0) cout << 4 << endl;
        else if (n == 1) cout << 6 << endl;
        else {
            ll _2n = qpow(2, n);// 2^n
            ll _3n = qpow(3, n); // 3^n
            ll Inv_2n = qpow(_2n, mod - 2); // 1/2^n

            ll res = Add(1, Add(_2n, Mul(Mul(2, _3n), Inv_2n)));
            cout << res << endl;
        }
    }
    return 0;
}

Function
Exam

---

Set1

官方题解有问题，就不贴了…

$cnt[i]$ 表示第 $i$ 个元素被留到最后的方案数

此时，被删去的元素有前面 $i-1$ 个 和 后面 $n-i$ 个

前面 $i-1$ 个元素里有操作1删去的，也有操作2随机删去的，但后面$n-i$个元素一定是被操作2随机删去的，

所以当且仅当 $n-i\le i-1$ 的时候，$i$ 才有可能被留下

题目规定操作1和操作2轮流操作，所以在前面 $i-1$ 个元素里挑选 $n-i$ 个，和后面 $n-i$ 个元素匹配，

且由于操作1每次都是删除整个集合里最小的元素，所以在前面$i-1$个元素里选出来的$n-i$个元素，一定是有序的，但后面的$n-i$个元素是无序的，所以这一部分的方案数有

$$C_{i-1}^{n-i}\times (n-i)!$$

再考虑前面剩下的 $i-1-(n-i)=2i-n-1$ 个元素两两配对，

我们知道如果n个元素里两两配对，会有
$$C_n^2\times C_{n-2}^2\times C_{n-4}^2...\times C_2^2$$$$=\frac{n(n-1)}{2}\times \frac{(n-2)(n-3)}{2}\times ...\times \frac{2\times 1}{2}=\frac{n!}{2^{\frac{n}{2}}}$$

其中两两分配的元素相同，但前后顺序不同（12|34、 34|12 这样的）有 $\frac{n}{2}!$ 种

但是在前面剩下的 $2i-n-1$ 个元素里，

分组后每一组里的两个数显然是小数为操作1删除的，大数为操作2删除的，是固定的顺序，

且所有组先后发生的顺序也是固定的，如果有分组如 12|34，显然3、4的删除不可能能发生在1、2的删除之前

所以这一部分的方案数有
$$\frac{(2i-n-1)!}{2^{\frac{2i-n-1}{2}}}\times \frac{1}{(\frac{2i-n-1}{2})!}$$

所以第 $i$ 个元素被留到最后的方案数为
$$cnt[i]=C_{i-1}^{n-i}\times (n-i)!\times \frac{(2i-n-1)!}{2^{\frac{2i-n-1}{2}}}\times \frac{1}{(\frac{2i-n-1}{2})!}$$

总方案数 $sum=\sum cnt[i]$

每个数被留下的概率 $p[i]=\frac{cnt[i]}{sum}$

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod = 998244353;
const int N = 5e6 + 10;
const int M = 5e6;
int n;

ll Mul(ll a, ll b) { return (a * b) % mod; }

ll Add(ll a, ll b) { return (a + b) % mod; }

ll Sub(ll a, ll b) { return (a - b + mod) % mod; }

ll qpow(ll a, ll b) {
    a %= mod;
    ll res = 1;
    while (b) {
        if (b & 1) res = Mul(res, a);
        b >>= 1;
        a = Mul(a, a);
    }
    return res;
}

ll fac[N], inv[N];
ll inv_2[N];// 1/(2^i)

void init(int n) {
    fac[0] = inv[0] = inv_2[0] = 1;
    int inv2 = qpow(2, mod - 2);
    for (int i = 1; i <= n; ++i) {
        fac[i] = Mul(fac[i - 1], i);
        inv_2[i] = Mul(inv_2[i - 1], inv2);
        inv[i] = Mul(inv[i - 1], qpow(i, mod - 2));
    }
}

ll C(int n, int m) {
    if (n < m || m < 0)return 0;
    return Mul(Mul(fac[n], inv[m]), inv[n - m]);
}

ll cnt[N];

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0), cout.tie(0);

    init(M);
    int T;
    cin >> T;
    for (int cs = 1; cs <= T; cs++) {

        cin >> n;
        ll sum = 0;
        for (int i = 1; i <= n; i++) {
            int tmp = i * 2 - n - 1;
            if (tmp < 0) cnt[i] = 0;
            else
                cnt[i] = Mul(Mul(C(i - 1, n - i), fac[n - i]),
                             Mul(Mul(fac[tmp], inv_2[tmp / 2]), inv[tmp / 2]));
            sum = Add(sum, cnt[i]);
        }
        
        ll inv_sum = qpow(sum, mod - 2);
        for (int i = 1; i <= n; i++) {
            cout << Mul(cnt[i], inv_sum) << (i == n ? "\n" : " ");
        }
    }
    return 0;
}

An Easy Matrix Problem