PTA满足条件的数字有多少？郑州大学实验作业（超时问题的一点理解）

   这道题就是一个个枚举，a1到b1之间的数字，然后判断x与a0的最大公约数是否为a1，以及求出x与b0的最大公约数，再用x与b0的乘积除以他俩的的最大公约数看结果是否等于b1。（ps：两个数的最大公约数与最小公倍数的乘积等于这两个数的乘积）如果是就计数加一。

   但这道题b1的最大值可以给到20亿，如果一个个枚举a1到b1之间的值很容易超时，在这里我也没想出能解决超时问题的有效算法，也在csdn上搜了搜，看到一个前辈写的题解。但我没理解，后来自己想了一节课时间才懂，特地想写一篇博客来记录一下这个方法的巧妙。

   首先我们要找的x他必须满足“b1%x==0”这个条件吧。所以我们可以在枚举时先把b1的约数找出来，然后再在b1的约数中间判断有几个数满足上述两个条件。然后找b1的约数的时候我们同样不可以在a1到b1之间一个个枚举。

下面一段是我对为什么枚举1到根号b1之间的数，就可以找到b1所有约数的理解

  我们可以在1到根号b1之间枚举，如果x是b1的约数，那么b1/x也一定是b1的约数。如果x不是b1的约数，那么b1/x也无法除尽，那么b1/x也自然不是b1的约数。那么这样我们首先已经找齐了1到根号b1之间b1的约数，也找了“一部分”根号b1到b1之间的约数。那这“一部分”是不是包含了所有根号b1到b1之间的约数呢，答案当然是是的。我们可以这样想，b1除以在1到根号b1之间的b1的约数的结果的范围在根号b1到b1之间，结果也都可以整除b1。那么也就是说b1除以在根号b1到b1之间他的约数的结果的范围在1到根号b1之间，而且结果也都可以整除b1。那么我们只需要把 在1到根号b1之间b1的约数，就可以通过一一对应关系找到根号b1到b1的约数。并且找全。因为如果没找全根号b1到b1之间的约数的话，那么肯定b1除以它的结果在1到根号b1之间并且这个结果可以整除b1，但我们现在已经把1到根号b1之间能整除b1的数找全了。所以我们通过枚举1到根号b1，就可以找到能整除b1所有的数。

  这让我想到了我们在判断一个数是不是素数时，只需要枚举2到根号这个数之间的数，看有没有能整除这个数的数就可以了，它也可以缩短代码所需要的时间。

  所以我们就是通过寻找枚举1到根号b1之间的数，来寻找b1的所有的约数这个方法来解决超时问题。

下面把我的代码，贴上，其实还有一些要注意的小问题，我到代后面再阐述一下。

#include<stdio.h>
int zw(int z,int w){  //用辗转相除法得到两个数的最大公约数 
	int p=w;
	while(z%w!=0){
		p=z%w;
		z=w;
		w=p;
	} 
	return p;
} 
int main (){
	int t;
	scanf("%d",&t);
	for(int j=1;j<=t;j++){
		int a0,a1,b0,b1,wz=0;
		scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
		for(int i=1;i*i<=b1;i++){   //这时如果用i<=sqrt(b1)的话还会超时/
			if(b1%i==0){            //可能开方耗费的时间更长吧。 
				if(zw(i,a0)==a1&&zw(b1/b0,b1/i)==1) wz++;
				int w=b1/i;         //这里的zw(b1/b0,b1/i)我在代码后面再说说我的理解吧 
				if(i==w) continue;
				if(zw(w,a0)==a1&&zw(b1/b0,b1/w)==1) wz++;
			}
		}
		printf("%d\n",wz);
	}
	return 0;
}

  至于那里的zw(b1/b0,b1/i)==1是判断x与b0的乘积除以他俩的的最大公约数的结果是否等于b1的。

  经过了一点变形，也就是zw(a,b)==k等价于zw(a/k,b/k)==1,也就是说zw(i,bo)==i*b0/b1等价于zw(b1/b0,b1/i)==1。但我用第一个条件来判断一些结果是不对的，我猜i*b0当两个过大时应该会超出int类型的范围吧，然而开long long又可能会超时，所以只能用等价替换的了。

  同时注意i如果与b1/i相等的话，说明我们找到了两个相等的因数，只用判断一次就可以了。

  这个题解就结束了，希望能对你有所帮助吧。

  还有这是萌新第一次用博客记录自己对题目的看法，如果对你有帮助可以点赞一下哦。因为点赞可以让我看到有人看到这篇文章了哈哈，也算是对我的一点点鼓励。感谢！