POJ2689——区间筛法——Prime Distance
http://poj.org/problem?id=2689
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 | /* 区间筛法适用于L, R范围比较大,但是区间长度较小时 套个模板然后取出最小最大就可以 */ /************************************************ * Author :Powatr * Created Time :2015-8-17 9:38:50 * File Name :POJ2689.cpp ************************************************/ #include <cstdio> #include <algorithm> #include <iostream> #include <sstream> #include <cstring> #include <cmath> #include <string> #include <vector> #include <queue> #include <deque> #include <stack> #include <list> #include <map> #include <set> #include <bitset> #include <cstdlib> #include <ctime> using namespace std; #define lson l, mid, rt << 1 #define rson mid + 1, r, rt << 1 | 1 typedef long long ll; const int MAXN = 1e6 + 1e4; const int INF = 0x3f3f3f3f; const int MOD = 1e9 + 7; const int N = 50010; bool prime[MAXN]; bool sieve[N]; int len; void segent_sieve(ll L, ll R) { len = R - L + 1; for ( int i = 0; i < len; i++) prime[i] = true ; if (1 - L >= 0) prime[1-L] = false ; //特判 for (ll i = 2; i*i < R; i++){ if (seive[i]){ for (ll j = max(1ll*2, (L-1+i)/i)*i; j <= R; j+=i) //(L-1+i/i)得到最接近L的i的倍数 prime[j-L] = false ; //偏移量,因为单纯j会导致数组越界,把范围弄到0到len上 } } } int main(){ //艾氏筛法 O(nlognlogn) for ( int i = 2; i < N; i++) sieve[i] = true ; for ( int i = 2; i*i < N; i++){ if (sieve[i]){ for ( int j = i*2; j < N; j+=i){ sieve[j] = false ; } } } int L, R; while (~ scanf ( "%d%d" , &L, &R)){ segent_seive(L, R); int mmax, mmin; int lmax, rmax, lmin, rmin; mmax = -1, mmin = 1 << 30; int t = -1; for ( int i = 0; i < len; i++){ if (prime[i]){ if (t >= 0){ if (mmax < i - t) mmax = i - t, lmax = L + t, rmax = L + i; if (mmin > i - t) mmin = i - t, lmin = L + t, rmin = L + i; t = i; } else t = i; } } if (mmax != -1 || mmin != 1 << 30) printf ( "%d,%d are closest, %d,%d are most distant.\n" , lmin, rmin, lmax, rmax); else printf ( "There are no adjacent primes.\n" ); } return 0; } |
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 智能桌面机器人:用.NET IoT库控制舵机并多方法播放表情
· Linux glibc自带哈希表的用例及性能测试
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· 新年开篇:在本地部署DeepSeek大模型实现联网增强的AI应用
· Janus Pro:DeepSeek 开源革新,多模态 AI 的未来
· 互联网不景气了那就玩玩嵌入式吧,用纯.NET开发并制作一个智能桌面机器人(三):用.NET IoT库
· 【非技术】说说2024年我都干了些啥