Goldbach

歌德巴赫猜想,暴力枚举。
事先打好1e5的素数表,然后用米勒罗宾判断枚举n-prime[i[是不是素数。
这里米勒罗宾判断素数网上都有模板的,具体就不解释啦(窝是不会告诉泥萌窝也不会的)

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <ctime>
#include <cmath>
using namespace std;
typedef unsigned long long  LL;const int N = 100000,M = 10000,count1 = 10;
int prime1[M], k = 0;
bool book[N];LL prime[6] = {2, 3, 5, 233, 331};
LL qmul(LL x, LL y, LL mod) { return (x * y - (long long)(x / (long double)mod * y + 1e-3) *mod + mod) % mod;
}
LL qpow(LL a, LL n, LL mod) {LL ret = 1;while(n) {if(n & 1) ret = qmul(ret, a, mod);a = qmul(a, a, mod);n >>= 1;}return ret;
}
bool Miller_Rabin(LL p) {if(p < 2) return 0;if(p != 2 && p % 2 == 0) return 0;LL s = p - 1;while(! (s & 1)) s >>= 1;for(int i = 0; i < 5; ++i) {if(p == prime[i]) return 1;LL t = s, m = qpow(prime[i], s, p);while(t != p - 1 && m != 1 && m != p - 1) {m = qmul(m, m, p);t <<= 1;}if(m != p - 1 && !(t & 1)) return 0;}return 1;
}
//以上均为米勒罗宾模板
void init()
{for(int i = 2; i < N; ++i){if(book[i]==false){prime1[k++] = i;}for(int j = 2; i*j <= N; ++j){if(book[i*j]==false)book[i*j] = true;}}
}int main()
{init();int t;LL n;cin>>t;while(t--){scanf("%llu",&n);int flag = 0;for(int i = 0; i < k; ++i){if(n-prime1[i]>=0&&Miller_Rabin(n-prime1[i])){printf("%d %llu\n",prime1[i],n-prime1[i]);flag = 1;break;}}}return 0;
}

2018 ACM-ICPC 宁夏预选赛网赛 B-Goldbach相关推荐

  1. 2018 ACM/ICPC 北京赛区网络赛 D 80 Days 线段树

    http://hihocoder.com/problemset/problem/1831?sid=1390457 描述 80 Days is an interesting game based on ...

  2. 解题报告(一)F、(2018 ACM - ICPC shenyang I)Distance Between Sweethearts(数学期望 + 乘法原理 + FWT)(4.5)

    繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量题解和代码,题目难度不一定按照题号排序,我会在每道题后面加上题目难度指数(1∼51 \sim 51∼5),以模板题难度 11 ...

  3. 2014 ACM/ICPC 北京赛区网络赛解题报告汇总

    首页 算法竞赛» 信息聚合 ONLINE JUDGE 书刊杂志 BLOG» 新闻故事» 招聘信息» 投稿须知 2014 ACM/ICPC 北京赛区网络赛解题报告汇总 九月 21st, 2014 | P ...

  4. 2011 ACM/ICPC 福州赛区网络赛解题报告

    第一次写网络赛的题解,福州赛区网络赛作为我第一年ACM最后一次网络赛酱油,画了一个很像逗号的句号.....好吧,还得为北京现场赛准备啊准备....... 这次酱油打的很犀利,貌似出第一题很快,之后节奏 ...

  5. 2018宁夏邀请赛网赛 G.Trouble of Tyrant(单调栈)

    题意: 有n个点,2n-3条边的图.点 1 到每个点有一条边,编号相邻的两个点有一条边.q次询问,每次询问一个增量d,问图中每条边都增加 d 后,1 到 n 的最短路是多少.增量独立,不累加. 1 & ...

  6. 2014 ACM/ICPC 鞍山赛区网络赛(清华命题)

    为迎接10月17号清华命题的鞍山现场赛 杭电上的题目 Biconnected(hdu4997)     状态压缩DP Rotate(hdu4998)    相对任一点的旋转 Overt(hdu4999 ...

  7. 2018 ACM/ICPC 沈阳站小结

    第二次现场赛,能拿金已经很开心了 热身赛 开幕式后的热身赛都是训练过的三个题,40min基本就解决了,剩下的时间就是在调试机器什么的 现场赛 上来队友先过了J的模拟,我看了看G题,发现可做就直接上了, ...

  8. 2018 ACM/ICPC 南京站小结

    第一次打现场赛,又是在半主场,其实挺激动的,但是结果不尽人意 热身赛 开幕式后的热身赛就打得有些问题,四个题都不难,我们都想到了做法,但因为机时不够+细节问题导致前两个小时一直没有过题,然后AC的三道 ...

  9. 2018 ACM ICPC Asia Regional - Seoul B.Cosmetic Survey

    参考大佬代码 题目大意 nnn个顾客,mmm种画妆品,每一个顾客会给mmm种化妆品一个值,这个值代表这个化妆品在他心中的排名,排名越小越喜欢,如果这个值为0说明最不喜欢这一种化妆品(值为0理解为无穷大 ...

最新文章

  1. python3.7.4安装教程桌面_Python 3.7.4 for Windows的安装
  2. Mac homebrew类似apt-get命令安装包
  3. jquery对象PHP转换,jquery对象和DOM对象如何相互转换?
  4. linux shell touch更新文件时间
  5. matlab中jet的例子,Python matplotlib.pyplot.jet()用法及代码示例
  6. GIS创新实践【课堂笔记】
  7. mysql hint 简书_MySQL
  8. html期末作品_声来不凡丨同学快来!这份期末考试卷由你来批改!
  9. 关于在Google Earth中动态加载地标问题
  10. 传奇服务器端地图链接在哪个文件夹,传奇客户端有新地图,怎样在服务端设置链接?让玩家进去?...
  11. 计算机中的标准差是哪个英语单词,标准差是什么意思
  12. echarts地图功能实现及坐标定位
  13. excel数据个数统计问题
  14. C语言采药动态规划,Vijos P1104 采药 动态规划
  15. 经纬度、平面坐标系转换方法
  16. React教程(五)——生命周期函数
  17. inno setup QQ音乐安装包界面
  18. 【C++】输入流和输出流
  19. 弹吉他节奏跟不上?试试用这个循序渐进的方法吧~
  20. ubantu 系统分区介绍

热门文章

  1. 几款常用LED显示屏驱动芯片介绍及选择VK1624
  2. 两寸标准证件照大小是多少?分享两种调整尺寸的方法
  3. 在KVM最小化搭建openstack平台 --快速部署openstack
  4. 数风流人物还看今朝|前后端分离微服务项目常用中间件以及指令
  5. 屏蔽各APP广告方法与广告路径(须Root)
  6. 计算机丢失quartz.dll什么意思,电脑系统提示没有找到quartz.dll该怎么办
  7. JAVA中接口存在的意义
  8. Unity 代码实现形成圆形及形成球形
  9. 业余选手常见的技术缺陷及解决方法(下)
  10. 当当网商品详情API接口(当当商品详情接口,当当商品评论接口,当当商品问答接口,当当抢购价接口,关键词搜索当当网商品接口)代码对接教程