题库练习 「USACO 2020 US Open Platinum」Exercise
← 上一题 下一题 →

A5903 | 「USACO 2020 US Open Platinum」Exercise

时间限制1s
内存限制512MB
通过 / 提交0/0

题目描述

**题目译自 [USACO 2020 US Open Contest, Platinum](http://usaco.org/index.php?page=open20results) Problem 2. [Exercise](http://usaco.org/index.php?page=viewproblem2&cpid=1033)**

Farmer John(又)想到了一个新的奶牛们的早操计划!

一如既往地,Farmer John 的 $N$ 头奶牛站成一排。左数第 $i$($1 \le i \le N$)头奶牛的编号为 $i$。FJ 让奶牛不断重复执行下列操作,直到奶牛们又回到一开始的站位为止:

- 给定一个 $1 \sim N$ 的排列 $A$,使得原位置为左数第 $i$ 个的奶牛的新位置为左数第 $A_i$ 个。

例如,如果 $A = (1, 2, 3, 4, 5)$,则奶牛们执行一次操作后就立刻回到了原站位。如果 $A = (2, 3, 1, 5, 4)$,则奶牛们需要执行 $6$ 次操作才能回到原站位。每次执行完的站位分别是:

- 第 $0$ 次后:$(1, 2, 3, 4, 5)$。
- 第 $1$ 次后:$(3, 1, 2, 5, 4)$。
- 第 $2$ 次后:$(2, 3, 1, 4, 5)$。
- 第 $3$ 次后:$(1, 2, 3, 5, 4)$。
- 第 $4$ 次后:$(3, 1, 2, 4, 5)$。
- 第 $5$ 次后:$(2, 3, 1, 5, 4)$。
- 第 $6$ 次后:$(1, 2, 3, 4, 5)$。

请你计算对于所有 $N!$ 种 $1 \sim N$ 的排列 $A$ 所需次数的乘积。

因为这个数可能非常大,所以请你输出答案对 $M$ 取模的结果,保证 $M$ 为质数。

下列来自 [KACTL](https://github.com/kth-competitive-programming/kactl/blob/master/content/various/FastMod.h) 的代码可能会对使用 C++ 语言的用户产生一定的帮助。此方法名为 [Barrett reduction](https://en.wikipedia.org/wiki/Barrett_reduction),它允许你使用更快的速度多次计算 $a \bmod b$,条件是 $b > 1$ 且为常数,但无法在程序编译期被得知。

```cpp
#include <bits/stdc++.h>
using namespace std;

typedef unsigned long long ull;
typedef __uint128_t L;
struct FastMod {
ull b, m;
FastMod(ull b) : b(b), m(ull((L(1) << 64) / b)) {}
ull reduce(ull a) {
ull q = (ull)((L(m) * a) >> 64);
ull r = a - q * b; // can be proven that 0 <= r < 2*b
return r >= b ? r - b : r;
}
};
FastMod F(2);

int main() {
int M = 1000000007; F = FastMod(M);
ull x = 10ULL*M+3;
cout << x << " " << F.reduce(x) << "\n"; // 10000000073 3
}
```

输入格式

从标准输入中读入数据。

仅一行两个正整数 $N, M$。

输出格式

输出到标准输出中。

一行一个数表示答案。

输入输出样例

输入 #1
5 1000000007
输出 #1
369329541
C++ 编辑器
输入
输出