麦森数(转)

 1 //形如2n-1的素数称为麦森数,这时n一定也是个素数。但反过来不一定,即如果n是个素数。2n-1不一定也是素数。
 2 
 3 #include<iostream>
 4  #include<cmath>
 5  #include<cstdio>
 6  #include<cstring>
 7  #define N 126
 8  using namespace std;
 9  int ans[N],anspow[N];
10  void mult(int ans[],int anspow[])
11  {
12    int i,j;
13    int c[N];
14    memset(c,0,sizeof(c));
15          for(i=0;i<N;i++)  
16          {  
17              for(j=0;j<N;j++)                  
18              {  
19                  if(i+j<N)//超出500位的部分不计算
20                  {
21                    c[i+j]+=ans[j]*anspow[i]; 
22                  }                              
23              }  
24             for(j=0;j<N-1;j++)  
25             {  
26                  if(c[j]>=10000) 
27                  {  
28                      c[j+1]+=c[j]/10000;//压4位  
29                       c[j]%=10000;              
30                  }                    
31             }  
32          }  
33          memcpy(ans,c,N*sizeof(int));  //复制函数
34  }
35  int main()
36  {
37    int P,i,j;
38    while(cin>>P)
39    {
40      memset(ans,0,sizeof(ans));
41      memset(anspow,0,sizeof(anspow));
42      printf("%d\n",(int)(P*log10(2)+1));
43      ans[0]=1;
44      anspow[0]=2;
45  /************关键部分计算2^P*******
46   2^p=(2^1)*(2^2)*(2^3)*(2^4)*(2^5)………………
47   简单说下:P=5 -----101(二进制)
48      p & 1 =1(最右边一位) -->>>ans=2 anspow=2
49      p>>1=110  anspow=2^2
50      p & 1 =0  此时表明2^3不存在 ans=2 anspow=2^4
51      p>>1=1       
52      p & 1 =1  ------>>>>> ans=2^5 
53      p>>1=0   ------结束    
54  ************************************/
55      while(P)
56      {
57          if( P & 1)
58              mult(ans,anspow);
59          P>>=1;
60              mult(anspow,anspow);
61      }
62      ans[0]--;//2^P的个位为2,4,6,8,故可以-1
63  /****************输出格式的控制************************/
64      for(i=124;i>=0;i--)
65      {
66          if(i%25==12)  
67          {
68             printf("%02d\n%02d",ans[i]/100,ans[i]%100);  
69          }
70          else  
71          {  
72             printf("%04d",ans[i]);
73             if(i%25==0) printf("\n");
74          }     
75      }
76  /***************************************************/
77    }
78  return 0;
79  }

 

 1 /*************麦森数****************/
 2  #include<iostream>
 3  #include<cstdio>
 4  #include<cmath>
 5  #define N 100   //压5位
 6  using namespace std;
 7  int ans[N];
 8  void  mult(int t)
 9  {
10    int i,temp,last=0;
11        for ( i=N-1; i>=0; i--)
12      { 
13          temp=(ans[i]<<t)+last;   //乘2^t,加进位 
14          last=temp/100000; 
15          ans[i]=temp-last*100000;  //temp%100000 
16      } 
17  }
18  void output()
19  {
20     int i;
21     for(i=1;i<=N;i++)
22     {     printf("%05d",ans[i-1]);
23          if (i%10==0)cout<<endl;   
24     } 
25  }
26  int main()
27  {
28    int P,times;
29    while(cin>>P)
30    {
31        memset(ans,0,sizeof(ans));
32        ans[N-1]=1;
33        cout<<(int)(P*log10(2)+1)<<endl;
34  /***********************关键部分*****************
35    2^P=2^(14*times)*2^(P%14) 用移位
36    之所以取14的原因 2^14*(99999)=1.638382*10^9 在(int)
37    范围,15以后都会超int,主要体现在mult()中      
38   **********************************************/
39        times=P/14;  //只能到14 15以后压缩时会超范围
40        while(times--)
41            mult(14);
42         mult(P%14);
43         --ans[99];
44        output();
45    }
46  return 0;
47  }

 

