The Embarrassed Cryptographer

题意

  • 给一个两个素数乘积(1e100)K, 给以个数L(1e6), 判断K的两个素数是不是都大于L

题解

  • 对于这么大的范围,素数肯定是要打表(可采用埃筛,欧拉筛,莫比乌斯筛);这里有别人模板
  • 简单的想法是遍历表中<L的素数去模K(即便对K分解也是如此办的)
  • 但K很大需要高精度取模,由于足够大需要转换成K,L都足够大需要转换成千进制

代码

素数打表

int isp[maxn];// isp[i]=0 i是素数
int su[maxn];// su[i]=p 第i个素数是p
int cnt;
void get_prime(){
mm0(isp);
cnt=0;
for(int i=2;i<=maxn;i++){
if(isp[i]==0){
su[++cnt]=i;
}
for(int j=1;su[j]<su[cnt]&&su[j]*i<=maxn;j++){
isp[su[j]*i]=1;// 非素数
if(i%su[j]==0) break;//一个优化,每个被删掉的合数都是以最小素因子标记,这样可以避免重复删掉合数
}
}
//确定能满足循环(找比L小的素数)结束
su[++cnt]=1e9+7;
}

高精度取模

基本原理:

(a+b)%c=(a%c+b%c)%c (1)

(ab)%c=((a%c)(b%c))%c (2)

对于(2) 式这样也成立 (ab)%c=((a%c)b)%c

证明可通过a=ck1+t1,b=ck2+t2 形式证明

所以如123%13 可以先求 1%13,然后求12%13,最后求123%13

我们这里 展示十进制情况:

int get_mod(int p){
int m=0;
for(int i=0;i<strlen(s);i++){
m=((m*10)%p+(s[i]-'0')%p)%p;
}
return m;
}
//值得一提的是 s[0]式最位,对这份代码进行改变时,需注意自己的存储哪一位时高位

题目的AC代码

#include <bits/stdc++.h>
typedef long long ll;
#define mm0(x) memset(x,0,sizeof(x))
#define mm1(x) memset(x,0x1f,sizeof(x))
#define mm2(x) memset(x,0x3f,sizeof(x))
#define mm3(x) memset(x,0xff,sizeof(x))
const int maxn = 1e6 + 5;
const int maxm = 1e5+5;
//author:fridayfang
//date:18.8月 21
//global varibles
//大数(高精度)取模+素数筛法
char s[102];//所有因子 K 10^100
int kt[34];
int isp[maxn];// isp[i]=0 i是素数
int su[maxn];// su[i]=p 第i个素数是p
int cnt;
void get_prime(){
mm0(isp);
cnt=0;
for(int i=2;i<=maxn;i++){
if(isp[i]==0){
su[++cnt]=i;
}
for(int j=1;su[j]<su[cnt]&&su[j]<=maxn/i;j++){
isp[su[j]*i]=1;// 非素数
if(i%su[j]==0) break;
}
}
//确定能满足循环(找比L小的素数)结束
su[++cnt]=1e9+7;
}
//get prime
int get_mod(int p,int po){
int m=0;
for(int i=po;i>=1;i--){
m=((m*1000)%p+(kt[i])%p)%p;
}
return m;
}
/*
int get_mod(int p){
int m=0;
for(int i=0;i<strlen(s);i++){
m=((m*10)%p+(s[i]-'0')%p)%p;
}
return m;
}
*/
void read(){
int L;
while(true){
scanf("%s",s);
scanf("%d",&L);
if(L==0) break;
//转换为千进制
int po=0;
int len=strlen(s);
int j;
for(j=len-1;j-2>=0;j=j-3){
kt[++po]=((s[j]-'0')+(s[j-1]-'0')*10+(s[j-2]-'0')*100);
}
po++;
kt[po]=0;
int mul=1;
while(j>=0){
kt[po]+=(s[j]-'0')*mul;
j--;
mul=mul*10;
}
//printf("kt%d %d\n",kt[1],kt[2]); int i=0;
bool good=true;
while(su[++i]<L){
//printf("mod %s %d %d\n",s,su[i],get_mod(su[i],po));
if(get_mod(su[i],po)==0){good=false;break;}
}
if(good) printf("GOOD\n");
else printf("BAD %d\n",su[i]);
}
} int main(){
get_prime();
read();
return 0;
}

