【数学分析】Codeforces Round #426(Div.2)C[The Meaningless Game]题解

题目概述

有一个游戏，过程如下：
1.有两个数 $a,b$ ，刚开始 $a=b=1$ 。
2.随便取一个正整数 $k$ ，使 (1)$a=a\times k,b=b\times k^2$ 或 (2)$a=a\times k^2,b=b\times k$ 。
问最后 $a,b$ 是否能够变成给出的 $A,B$ 。

解题报告

其实我们只需要考虑 $k$ 是质因子就行了，其他情况由质因子组成即可。
那么考虑 $A,B$ 的唯一分解式，设第 $i$ 个素数为 $p_i$ ， $A_i$ 表示 $A$ 中 $p_i$ 的指数， $B_i$ 表示 $B$ 中 $p_i$ 的指数。
对于 $\forall i$ ，假设 $x$ 表示(1)操作的次数， $y$ 表示(2)操作的次数，那么可以得出一个方程组：
$x+2y=A_i$
$2x+y=B_i$
即：
$x=\frac{2B_i-A_i}{3}$
$y=\frac{2A_i-B_i}{3}$
所以有解的充要条件是 $A_i\le 2B_i,B_i\le 2A_i,3|(2B_i-A_i),3|(2A_i-B_i)$ 。
转化一下变为： $A|B^2,B|A^2$,$A^2\over B$ 的立方根是整数, $B^2\over A$ 的立方根是整数。
那么就很好判断了。

示例程序

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long LL;

int te;LL A,B;

bool check(LL A,LL B)
{
    if (A*A%B) return false;
    LL L=1,R=1000000,num=A*A/B;
    while (L<=R)
    {
        LL mid=L+(R-L>>1);
        if (mid*mid*mid>=num) R=mid-1; else L=mid+1;
    }
    if (L*L*L!=num) return false;
    if (B*B%A) return false;
    L=1;R=1000000;num=B*B/A;
    while (L<=R)
    {
        LL mid=L+(R-L>>1);
        if (mid*mid*mid>=num) R=mid-1; else L=mid+1;
    }
    if (L*L*L!=num) return false;
    return true;
}
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    scanf("%d",&te);
    while (te--)
    {
        scanf("%lld%lld",&A,&B);
        if (check(A,B)) printf("Yes\n"); else printf("No\n");
    }
    return 0;
}