【模拟赛】【数论】2021.8.12.D

题目描述

给出n，求所有x，使得x小于n，且$x^2\%n=1$
将所有x从小到大输出

思路

$x^2\%n=1$
$(x^2-1)\%n=0$
$x^2-1=kn$
$(x-1)(x+1)=kn$
然后设
$(x+1)=k_1\ast n_1$
$(x-1)=k_2\ast n_2$
$k_1\ast k_2=k$
$n_1\ast n_2=n$
我们枚举$n_1$并记为$a$
那么
$n_2=n/a$
$(x+1)=k_1\ast a$
我们再枚举$k_1$，记作$b$
则$(x+1)=a\ast b$
即$x=a\ast b-1$
将此式代入$(x-1)=k_2\ast n_2$
得$(a\ast b-1)=k_2\ast (n/a)$
若$(a\ast b-2)$为$(n/a)$的倍数
则$x=a\ast b-1$为一组解
其余参考大爷博客

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#define ll long long
using namespace std;

ll Ans[5000250];
ll n, m, t;

void init(ll n1)
{
	for(ll k1 = 0; k1 <= n / n1; ++k1)
	{
		m = n1 * k1 - 2;//处理(x-1)和(x+1)
		if(m % (n / n1) == 0)Ans[++t] = n1 * k1 - 1;
		m = n1 * k1 + 2;
		if(m % (n / n1) == 0)Ans[++t] = n1 * k1 + 1;
	}
}

int main()
{
	scanf("%lld", &n);
	if(n == 1){printf("None\n"); return 0;}
	for(int i = 1; i <= (ll)sqrt(n); ++i)
		if(!(n % i))init(n / i);
	sort(Ans + 1, Ans + t + 1);
	t = 0;
	while(Ans[++t] <= 0);t--;
	while(Ans[++t] && Ans[t] <= n)if(Ans[t] != Ans[t - 1])printf("%lld\n", Ans[t]);
	return 0;
	
}