<NOIP> 25 . P1217 [USACO1.5]回文质数 Prime Palindromes

最新推荐文章于 2020-09-21 10:58:09 发布
最新推荐文章于 2020-09-21 10:58:09 发布
阅读量312 收藏
点赞数
分类专栏： OJ 文章标签： NOIP OJ c++ 源代码 编程
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u010202481/article/details/74160146
版权

题解：这是洛谷的第25道题目，关于这个题，我的方法有点Low，就是第一构造回文数，第二素性测试（虽然AC了，但是感觉好丢人）。

源代码

#include <iostream>
#include <stdlib.h>
#include <math.h>
#include <vector>

using namespace std;

#define  DATA_TYPE    unsigned long long

DATA_TYPE The_Multiplication_Operator(DATA_TYPE a, DATA_TYPE b, DATA_TYPE mod)
{
    DATA_TYPE ans = 0;
    while (b)
    {
        if (b & 1)
        {
            b--;
            ans = (ans + a) % mod;
        }
        b /= 2;
        a = (a + a) % mod;
    }
    return ans;
}

DATA_TYPE Addition_Operator(DATA_TYPE a, DATA_TYPE b, DATA_TYPE mod)
{
    DATA_TYPE ans = 1;
    while (b)
    {
        if (b & 1)
        {
            ans = The_Multiplication_Operator(ans, a, mod);
        }

        b /= 2;
        a = The_Multiplication_Operator(a, a, mod);
    }
    return ans;
}

bool Inspection(DATA_TYPE a, DATA_TYPE n)
{
    DATA_TYPE tem = n - 1;
    int j = 0;
    while (tem % 2 == 0)
    {
        tem /= 2;
        j++;
    }

    DATA_TYPE x = Addition_Operator(a, tem, n);
    if (x == 1 || x == n - 1)
        return true;
    while (j--)
    {
        x = The_Multiplication_Operator(x, x, n);
        if (x == n - 1) return true;
    }

    return false;
}

bool Miller_Rabin(DATA_TYPE Target_Number)
{
    if (Target_Number == 2)
        return true;
    if (Target_Number < 2 || Target_Number % 2 == 0)
        return false;

    for (int i = 1; i <= 20; i++)
    {
        DATA_TYPE a = rand() % (Target_Number - 2) + 1;
        if (!Inspection(a, Target_Number))
            return false;
    }

    return true;
}

int main()
{
    int number = 0;
    unsigned long long a, b, c;
    cin >> a >> b;

    c = b;
    while (c != 0)
    {
        c /= 10;
        number++;
    }

    switch (number)
    {
    case 1:
    {
        if(a<=5)
            cout << "5" << endl;
        if(a<=7)
            cout << "7" << endl;

        break;
    }
    case 2:
    {
        if (a <= 5)
            cout << "5" << endl;
        if (a <= 7)
            cout << "7" << endl;
        if(a<=11)
            cout << "11" << endl;

        break;
    }
    case 3:
    {
        if (a <= 5)
            cout << "5" << endl;
        if (a <= 7)
            cout << "7" << endl;
        if (a <= 11)
            cout << "11" << endl;

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 100 * d1 + 10 * d2 + d1;
                if (palindrome>=a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        break;
    }
    case 4:
    {
        if (a <= 5)
            cout << "5" << endl;
        if (a <= 7)
            cout << "7" << endl;
        if (a <= 11)
            cout << "11" << endl;

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 100 * d1 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 1000 * d1 + 100 * d2 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        break;
    }
    case 5:
    {
        if (a <= 5)
            cout << "5" << endl;
        if (a <= 7)
            cout << "7" << endl;
        if (a <= 11)
            cout << "11" << endl;

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 100 * d1 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 1000 * d1 + 100 * d2 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    long long palindrome = 10000 * d1 + 1000 * d2 + 100 * d3 + 10 * d2 + d1;
                    if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                        cout << palindrome << endl;
                }
            }
        }

        break;
    }
    case 6:
    {
        if (a <= 5)
            cout << "5" << endl;
        if (a <= 7)
            cout << "7" << endl;
        if (a <= 11)
            cout << "11" << endl;

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 100 * d1 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 1000 * d1 + 100 * d2 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    long long palindrome = 10000 * d1 + 1000 * d2 + 100 * d3 + 10 * d2 + d1;
                    if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                        cout << palindrome << endl;
                }
            }
        }
        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    long long palindrome = 100000 * d1 + 10000 * d2 + 1000 * d3 + 100 * d3 + 10 * d2 + d1;
                    if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                        cout << palindrome << endl;
                }
            }
        }

        break;
    }
    case 7:
    {
        if (a <= 5)
            cout << "5" << endl;
        if (a <= 7)
            cout << "7" << endl;
        if (a <= 11)
            cout << "11" << endl;

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 100 * d1 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 1000 * d1 + 100 * d2 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    long long palindrome = 10000 * d1 + 1000 * d2 + 100 * d3 + 10 * d2 + d1;
                    if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                        cout << palindrome << endl;
                }
            }
        }
        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    long long palindrome = 100000 * d1 + 10000 * d2 + 1000 * d3 + 100 * d3 + 10 * d2 + d1;
                    if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                        cout << palindrome << endl;
                }
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    for (int d4 = 0; d4 <= 9; d4++)
                    {
                        long long palindrome = 1000000 * d1 + 100000 * d2 + 10000 * d3 + 1000 * d4 + 100 * d3 + 10 * d2 + d1;
                        if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                            cout << palindrome << endl;
                    }
                }
            }
        }

        break;
    }

    case 8:
    case 9:
    {
        if (a <= 5)
            cout << "5" << endl;
        if (a <= 7)
            cout << "7" << endl;
        if (a <= 11)
            cout << "11" << endl;

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 100 * d1 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                long long palindrome = 1000 * d1 + 100 * d2 + 10 * d2 + d1;
                if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                    cout << palindrome << endl;
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    long long palindrome = 10000 * d1 + 1000 * d2 + 100 * d3 + 10 * d2 + d1;
                    if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                        cout << palindrome << endl;
                }
            }
        }
        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    long long palindrome = 100000 * d1 + 10000 * d2 + 1000 * d3 + 100 * d3 + 10 * d2 + d1;
                    if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                        cout << palindrome << endl;
                }
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    for (int d4 = 0; d4 <= 9; d4++)
                    {
                        long long palindrome = 1000000 * d1 + 100000 * d2 + 10000 * d3 + 1000 * d4 + 100 * d3 + 10 * d2 + d1;
                        if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                            cout << palindrome << endl;
                    }
                }
            }
        }

        for (int d1 = 1; d1 <= 9; d1 += 2)
        {
            for (int d2 = 0; d2 <= 9; d2++)
            {
                for (int d3 = 0; d3 <= 9; d3++)
                {
                    for (int d4 = 0; d4 <= 9; d4++)
                    {
                        long long palindrome = 10000000 * d1 + 1000000 * d2 + 100000 * d3 + 10000 * d4 + 1000 * d4 + 100 * d3 + 10 * d2 + d1;
                        if (palindrome >= a && palindrome <= b && Miller_Rabin(palindrome))
                            cout << palindrome << endl;
                    }
                }
            }
        }

        break;
    }

    default:
        break;
    }
    
    system("pause");
    return 0;
}


