题意:
Please calculate the quantity of different functions
f
satisfying that
f(i)=bf(ai)
for each
i
from
0
to
n−1
.
(嗯...就是这样)
思路:
首先可以将 a 数组划分成若干个循环节,每个节是相互独立的,取其中一个节来看,令其长度为 l
显然,由定义,该节就决定了 bx = y, by = z, ... bxx = x, 的这个循环节有多长,
再看 b,显然我们应该将 b 的循环节 (设长度为 l0) 对应赋给上述循环节,至于长度呢,
1. l0 == l,显然可行,总共有 l 种排列
2. l0 == 1, 显然可行,注意到,此时第一个所赋的值决定了这 l 个式子的值,
即,若赋了 bx = x, 则接下来 l - 1 个式子必然也都是 bx = x, 因为第一个式子决定了 y = x, 继续下去, z = y, 以此类推,
故对每种赋值方法,有且只有 1 种排列,
3. 根据第二点就很容易想到了,l 的所有约数长度的循环节都是可以的,并且第一整段 l0 长度的式子决定了后面若干个整段 l0 长度的式子(和第一整段均一模一样),
故,一旦选定了b的某个循环节,且一旦确定了第一整段的排列,那么长度为 l 的整段的排列都唯一确定了,
而第一整段的排列有 l0 种,
答案就是 ∏i=1k∑j∣lij⋅calj ,其中 k 是置换 a 中循环节的个数, li 表示置换 a 中第 i 个循环节的长度, calj 表示置换 b 中长度为 j 的循环节的个数。
AC代码如下:
#include <cstdio>
#include <cmath>
#include <cstring>
#define mod 1000000007
#define maxn 100010
typedef long long LL;
int kas, a[maxn], b[maxn], n, m;
LL tot[maxn];
bool in[maxn];
void findcycle(int* a, int n) {
for (int i = 0; i < n; ++i) {
if (!in[i]) {
in[i] = true;
int ii = a[i], cnt = 1;
while (ii != i) {
in[ii] = true;
ii = a[ii];
++cnt;
}
++tot[cnt];
}
}
}
void work() {
for (int i = 0; i < n; ++i) scanf("%d", &a[i]);
for (int i = 0; i < m; ++i) scanf("%d", &b[i]);
memset(tot, 0, sizeof(tot));
memset(in, 0, sizeof(in));
findcycle(b, m);
memset(in, 0, sizeof(in));
LL ans = 1;
for (int i = 0; i < n; ++i) {
if (!in[i]) {
in[i] = true;
int ii = a[i], cnt = 1;
while (ii != i) {
// printf("%d ", ii);
in[ii] = true;
ii = a[ii];
++cnt;
}
// printf("\n");
// printf("%d %d\n", i, cnt);
LL mul = 0;
for (int j = 1; j <= cnt; ++j) {
if (cnt % j == 0) {
mul += (tot[j] * j) % mod;
mul %= mod;
}
}
ans *= mul;
ans %= mod;
}
}
printf("Case #%d: %lld\n", ++kas, ans % mod);
}
int main() {
while (scanf("%d%d\n", &n, &m) != EOF) work();
return 0;
}
怎么说呢, 虽然有点马后炮...
然而,这道题应该还是挺好做的
主要是当时太浮躁,
一来认为自己做不出来
二来被前面 02 WA了太多次很丧气
三来感觉当时做得人并不太多(也不知哪来的错觉...
四来就是其它一些琐碎的原因了...
其实主要吧...还就是不想做(叹气
结束之后定下来想了一想,很快就想得八九不离十了...
嗯...挺遗憾的......
接下来都要振作起来啊