2021-08-16 SSL 模拟赛 T3 【约数】【思维】【数学】

2021-08-16 SSL 模拟赛 T3

题目大意：

求出$1-n-1$ 中所有满足$\frac{n-0.5x}{n-x}=k,k\in N\ast$ 的 $x$。

思路：

我们来转化一下这个公式，让它变得好看一点：
$$\frac{n-0.5x}{n-x}=\frac{2n-x}{2n-2x}=k$$
移项得到：
$$2n-x=n+(n-x)=k(2n-2x)=2k(n-x)$$
然后左边只保留 $n$ ：$$n=(2k-1)(n-x)$$
所以就有$(n-x)$是 $n$ 的约数。
然后我们可以在$O(\sqrt{n})$ 的时间复杂度内求出 $n$ 的所有约数并且求出$x$
注意$n/(n-x)=2k-1$ $2k-1$ 是奇数，这是另一个约束条件。

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cmath>
#define r register
#define rep(i,x,y) for(r ll i=x;i<=y;++i)
#define per(i,x,y) for(r ll i=x;i>=y;--i)
using namespace std;
typedef long long ll;
const ll V=1e7+10;
ll n,top,a[V],now;
ll ans[V],cnt;
void cut()
{
	rep(i,1,now)
	 if(n%i==0)
	 {
	 	a[++top]=i;
	 	ll x=n/i;
	 	if(i!=x) a[++top]=n/i;
	 }
}
void put(ll x)
{
	if(x<0) putchar('-'),x*=-1;
	if(x>9) put(x/10);
	putchar(x%10+48);
}
int main()
{
	scanf("%lld",&n);
	now=sqrt(n);
	cut();
	sort(a+1,a+top+1,greater<ll>());
	rep(i,1,top)
	 if(n-a[i]>0)
	  if((n/a[i])&1)
	   ans[++cnt]=n-a[i];
	put(cnt);
	putchar(' ');
	rep(i,1,cnt)
	 put(ans[i]),putchar(' ');
	return 0;
}