I love permutation

最新推荐文章于 2021-07-30 16:28:58 发布

原创最新推荐文章于 2021-07-30 16:28:58 发布 · 229 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

补题专栏收录该内容

42 篇文章

订阅专栏

博客探讨了在给定奇素数p和整数a下，生成序列A的逆序对对数模2的结果。通过数学证明，序列A是一个排列，并通过分析逆序对的奇偶性，利用扩展欧拉定理简化计算，最终确定逆序对数量的奇偶性。代码实现中涉及高精度快速幂运算，讨论了不同编程语言处理大数据的挑战。

题目链接

题意:给定一个奇素数 $p$ 和一个小于 $p$ 的整数 $a$ ,按如下规则生成一序列 $A$ : $A_i = a*i%p$ , $1$ $\leq$ $i$ $<$ $p$ ,问A序列的逆序对对数模 $2$ 的结果

解析:
①可以证明A是一个又 $1$ $\cdots$ $p - 1$ 构成的排列:不妨假定在取模后 $A_i=A_j$ ,且 $i > j$ 由定义有在取模之前 $A_i=k_1*p+c=a*i$ , $A_j=k_2*p+c=a*j$ ,则 $A_i-B_j=(k_1-k_2)*p=a*(i-j)$ ,移项有 $i - j =$ $(k_1-k_2)*p \over a$ $< p$ ,由于 $i - j$ 为整数且 $a, p$ 互质,所以 $\over a$ 为整数,由 $(k_1-k_2)*p \over a$ $< p$ 得 $\over a$ $\leq$ $1$ ,若 $\over a$ $= 0$ ,则有取模之前 $A_i=B_j$ ,与 $a * i$ $\neq$ $a * j$ 矛盾,若 $\over a$ $= 1$ ,则有 $a * i$ $\geq$ $a * p + c$ ,与 $1$ $\leq$ $i$ $<$ $p$ 矛盾,故不存在取模后 $A_i=A_j$ 的情况
②当 $i < j$ 时,如果 $A_i,A_j$ 构成逆序对,则 $(A_i-A_j) \over (i-j)$ $< 0$ ，反之 $(A_i-A_j) \over (i-j)$ $＞ 0$ ,所以判断一个序列中的逆序对对数的奇偶性只需要判断 $\prod_{0<i<j<p}$ $(A_i-A_j) \over (i-j)$ 的正负即可,当A是一个排列时, $\prod_{0<i<j<p}$ $(A_i-A_j) \over (i-j)$ $1)^x$ ,因为 $i,j和A_i,A_j$ 都是从一个由 $1$ $\cdots$ $p - 1$ 的排列中取数,所以 $\prod_{0<i<j<p}$ $A_i-A_j)$ 和 $\prod_{0<i<j<p}$ $(i - j)$ 的绝对值是相等的
③带入题目A的定义式化简2中式子, $\prod_{0<i<j<p}$ $(A_i-A_j) \over (i-j)$ $=$ $\prod_{0<i<j<p}$ $\over (i-j)$ $=$ $\prod_{0<i<j<p}$ $\over (i-j)$ $* a$ $=$ $a^{p*(p-1) \over 2}$ ,由扩展欧拉定理可化为 $2^{p-1 \over 2}$ (由于显示原因上述运算没加上取模),所以我们只需要判断 $a^{p*(p-1) \over 2}$ 是否为1即可,若为1,则说明逆序对数为偶数个

吐槽一波数据范围真的毒,C++好像只能__int128过?(也可能我的高精太慢了?),好奇java选手怎么过的,试着写了一次java大数然后T了

AcCode:

#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<vector>
#include<cstdio>

#define int long long
#define iter vector<int>::iterator

using namespace std;

const int base = 100000000;
const int N = 300;


int quick_pow(__int128 a, __int128 b, __int128 p)
{
	__int128 ans = 1;
	while (b)
	{
		if (b & 1)
			ans = ans * a % p;
		a = a * a % p;
		b >>= 1;
	}
	return ans;
}
inline int mul(int a, int b, int p) {
	int ans = 0;
	while (b) {
		if (b & 1) ans = ans + a % p;
		a = a + a % p;
		b >>= 1;
	}
	return ans;
}

signed main() {
	int t; scanf("%lld", &t);
	while (t--) {
		int a, p; scanf("%lld %lld", &a, &p);
		int pw = (p - 1) / 2;
		int ans = quick_pow(a, pw, p);
		if (ans == 1) printf("0\n");
		else printf("1\n");
	}
}