需要对高精度和数论进行补题

HDU 2303 The Embarrassed Cryptographer的更多相关文章

  1. POJ 2635 The Embarrassed Cryptographer

    大数取MOD... The Embarrassed Cryptographer Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 1 ...

  2. (POJ2635)The Embarrassed Cryptographer(大数取模)

    The Embarrassed Cryptographer Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 13041 Accep ...

  3. POJ2635——The Embarrassed Cryptographer(高精度取模+筛选取素数)

    The Embarrassed Cryptographer DescriptionThe young and very promising cryptographer Odd Even has imp ...

  4. poj2635The Embarrassed Cryptographer(同余膜定理)

    The Embarrassed Cryptographer Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 15069   A ...

  5. [ACM] POJ 2635 The Embarrassed Cryptographer (同余定理,素数打表)

    The Embarrassed Cryptographer Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 11978   A ...

  6. POJ 2635 The Embarrassed Cryptographer (千进制,素数筛,同余定理)

    The Embarrassed Cryptographer Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 15767   A ...

  7. The Embarrassed Cryptographer(高精度取模+同余模定理)

    Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 11435   Accepted: 3040 Description The ...

  8. POJ - 2635 E - The Embarrassed Cryptographer

    The young and very promising cryptographer Odd Even has implemented the security module of a large s ...

  9. POJ2635The Embarrassed Cryptographer(大数取余+素数筛选+好题)

    题目链接 题意:K是由两个素数乘积,如果最小的素数小于L,输出BAD最小的素数,否则输出GOOD 分析 素数打表将 L 大点的素数打出来,一定要比L大,然后就开始枚举,只需K对 素数 取余 看看是否为 ...

随机推荐

  1. 路飞学城Python-Day136

    列举Http请求中常见的请求方式 根据HTTP标准,HTTP请求可以使用多种请求方法. HTTP1.0定义了三种请求方法: GET, POST 和 HEAD方法. HTTP1.1新增了五种请求方法:O ...

  2. 基于fullpage的自动播放,手动播放,暂停页面的功能

    功能如下: 1.默认加载方式为“自动播放 ”方式,即从第1屏至第5屏 页面循环加载显示,每屏每次仅显示1个页面,页面间停留时间为“10”秒2.手动播放过程中,按数字键“1”-“5”,将直接切到指定页面 ...

  3. PostgreSQL 安装配置 (亲测可用)

    转自:http://blog.csdn.net/jesseyoung/article/details/41348835 受作者博客限制,请访问上面的链接 ---------- 下面是另一个转载 --- ...

  4. [网络流24题] 最长k可重区间集问题 (费用流)

    洛谷传送门 LOJ传送门 很巧妙的建图啊...刚了$1h$也没想出来,最后看的题解 发现这道题并不类似于我们平时做的网络流题,它是在序列上的,且很难建出来二分图的形. 那就让它在序列上待着吧= = 对 ...

  5. centos solr 部署到 tomcat 上

    一.安装 java1.7 环境 链接:https://pan.baidu.com/s/1ti6j9jD-RwUN5xl3bc3ZDw 密码:oc9a 二.下载 tomcat 并解压 链接:https: ...

  6. selenium+java实现浏览器前进、后退和刷新

  7. 2015 Multi-University Training Contest 10 hdu 5411 CRB and Puzzle

    CRB and Puzzle Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)To ...

  8. Elasticsearch 三种分页方式

    from + size 浅分页 "浅"分页可以理解为简单意义上的分页.它的原理很简单,就是查询前20条数据,然后截断前10条,只返回10-20的数据.这样其实白白浪费了前10条的查 ...

  9. MySQL 函数大全及用法示例

    1.字符串函数ascii(str)   返回字符串str的第一个字符的ascii值(str是空串时返回0)  mysql> select ascii('2');  -> 50  mysql ...

  10. TNS-12557: TNS:protocol adapter not loadable TNS-12560: TNS:protocol adapter error

    Description: Oracle 10.2 on hpux 11.23 PA. When i try to start listener i go the next errors: Error ...