这道题太精彩了!

我一开始想直接一波暴力算,然后叫上去只有50分,50分超时

然后我改成万位制提高运算效率,还是只有50分

然后我丧心病狂开long long用10的10次方作为一位,也就是100亿进制

去做,然后交上去多过了一个点,60分

附上丧心病狂的代码

#include<cstdio>
#include<cctype>
#include<algorithm>
#define REP(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std; typedef long long ll;
const int MAXN = 1123456 / 10;
const ll base = 1e10;
ll a[MAXN];
int len, n; void cal()
{
_for(i, 1, len) a[i] <<= 1;
_for(i, 1, len)
{
a[i+1] += a[i] / base;
a[i] %= base;
}
if(a[len+1])
{
len++;
a[len+1] += a[len] / base;
a[len] %= base;
}
} int main()
{
scanf("%d", &n);
len = 1; a[1] = 2;
REP(i, 1, n) cal(); int p = 1;
a[p] -= 1;
while(a[p] < 0)
{
a[p + 1]--;
a[p] += base;
p++;
}
if(a[len] == 0) len--; ll t = a[len], num = 0;
while(t) t /= 10, num++; printf("%d\n", len * 10 - 10 + num);
for(int i = 50, t = 0; i >= 1; i--, t++)
{
if(t == 5) puts(""), t = 0;
printf("%010lld", a[i]);
} return 0;
}

然后我最后还是看了题解

然后看到算位数真的是给折服的

还有这种操作???

由于10的x次方的位数是x+1

然后可以把2的p次方转化成类似10的x次方来算位数

然后这个时候就用到了log

一波骚操作可以得出位数就是p * lg(2) + 1

牛逼!

然后算前500位。

我竟然没看出这是个快速幂??????

用高精度做快速幂还是头一次。

不过稍微改一下就好了。

#include<cstdio>
#include<cctype>
#include<cstring>
#include<cmath>
#include<algorithm>
#define REP(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std; typedef long long ll;
const int MAXN = 1123;
struct bignum
{
ll s[MAXN]; int len;
bignum() { memset(s, 0, sizeof(s)); len = 1; }
}; bignum operator * (const bignum& a, const bignum& b)
{
bignum c;
c.len = min(a.len + b.len - 1 , 500);
_for(i, 1, a.len)
_for(j, 1, b.len)
{
c.s[i+j-1] += a.s[i] * b.s[j];
c.s[i+j] += c.s[i+j-1] / 10;
c.s[i+j-1] %= 10;
}
if(c.s[c.len+1] && c.len < 500) c.len++;
return c;
} int main()
{
int b;
scanf("%d", &b);
printf("%d\n", (int)(b * log10(2)+ 1)); bignum res; res.s[1] = 1;
bignum a; a.s[1] = 2;
while(b)
{
if(b & 1) res = res * a;
b >>= 1;
a = a * a;
}
res.s[1] -= 1; for(int i = 500; i >= 1; i--)
{
if(i != 500 && i % 50 == 0) puts(""); //这个写法很简便
printf("%d", res.s[i]);
}
puts(""); return 0;
}