转载于:https://www.cnblogs.com/hduacm/archive/2012/08/22/2650796.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/474981.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

LeetCode 616. 给字符串添加加粗标签(Trie树)

文章目录1. 题目2. 解题1. 题目 给一个字符串 s 和一个字符串列表 dict &#xff0c;你需要将在字符串列表中出现过的 s 的子串添加加粗闭合标签 <b> 和 </b> 。 如果两个子串有重叠部分&#xff0c;你需要把它们一起用一个闭合标签包围起来。 同理&#xff0c;如…

var与dynamic区别

注意&#xff1a;var与dynamic这两个关键字&#xff0c;只是看起来很相 似&#xff0c;仅此而已&#xff01; var表示“变量的类型是在编译时决定的”&#xff0c;但是dynamic表 示“变量的类型是在运行时决定的”。因此&#xff0c;dynamic与var具有截然不同的含义。 var让 …

LeetCode 158. 用 Read4 读取 N 个字符 II

文章目录1. 题目2. 解题1. 题目 给你一个文件&#xff0c;并且该文件只能通过给定的 read4 方法来读取&#xff0c;请实现一个方法使其能够读取 n 个字符。 注意&#xff1a;你的 read 方法可能会被调用多次。 read4 的定义&#xff1a; 参数类型: char[] buf返回类型: int …

小白学数据分析-----数据指标 累计用户数的使用

小白学数据分析--数据指标累计用户数的使用 累计用户数是指注册用户数的累计&#xff0c;即可以认为是新用户的累计。在一般的数据统计中&#xff0c;我们基本上都会涉及到这个指标&#xff0c;且这个指标是逐渐累加的&#xff0c;比如&#xff1a; 时间 注册…

LeetCode 751. IP 到 CIDR(贪心)

文章目录1. 题目2. 解题1. 题目 给定一个起始 IP 地址 ip 和一个我们需要包含的 IP 的数量 n&#xff0c;返回用列表&#xff08;最小可能的长度&#xff09;表示的 CIDR块的范围。 CIDR 块是包含 IP 的字符串&#xff0c;后接斜杠和固定长度。例如&#xff1a;“123.45.67.8…

ORA-23616:执行块5失败

解决方案&#xff1a; 1. 查看错误信息 STRMADMning>select error_number,error_message 2 from dba_recoverable_script_errors 3 where script_idD74179203F11445D8F3F3F77C0749A1D 4 and block_num5; 2. 根据错误信息进行排错 3. 根据实际情况进行执行/回滚/清除…

LeetCode 308. 二维区域和检索 - 可变(前缀和)

文章目录1. 题目2. 解题1. 题目 给你一个 2D 矩阵 matrix&#xff0c;请计算出从左上角 (row1, col1) 到右下角 (row2, col2) 组成的矩形中所有元素的和。 上述粉色矩形框内的&#xff0c;该矩形由左上角 (row1, col1) (2, 1) 和右下角 (row2, col2) (4, 3) 确定。其中&am…

实战 SQL Server 2008 数据库误删除数据的恢复 (转)

今天有个朋友很着急地打电话给我&#xff0c;他用delete语句误删除了SQL Server 2008数据库中两个表中的所有记录&#xff0c;而这个数据库之前没有任何备份。让我帮他解决一下&#xff0c;不然他要赔偿客户很多钱。 SQL Server中误删除数据的恢复本来不是件难事&#xff0c;从…

在线销售数据分析–人货场三维分析角度

文章目录一、数据来源及理解二、分析思路三、数据处理数据预处理数据清洗数据转换四、数据描述性统计五、三维分析-人用户质量分析用户类别分析DM(管理者)排名分析六、三维分析-货销售金额及销量分布情况商品退货率七、三维分析-场城市区域八、总结一、数据来源及理解 此次分析…

LeetCode 348. 判定井字棋胜负(计数)

