LightOJ 1341

最新推荐文章于 2021-08-22 22:45:49 发布

敬宅神社

最新推荐文章于 2021-08-22 22:45:49 发布

阅读量344

点赞数

CC 4.0 BY-SA版权

分类专栏：数论

本文链接：https://blog.csdn.net/qq_20743237/article/details/45848575

数论专栏收录该内容

16 篇文章

订阅专栏

题意：

一个数可以表示 N 对数相乘（x1，x2）

输入另一个数b 求有多少对x1 x2 俊大于等于 b

求标准分解式的指数因子数为 II（ai + 1）（表示不出数学符号.............）

代码：

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
using namespace std;
typedef long long LL;
const int maxn = 1e6 +131;
int vis[maxn];
LL Primes[maxn];
int Count,ans;
LL a,b;
void INT() { //质数表  
	Count = 0;
	for(int i = 2; i <= maxn; ++i) 
	{
		if(!vis[i]) {
			Primes[Count++] = i;
			for(int j = i; j <= maxn; j+= i) 
			vis[j] = 1;
		}
	}
}

void DFS(LL dep, LL x) { //构造解答树，求出小于b的有多少个a 的因数 
	ans--;
	for(int i = dep; i < Count; ++i) 
	{
		if(x * Primes[i] < b) 
		{
			if( a  % (x* Primes[i]) == 0) 
			DFS(i,x*Primes[i]);
		}
		else return ;
	}
}

int main() {
	INT();
	int T;
	scanf("%d",&T);
	for(int kase = 1; kase <= T; ++kase) {
		scanf("%lld%lld",&a,&b) ;
		if(b*b >= a) {
			printf("Case %d: 0\n",kase);
			continue;
		}
		ans = 1;
		LL tmp = a;
		for(int i = 0; i < Count && Primes[i] * Primes[i] <= tmp; ++i ) {
			if(tmp % Primes[i] == 0) {
				LL tot = 0;
				while(tmp % Primes[i] == 0) {
					tot++;
					tmp /= Primes[i];
				}
				ans *= (tot+1);
			} 
		}
		
		if(tmp > 1) ans <<= 1; //如果还有剩余的质数，则因子数还要 * 2 
		ans >>= 1;//求有多少对
		DFS(0,1); if(b == 1) ans ++;
		printf("Case %d: %lld\n",kase, ans);
	}
	return 0;
}