通过观察发现:
1. 1肯定不符合题目给定的条件。
2. 对于每个质数i,一定可以和i*i组成满足题意的数对。
3. 对于非质数,一定可以和其他数字组成满足题意的数对。
显然,我们需要筛出所有质数。由于题目的范围太大,所以我们可以预处理出来答案,对于每次询问,直接输出答案即可。预处理可以使用差分数组来完成。
注意一定要使用快读。
时间复杂度为O(n+t)
#include<bits/stdc++.h>
#define x first
#define y second
#define PB push_back
#define mst(x,a) memset(x,a,sizeof(x))
#define all(a) begin(a),end(a)
#define rep(x,l,u) for(ll x=l;x<u;x++)
#define rrep(x,l,u) for(ll x=l;x>=u;x--)
#define sz(x) x.size()
#define IOS ios::sync_with_stdio(false);cin.tie(0);
using namespace std;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef pair<long,long> PLL;
typedef pair<char,char> PCC;
typedef long long ll;
const int N=1e6+10;
const int M=1e6+10;
const int INF=0x3f3f3f3f;
const int MOD=1e9+7;
/*独立思考
不要看测试样例
找性质
试着证明
写完后不盲目交,检查*/
int ans[N],prime[N];
bool st[N];
void solve(){
int n;scanf("%d",&n);
printf("%d\n",ans[n]);
}
int main(){
//IOS;
//freopen("test.txt", "r", stdin);
//freopen("test.txt", "w", stdout);
int cnt=0;
rep(i,2,N+1){
if(!st[i]){
prime[cnt++]=i;
}
for(int j=0;prime[j]<=N/i;j++){
st[i*prime[j]]=true;
if(i%prime[j]==0) break;
}
}
ans[1]=1;
rep(i,2,N+1){
if(st[i]) continue;
ans[i]++;
if(i*i<=N) ans[i*i]--;
}
rep(i,1,N+1){
ans[i]+=ans[i-1];
}
int t;scanf("%d",&t);
while(t--)
solve();
return 0;
}
原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。
如有侵权,请联系 cloudcommunity@tencent.com 删除。
原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。
如有侵权,请联系 cloudcommunity@tencent.com 删除。