【JZOJ5773】简单数学题【数论,数学】

题目大意:

题目链接:https://jzoj.net/senior/#main/show/5773
题目图片:
https://www.z4a.net/images/2018/09/22/Screenshot.md.png
https://www.z4a.net/images/2018/09/22/Screenshot-1.png
给出NN,求N12TNT\frac{N-\frac{1}{2}T}{N-T}值为正整数时TT的所有正整数解(0<T<N0<T<N)。


思路:

40分做法:

N106N\leq 10^6
直接暴力枚举TT,输出符合要求的即可。

100分做法:

首先,我们要求的是N12TNT\frac{N-\frac{1}{2}T}{N-T}
x=NTx=N-T,那么就有
N12Tx\frac{N-\frac{1}{2}T}{x}
我们知道,T=NN+T=N(NT)T=N-N+T=N-(N-T),所以
N12[N(NT)]x\frac{N-\frac{1}{2}[N-(N-T)]}{x}
用乘法分配律脱括号
N12N+12(NT)x\frac{N-\frac{1}{2}N+\frac{1}{2}(N-T)}{x}
拆成两半
N12Nx+12(NT)x\frac{N-\frac{1}{2}N}{x}+\frac{\frac {1}{2}(N-T)}{x}
左边可以简化
12Nx+12(NT)x\frac{\frac{1}{2}N}{x}+\frac{\frac {1}{2}(N-T)}{x}
一开始设了x=NTx=N-T,右边有一项就是NTN-T,所以我们可以把它简化成
12Nx+12xx\frac{\frac{1}{2}N}{x}+\frac{\frac {1}{2}x}{x}
右边约分得
12Nx+12\frac{\frac{1}{2}N}{x}+\frac {1}{2}
我们知道,这个式子的值必须是个整数。我们设答案为整数kk
12Nx+12=k\frac{\frac{1}{2}N}{x}+\frac {1}{2}=k
等号两边同时乘2得
Nx+1=2k\frac{N}{x}+1=2k
那么由于2k2k和1都是整数,所以就有Nx\frac{N}{x}也是整数
Nx\frac{N}{x}为整数时,xx就一定是NN的因数,所以NTN-T就一定是NN的因数!
那么我们就枚举NN的因数d[i]d[i],当d[i]+12\frac{d[i]+1}{2}为正整数(即d[i]d[i]是奇数)时,Nd[i]N-d[i]就是一个合法的解!


代码:

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

ll n,a[N],d[N],sum,ans;

int main()
{
	cin>>n;
	for (ll i=1;i<=sqrt(n);i++)  //求出n的所有因数
	 if ((n%i)==0)
	 {
	 	if (i*i==n) d[++sum]=i;
	 	else
	 	{
	 		d[++sum]=i;
	 		d[++sum]=n/i;
	 	}
	 }
	for (ll i=1;i<=sum;i++)
	{
		if (((n/d[i])%2)&&n-d[i])  //符合要求
		 a[++ans]=n-d[i];
	}
	sort(a+1,a+1+ans);
	cout<<ans;
	for (int i=1;i<=ans;i++)
	 cout<<' '<<a[i];
	return 0;
}
原文地址:https://www.cnblogs.com/hello-tomorrow/p/11998572.html