第13届景驰-埃森哲杯广东工业大学ACM程序设计大赛 F 等式（求因子个数）

最新推荐文章于 2019-05-08 15:43:52 发布

原创最新推荐文章于 2019-05-08 15:43:52 发布 · 283 阅读

0 ·

CC 4.0 BY-SA版权

数学类专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一种高效的算法来解决给定正整数n时，求解1/x+1/y=1/n（x≤y）的正整数解的数量问题。通过数学转换将问题转化为求n²的因子个数，并给出了一段C++代码实现。

链接：https://www.nowcoder.com/acm/contest/90/F
来源：牛客网

题目描述

给定n，求1/x + 1/y = 1/n （x<=y）的解数。（x、y、n均为正整数）

输入描述:

在第一行输入一个正整数T。
接下来有T行，每行输入一个正整数n，请求出符合该方程要求的解数。
（1<=n<=1e9）

输出描述:

输出符合该方程要求的解数。

输入

输出

1
5
181

解：易得 x=n*y/(y-n) 设 y=y1+n; 则 x=n*n/y1 +n 因为 n,x,y1都为整数则y1是n*n的因子

那么问题变成了求n^2的因子个数，比如 n=4 ，n*n=16;

y1 y x

1 5 20

2 6 12

4 8 8

8 12 6

16 20 5

显然答案为ans=（（n^2）的因子数+1）/2;

设 p1^a1 *p2^a2 * ..*pn ^ an=n

那么 (a1+1)*(a2+1)*..*(an+1)为n的因子数

显然 n*n的因子数为 (2*a1+1)*(2*a2+1)*..*(2*an+1)

#include <bits/stdc++.h>
using namespace std;
int main()
{
    int T,n;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        int res=1;
        for(int i=2;i*i<=n;i++)
            if(n%i==0)
            {
                int cnt=0;
                while(n%i==0)
                {
                    cnt++;
                    n/=i;
                }
                res*=(cnt*2+1);
            }
        if(n>1) res*=3;
        printf("%d\n",(res+1)/2);
    }
    return 0;
}