【2015广东工业大学新生赛B】【数论 素因子拆分】小刀带你们玩游戏 lcm（A,B）==C下B的个数

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<ctype.h>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre(){freopen("c://test//input.in","r",stdin);freopen("c://test//output.out","w",stdout);}
#define MS(x,y) memset(x,y,sizeof(x))
#define MC(x,y) memcpy(x,y,sizeof(x))
#define MP(x,y) make_pair(x,y)
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1,class T2>inline void gmax(T1 &a,T2 b){if(b>a)a=b;}
template <class T1,class T2>inline void gmin(T1 &a,T2 b){if(b<a)a=b;}
const int N=0,M=0,Z=1e9+7,ms63=1061109567;
int casenum,casei;
int A,C;
const int TOP=32000;
bool e[TOP+10];
int p[TOP/5],pnum=0;
void prime()
{
	for(int i=2;i<=TOP;++i)if(e[i]==0)
	{
		p[++pnum]=i;
		for(int j=i;j<=TOP;j+=i)e[j]=1;
	}
}
pair<int,int>a[12];
int main()
{
	prime();
	scanf("%d",&casenum);
	for(casei=1;casei<=casenum;++casei)
	{
		scanf("%d%d",&A,&C);
		int G=C/A;
		int n=0;
		for(int i=1;p[i]*p[i]<=A;++i)if(A%p[i]==0)
		{
			a[++n]=MP(p[i],0);
			while(A%p[i]==0)
			{
				A/=p[i];
				++a[n].second;
			}
		}
		if(A>1)a[++n]=MP(A,1);
		int ans=1;
		for(int i=1;i<=n;++i)
		{
			if(G%a[i].first)//如果不含有这个素因子，那么B中a[i].first的个数可以是[0,a[i].second]
			{
				ans*=a[i].second+1;
			}//否则B中a[i].first的个数是确定而且唯一的
		}
		printf("%d\n",ans);
	}
	return 0;
}
/*
【trick&&吐槽】
1，这个题意有漏洞——
"规则是这样的：两个人以相同的速度往前跳跃，但是每次跳跃距离不同"
于是判定A!=B？判定就错啦>_< 我们应该自然忽略这句话！理解为A与B可能不同。

2，预处理素数加速，我的速度果然是最快的，棒棒棒！

【题意】
给你2个数A和C，都是1e9范围的正整数且保证C是A的倍数
问你，存在有多少个B，使得lcm（A，B）==C。

【类型】
数论

【分析】
lcm（A，B）==C，我们设G=C/A。

事实上，C可以写成（p1^w1）*（p2^w2）*（...）*（pn^wn）的形式，px代表互不相同的素数。
那么显然，对于pi^wi，要不就是A中恰好含有wi个pi，要不就是B中含有wi个pi。

于是，我们对A做因子拆分。
如果A中恰好含有wi个pi，也就是G中不再含有pi，那么B中可以含有[0,wi]个pi，于是ans要*=(wi+1)
如果A中不是恰好含有wi个pi，个数要少，也就是G中依然含有pi，那么B中一定含有wi个pi。
也就是说，G中只需要看A中有的素因子它是否有即可，然后这道题就做完喽。

【时间复杂度&&优化】
O（Tsqrt(A)）

*/