夜里的柚子
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
洛谷题目目录
zqhf123的博客
02-07 25万+
#1 新手村 关卡1-1 洛谷的第一个任务 超级玛丽游戏（洛谷-P1000）：点击这里 A+B Problem（洛谷-P1001）：点击这里 小玉买文具（洛谷-P1421）：点击这里 小鱼的游泳时间（洛谷-P1425）：点击这里 关卡1-2 顺序与分支 小玉家的电费（洛谷-P1422）：点击这里 不高兴的津津（洛谷-P1085）：点击这里 津津的储蓄计划（洛谷-P1089）：点击这里 买铅笔（洛谷...
P1217 [USACO1.5]回文质数 Prime Palindromes
，，，，，，，，，的博客
07-23 157
题目描述 因为 151 既是一个质数又是一个回文数（从左到右和从右到左是看一样的），所以 151 是回文质数。 写一个程序来找出范围[a,b](5≤a<b≤100,000,000)[a,b](5≤a<b≤100,000,000)( 一亿)间的所有回文质数。 输入输出格式 输入格式： 第 1 行: 二个整数 a 和 b . 输出格式： 输出一个回文质数的列表，一行一...
洛谷题解 P1015回文数（高精度+模拟）
Enternalwiser的博客
09-21 900
洛谷 P1015回文数
洛谷3258：[USACO2012 MAR]Flowerpot 花盆——题解
weixin_33747129的博客
03-12 156
https://www.luogu.org/problemnew/show/P2698#sub 老板需要你帮忙浇花。给出N滴水的坐标，y表示水滴的高度，x表示它下落到x轴的位置。 每滴水以每秒1个单位长度的速度下落。你需要把花盆放在x轴上的某个位置，使得从被花盆接着的第1滴水开始，到被花盆接着的最后1滴水结束，之间的时间差至少为D。 我们认为，只要水滴落到x轴上，与花盆的边沿对齐，就认为...
USACO1.5]回文质数 Prime Palindromes
LjcoderDCC的博客
12-05 4130
题目描述 因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。 写一个程序来找出范围[a,b](5 输入输出格式 输入格式： 第 1 行: 二个整数 a 和 b . 输出格式： 输出一个回文质数的列表，一行一个。 输入输出样例 输入样例#1： 5 500 输出样例#1： 5 7 11
[USACO1.5]回文质数 Prime Palindromes
Bond0909的博客
04-25 366
题目描述因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。写一个程序来找出范围[a,b](5 &lt;= a &lt; b &lt;= 100,000,000)( 一亿)间的所有回文质数;输入输出格式输入格式：第 1 行: 二个整数 a 和 b .输出格式：输出一个回文质数的列表，一行一个。如果直接进行枚举，则显然会超时，显然需要进行剪枝。首先，可以通...
usaco_A.rar_NOI_NOIP_usaco
09-23
【标题】"USACO_A.rar_NOI_NOIP_usaco" 涉及的主要知识点是USACO（美国计算机奥林匹克）的编程题库，这部分内容是针对NOI（全国青少年信息学奥林匹克）和NOIP（全国青少年信息学奥林匹克联赛）的训练资料，适合对ACM...
extend-noip：延长noip.com上的dns到期日期
02-25
扩展noip 在noip.com上延长DNS过期日期 在noip.com上自动延长dns /域名的有效期 分叉此仓库 设置结果回购Secrets ...python -m extend-noip -u your_noip_username -p password 或直接使用extend-no
NOIP-Junior.rar
07-20
《全国青少年信息学奥林匹克联赛（NOIP）普及组历年复赛试题详解与解析》 全国青少年信息学奥林匹克联赛（National Olympiad in Informatics in Provinces，简称NOIP）是由中国计算机学会主办的一项针对中学生的...
NOIP普及组近5年NOIP试题分析.ppt
07-21
NOIP普及组近5年NOIP试题分析，全国青少年信息学奥林匹克竞赛
NOIP模板1.zip
05-12
7. **大质数.cpp**：这个文件可能包含了生成大质数的方法，如 Miller-Rabin素数测试或者AKS素数测试，用于处理需要大整数质数的题目。 这些模板覆盖了信息学竞赛中的核心算法，包括数值计算、图论、字符串处理和...
1254. 【USACO题库】1.5.2 Prime Palindromes回文质数 (Standard IO)
happyfor_me的博客
08-13 465
很多同志在看到这道题时不是傻就是呆了，这么大的数据！这么恶心的数据！其实这题可以打表--------先用超时的程序输出正确答案到文件里，我好机智把文件里全部东东放到一个数组里，然后用repeat循环，直到ai大于给定的b为止。程序如下： Var sheet : array [1..779] of longint = (5,7,11,101,131,151,181,191,313,353,3
回文质数 Prime Palindromes
wzw-yali的博客
07-12 570
题目描述 因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。 写一个程序来找出范围[a,b](5 输入输出格式 输入格式： 第 1 行: 二个整数 a 和 b . 输出格式： 输出一个回文质数的列表，一行一个。 输入输出样例 输入样例#1： 5 500 输出样例#1： 5 7 11 101 131 151 181 1
[P1217]回文质数
今天的北宅依旧没有回家
06-09 403
原题链接这次的数据还是挺亲切的 最暴力的写法拿了8个点最后一个点怎么也不过 然后看题解 Duang 忽然想起来这题讲过因为4/6/8位是没有回文质数的 [4/6/8位的回文数都是11的倍数] 所以只需要对其他位数的回文数打表 判断是否为质数但是又忘了还有上下界 所以又白交了一次总之还是过了啦#include<iostream> #include<cstring> #include<a
洛谷 P1217 [USACO1.5]回文质数 Prime Palindromes
邵光亮的博客
05-29 704
题目描述 因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。 写一个程序来找出范围[a,b](5 <= a < b <= 100,000,000)( 一亿)间的所有回文质数; 输入输出格式 输入格式： 第 1 行: 二个整数 a 和 b . 输出格式： 输出一个回文质数的列表，一行一个。 输入输...
洛谷 P1217
qq670471201的博客
07-31 1275
题目链接如下https://www.luogu.org/problemnew/show/P1217 找出a到b 是回文数并且是素数的数。 这还不简单 前两天刚学的线性筛素数。卡卡卡卡一顿猛敲。最后一个测试点TLE了。我无语了。。。。。。。 改吧。 自古评论区里出人才。 判断素数也有几个小结论。如下： 1.偶数位数回文数（除11）必定不是质数（自行百度），所以只要运行到100000...
回文质数
青耕
07-20 4176
回文质数 本题地址： http://www.luogu.org/problem/show?pid=1217题目描述 因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。 写一个程序来找出范围[a,b](5 <= a < b <= 100,000,000)( 一亿)间的所有回文质数; 输入输出格式输入格式： 第 1 行: 二个整数 a 和 b .输出格
<NOIP> 4 . P1425 小鱼的游泳时间
Charmain
06-29 3382
题解：这是洛谷的第四题，考察的是对一个标准时间的时间差的计算，有兴趣的话，可以写一个关于时间的函数或者模板类。 源代码： #include //#include using namespace std; int main() { long long a, b, c, d, e, f; cin >> a >> b >> c >> d;
P1015 [NOIP1999 普及组] 回文数
最新发布
03-14
回文数是指从左到右和从右到左读起来都一样的数字。现在给定一个数字，你需要判断它是否是回文数。 解题思路： 将数字转换为字符串，然后判断字符串是否是回文串即可。 代码实现： #include <iostream> #include...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值