文章目录1. 题目2. 解题1. 题目 请在 n n 的棋盘上&#xff0c;实现一个判定井字棋&#xff08;Tic-Tac-Toe&#xff09;胜负的神器&#xff0c;判断每一次玩家落子后&#xff0c;是否有胜出的玩家。 在这个井字棋游戏中&#xff0c;会有 2 名玩家&#xff0c;他们将轮流在棋…

C中几组指针

具体关于c指针说明可参考前面两篇文章。C中指针详解和C中复杂类型声明。 1、二维数组 下面就三种二维数组进行说明。 1: int **Ptr; 2: int *Ptr[ 5 ]; 3: int ( *Ptr )[ 5 ]; 以上三例都是整数的二维数组&#xff0c;都可以用形如 Ptr[ 1 ][ 1 ] 的方式访问其内容&#xff1b;…

LeetCode 694. 不同岛屿的数量(BFS/DFS+set)

文章目录1. 题目2. 解题2.1 BFS2.2 DFS1. 题目 给定一个非空01二维数组表示的网格&#xff0c;一个岛屿由四连通&#xff08;上、下、左、右四个方向&#xff09;的 1 组成&#xff0c;你可以认为网格的四周被海水包围。 请你计算这个网格中共有多少个形状不同的岛屿。 两个岛…

xcode 4.2开发 ——navigation controller 添加按钮

添加一个按钮 UIBarButtonItem *rightButton [[UIBarButtonItem alloc] initWithTitle:"Add" style:UITabBarSystemItemContacts target:self action:selector(clickSettings:)]; self.navigationItem.rightBarButtonItem rightButton; 新版本xcode不用release 添…

数据分析-书籍整理(二)

业务书籍 《数据化管理-洞悉零售及电子商务》讲解了关于零售和电商的一些知识&#xff0c;有很多实用案例。很有借鉴意义。 《游戏数据分析实战》游戏各个阶段的数据分析&#xff0c;方法&#xff0c;数据来源&#xff0c;案例等。 《增长黑客》这本书我看了两遍&#xff0c…

LeetCode 1516. Move Sub-Tree of N-Ary Tree(DFS)

文章目录1. 题目2. 解题1. 题目 Given the root of an N-ary tree of unique values, and two nodes of the tree p and q. You should move the subtree of the node p to become a direct child of node q. If p is already a direct child of q, don’t change anything. …

空气质量指数(AQI)分析与预测

文章目录前言一、背景Ⅰ数据来源Ⅱ 数据背景Ⅲ 分析目的二、数据探索性分析Ⅰ 数据类型Ⅱ 描述性统计Ⅲ 数据预处理1.缺失值处理1&#xff09;删除缺失值2&#xff09;均值/中值填充2.异常值处理3.重复值处理三、数据分析Ⅰ空气质量排名Ⅱ 全国空气质量1.全国空气质量等级统计2…

如果程序跑着跑着就崩溃了,查看内存

程序跑的时间长就发生中断。 启动任务管理器&#xff0c;选择“进程”选项卡&#xff0c;查看程序的内存。 发现程序运行的过程中内存不断增长。 那么原因可能是运行过程中不断地分配了新的内存而没有释放。 查找的方法是注释代码&#xff0c;然后看内存是不是还在增长&#xf…

LeetCode MySQL解题目录

已完成的 LeetCode MySQL 数据库题目。点击查看我的 LeetCode 算法解题目录。 已解决 123/123 - 简单 54 中等 51 困难 18 前置入门学习 MySQL 基本查询、条件查询、投影查询MySQL 排序、分页查询、聚合查询MySQL 多表查询、连接查询&#xff08;内连接、外连接&#xff09;…

Anaconda3使用过程中遇到的问题

文章目录数据可视化中文、负号显示Anaconda使用pyecharts安装方法方法一&#xff1a;pip install pyecharts安装方法二&#xff1a;清华镜像安装方法三&#xff1a;离线安装成功安装后遇到的问题成功安装却不能使用不能使用某一个图&#xff1a;bar\grid决策树sklearn.exceptio…