【数论+技巧】神奇的Noip模拟试题第二试 T1 素数统计

1.      素数统计

(pcount.pas/.c/.cpp)

【问题描述】

      小tan的老师揣谙戈给同学们布置了一道题,要求统计给定区间内素数的个数。“这不是很简单吗?”小tan忍不住说。揣谙戈冷笑一下说:“等你们看到题目就知道了。”便转身离去。

      果然,小tan被那极大的区间吓怕了,现在是你拯救她的时候。

【输入】

      输入文件名为pcount.in。

      输入一行两个正整数a和b,表示给定区间为[a,b]。

【输出】

      输出文件名为pcount.out。

      输出一个整数,表示区间内素数数量。

 

【输入输出样例】

pcount.in

pcount.out

1 17

7

 

【数据范围】

      对于30%的数据,有n<m≤1,000;

      对于60%的数据,有n<m≤1,000,000;

      对于100%的数据,有n<m<2^31,m-n≤1,000,000。

 

解题报告:

  看到这么大的数据,果然是被吓怕了。前一天才学的数论,这一次考试就考了两道题,虽然我只做了两道题。先说说这道题吧,这道题自己做的时候心想,反正都得不全分,干脆就把数组定小一点,得个六七十分吧。于是开始暴力筛数法,把1~m和1~n的素数都筛出来,然后相减,只不过要注意n是素数的情况,ans--;于是得了70。

  正解:这么大的数组肯定不行,所以需要下标(数组)平移,把n~m之间的及大于sqrt(m)的平移到0~m-n的下标。然后扫一遍这个区间,把ans加上。但要注意特殊情况:n==1的时候,ans--;

  值得注意的是:循环时把 j 定义为long long,因为当m很大时,i * i+ i 就超出了int的范围,这个时候会变为负数,然后无限循环,最后一组数据就超时了。

  代码附下:

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cmath>
 4 using namespace std;
 5 long long n,m,ans;
 6 const int maxn=1000005;
 7 bool vis1[maxn],vis2[maxn];
 8 int main()
 9 {
10     freopen("pcount.in","r",stdin);
11     freopen("pcount.out","w",stdout);
12     cin>>n>>m;
13     long long mid=sqrt(m)+1;
14     for (long long i=2;i<=mid;i++)
15     {
16         if (!vis1[i])
17         {
18             for (long long j=i*i;j<=mid;j+=i)//long long
19                 vis1[j]=1;
20             for (long long j=n/i*i;j<=m;j+=i)
21              if (j>i&&j>=n) vis2[j-n]=1;
22         }    
23     }
24     for (int i=0;i+n<=m;i++)
25        if (!vis2[i]) ans++;
26     if (n==1) ans--;
27     cout<<ans;
28     return 0;
29 }

 

转载于:https://www.cnblogs.com/lx0319/p/5690378.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值