[ABC268G] Random Student ID Solution

更好的阅读体验戳此进入

题面

给定 $n$$n$ 名学生的名字 $S_i$$S_i$，保证不重复。现在随机一个 $a\to z$$a \rightarrow z$ 的排列作为字典序的偏序关系，求每个同学排名的期望，对 $998244353$$998244353$ 取模。

Solution

这一题可以说确实很高妙，确实没想到还可以这么处理。

第一眼以为是一些类似 期望DP 的东西，然后想了半天也没想到什么复杂度正确的做法，实际上这道题是从期望的意义去找性质。

首先我们考虑对于一个串 $S$$S$ 的期望排名，若存在串 $S$$S$ 的前缀串 $T$$T$，那么显然串 $T$$T$ 无论字典序的偏序关系如何都会在 $S$$S$ 之前，所以一定会贡献 $1$$1$。

而若存在 $S$$S$ 是 $T$$T$ 的前缀，那么显然 $T$$T$ 在 $S$$S$ 之后，则贡献一定为 $0$$0$。

对于其它的情况，显然会由 $\mathtt{\text{LCP}}$$\texttt{LCP}$ 之后的第一个不同的字符决定，而对于这种显然在所有字典序关系中两者的偏序关系的概率是相同的，所以此时的贡献为 ${\displaystyle \frac{1}{2}}$$\dfrac{1}{2}$。那么对于字符串 $S$$S$，若其前缀串有 $x$$x$ 个，满足有 $S$$S$ 前缀的串有 $y$$y$ 个，那么答案即为：

然后对于求前缀，显然可以通过 Trie树 解决，同时也可以考虑一些更简单的方法：

开一个 unordered_set，往里面丢每个串的哈希值，然后对于每个串的每个前缀的哈希在里面查一下然后记录即可，最终复杂度 $O(\sum |S_i|)$$O(\sum \lvert S_i \rvert)$。

Code

xxxxxxxxxx
79
1
#define _USE_MATH_DEFINES
2
#include <bits/stdc++.h>
3

4
#define PI M_PI
5
#define E M_E
6
#define npt nullptr
7
#define SON i->to
8
#define OPNEW void* operator new(size_t)
9
#define ROPNEW void* Edge::operator new(size_t){static Edge* P = ed; return P++;}
10

11
using namespace std;
12

13
mt19937 rnd(random_device{}());
14
int rndd(int l, int r){return rnd() % (r - l + 1) + l;}
15
bool rnddd(int x){return rndd(1, 100) <= x;}
16

17
typedef unsigned int uint;
18
typedef unsigned long long unll;
19
typedef long long ll;
20
typedef long double ld;
21

22
#define MOD (998244353ll)
23
#define BASE (31ll)
24

25
template < typename T = int >
26
inline T read(void);
27

28
int N;
29
int cntx[510000], cnty[510000];
30
ll inv_2;
31
string S[510000];
32
unordered_map < unll, int > idx;
33

34
ll qpow(ll a, ll b){
35
    ll ret(1), mul(a);
36
    while(b){
37
        if(b & 1)ret = ret * mul % MOD;
38
        b >>= 1;
39
        mul = mul * mul % MOD;
40
    }return ret;
41
}
42

43
int main(){
44
    inv_2 = qpow(2, MOD - 2);
45
    N = read();
46
    for(int i = 1; i <= N; ++i){
47
        cin >> S[i];
48
        unll cur(0);
49
        for(auto c : S[i])(cur *= BASE) += c;
50
        idx.insert({cur, i});
51
    }
52
    for(int i = 1; i <= N; ++i){
53
        unll cur(0);
54
        for(auto c : S[i]){
55
            (cur *= BASE) += c;
56
            if(idx.find(cur) != idx.end())++cntx[i], ++cnty[idx[cur]];
57
        }cntx[i]--, cnty[i]--;
58
    }
59
    for(int i = 1; i <= N; ++i)
60
        printf("%lld\n", (cntx[i] + (N - cntx[i] - cnty[i] - 1) * inv_2 % MOD + 1) % MOD);
61
    fprintf(stderr, "Time: %.6lf\n", (double)clock() / CLOCKS_PER_SEC);
62
    return 0;
63
}
64

65
template < typename T >
66
inline T read(void){
67
    T ret(0);
68
    int flag(1);
69
    char c = getchar();
70
    while(c != '-' && !isdigit(c))c = getchar();
71
    if(c == '-')flag = -1, c = getchar();
72
    while(isdigit(c)){
73
        ret *= 10;
74
        ret += int(c - '0');
75
        c = getchar();
76
    }
77
    ret *= flag;
78
    return ret;
79
}

UPD

update-2023_01_18 初稿