[蓝桥杯 2023 国 B] 双子数

【问题描述】
  若一个正整数 x 可以被表示为 p^2 * q^2 ，其中 p、q 为质数且 p , q，则 x 是一个 “双子数”。请计算区间 [2333, 23333333333333] 内有多少个 “双子数”？

【答案提交】
  这是一道结果填空的题，你只需要算出结果后提交即可。本题的结果为一个整数，在提交答案时只填写这个整数，填写多余的内容将无法得分。

【思路】
由于23333333333333的平方根约为5e6，只需先求出5e6内的素数（欧拉筛和埃氏筛法都可），再根据题目要求计算乘积即可。
需注意剪枝操作，不然可能爆long long
结果为：947293

【代码】

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll maxn=1e7+10; 
ll cnt;
string s;
ll isprime[maxn],prime[maxn];
void f(){
	for (ll i=2;i<=maxn;i++){
		if (isprime[i]==0) prime[++cnt]=i;
		for (ll j=1;j<=cnt&&i*prime[j]<=maxn;j++){
			isprime[i*prime[j]]=1;
			if (i%prime[j]==0) break;
		}
	}
}
int main(){
  ios::sync_with_stdio(false);
  ll n,l2,t,l1,l3,l4,h=0,m,k,cnt1=0,cnt2=0,cnt3=0,flag=0,q,ma=-1;
  f();
  for (int i=1;i<=cnt;i++){
	if (prime[i]*prime[i]*prime[i+1]*prime[i+1]>23333333333333) break;  //避免爆掉long long
	for (int j=i+1;j<=cnt;j++){
		if (prime[i]*prime[i]*prime[j]*prime[j]<2333) continue;
		else if (prime[i]*prime[i]*prime[j]*prime[j]>23333333333333) break;
		cnt1++;
	}
  }
  cout<<cnt1;
  return 0;
}