洛谷 P1045 麦森数 (快速幂+高精度+算位数骚操作)的更多相关文章

  1. 洛谷 P1045 麦森数

    题目描述 形如2^{P}-1的素数称为麦森数,这时P一定也是个素数.但反过来不一定,即如果P是个素数,2^{P}-1不一定也是素数.到1998年底,人们已找到了37个麦森数.最大的一个是P=30213 ...

  2. 洛谷P1045 麦森数

    题目描述 形如2^{P}-12 ​P ​​ −1的素数称为麦森数,这时PP一定也是个素数.但反过来不一定,即如果PP是个素数,2^{P}-12 ​P ​​ −1不一定也是素数.到1998年底,人们已找 ...

  3. NOIP2003 普及组 洛谷P1045 麦森数 (快速幂+高精度)

    有两个问题:求位数和求后500位的数. 求位数:最后减去1对答案的位数是不影响的,就是求2p的位数,直接有公式log10(2)*p+1; 求后500位的数:容易想到快速幂和高精度: 1 #includ ...

  4. P1045麦森数

    P1045麦森数 #include<iostream> #include <cmath> #include <cstring> const int maxn = 1 ...

  5. 洛谷 P1226 【模板】快速幂||取余运算

    题目链接 https://www.luogu.org/problemnew/show/P1226 题目描述 输入b,p,k的值,求b^p mod k的值.其中b,p,k*k为长整型数. 输入输出格式 ...

  6. 洛谷试炼场-简单数学问题-P1045 麦森数-高精度快速幂

    洛谷试炼场-简单数学问题 B--P1045 麦森数 Description 形如2^P−1的素数称为麦森数,这时P一定也是个素数.但反过来不一定,即如果PP是个素数,2^P-1 不一定也是素数.到19 ...

  7. 【题解】[P1045] 麦森数

    题目 题目描述 形如2^P-1的素数称为麦森数,这时P一定也是个素数.但反过来不一定,即如果P是个素数,2^P-1 不一定也是素数.到1998年底,人们已找到了37个麦森数.最大的一个是P=30213 ...

  8. 洛谷P1226 【模板】快速幂||取余运算

    题目描述 输入b,p,k的值,求b^p mod k的值.其中b,p,k*k为长整型数. 输入输出格式 输入格式: 三个整数b,p,k. 输出格式: 输出“b^p mod k=s” s为运算结果 S1: ...

  9. 洛谷P1313 计算系数【快速幂+dp】

    P1313 计算系数 题目描述 给定一个多项式(by+ax)^k,请求出多项式展开后x^n*y^m 项的系数. 输入输出格式 输入格式: 输入文件名为factor.in. 共一行,包含5 个整数,分别 ...

随机推荐

  1. Python——Day4(基础知识练习二)

    # 1.请用代码实现:利用下划线将列表的每一个元素拼接成字符串.# li = ['alex','eric','rain']# li = ['alex','eric','rain']# li2 = &q ...

  2. iOS开发——自动填充短信验证码

    苹果在iOS 12,改进了一个很人性化的小细节.在做短信验证码功能的时候,自动获取短信中的验证码,然后点击填充即可.不用再向之前那样麻烦,自己看到弹出的短信信息后,死记硬背,再一个个敲上去,害怕背错了 ...

  3. pycharm 2018 3.4 for mac破解

    使用pycharm的小伙伴都知道,pycharm分为社区版和专业版,这里具体区别不作过多介绍.本文带大家安装mac版的2018 pycharm 3.4 1.去官网下载pycharm 3.4 for m ...

  4. thinkPHP利用ajax异步上传图片并显示、删除

    近来学习tp5的过程中,项目中有个发帖功能,选择主题图片.如下: 利用原始的文件上传处理,虽然通过原始js语句能实时显示上传图片,但是这样的话会涉及很多兼容问题.使用ajax技术,实现选择性删除所选图 ...

  5. 2019-03-20 Python爬取需要登录的有验证码的网站

    当你向验证码发起请求的时候,就有session了,记录下这次session 因为每当你请求一次验证码 或者 请求一次登录首页,验证码都在变动 验证码的链接可能不是固定的,可能需要GET/POST请求, ...

  6. 自备LocalDateTime工具类

    package cn.zytao.taosir.common.utils; import java.time.Instant; import java.time.LocalDate; import j ...

  7. 关于@SuppressWarnings("unchecked")注解

    解释一: 屏蔽某些编译时的警告信息         在强制类型转换的时候编译器会给出警告        加上程序代码        @SuppressWarnings("unchecked& ...

  8. apk去广告工具(利用apktool去除apk文件里的广告)

    基本知识 apk安装包的文件结构 以知名桌面软件“LauncherPro”为例,apk安装包文件目录: 文件目录如下: - META-INF - res - anim - color - drawab ...

  9. 使用Modernizr检测支持CSS3

    使用Modernizr检测支持CSS3 如果支持某个属性,会增加一个class,名字就是该属性: 不支持,名字是no-某属性 还提供了一个全局Modernizr对象,使用如下: <script ...

  10. UVA 1329 Corporative Network【并查集】

    题目链接: https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_proble ...