Product 1 Modulo N
题意:
给一个n,(2<=n<=1e5),找到排列[1,2,…,n-1]的最长子序列,使得它元素的乘积模n意义下为1.
输出子序列长度以及子序列。
思路:
考虑1~n-1中那些元素可以选,如果x与n不互质,即gcd(x, n) == gcd(x % n, n) != 1.所以x % n != 1.因此子序列元素的乘积得与n互质,所以考虑选择1~n-1中与n互质的元素,gcd(x, n) == gcd(x % n, n) == 1就符合题意,如果最后x % n不为1,那只要把x%n这个元素去除即可。证明:
设x % n == p, 则除了p之外的乘积为ans, ans * p % n == p, <=> ans % n == 1.
代码:
/*************************************************************************
> File Name: a.cpp
> Author: Beans
> Mail: 3112748286@qq.com
> Created Time: 2023/5/25 11:07:52
************************************************************************/
#include <iostream>
#include <algorithm>
#define int long long
#define endl '\n'
using namespace std;
const int maxn = 3e5 + 7;
int a[maxn];
int gcd(int a, int b){
return b ? gcd(b, a % b) : a;
}
signed main(){
int n;
cin >> n;
int ans = 1;
for(int i = 1; i <= n - 1; i ++ )
if(gcd(i, n) == 1){
a[i] = 1;
ans = (ans * i) % n;
}
if(ans != 1) a[ans] = 0;
ans = 0;
for(int i = 1; i <= n - 1; i ++ )
if(a[i]) ans ++ ;
cout << ans << endl;
for(int i = 1; i <= n - 1; i ++ )
if(a[i]) cout << i << ' ';
cout << endl;
}