代码改变世界

求数n的最大素因子是第几个素数

2012-12-27 14:30  coodoing  阅读(1772)  评论(0编辑  收藏  举报

问题描述

每个数都可以分解成质因数的乘积。现在求一个数n的最大素因子是素数集合中的第几个素数。n小于1000000。

解决方法

利用筛选法获得,具体可见求质数算法的N种境界- 试除法和初级筛法。 图解说明如下:

OQAAAMQ7u67PV8k_Vp7UpVJxenDnJ_h7kdt6BKZFbh75EuvIKjjoEVYbbESQ4Kr72tv2m1H3X9ktmzHsLQyI-XbpxC8A15jOjL4ByeTggK6hVvMCrhy5D6R0ydSL

1、具体过程

对于1000000内的数字而言,利用筛选法计算其中素数,并标记出每个素数对应的位置信息。利用prime数组保存范围内的素数;利用position数组保存每个数所对应的位置。

   1: public static final int NUM = 1000001; // 100万个数
   2: static int primeNum = 80000; // 根据素数定理估算得到
   3: static int prime[] = new int[primeNum]; // 素数
   4: static int position[] = new int[NUM]; // 每个素数对应的位置

2、代码

   1: public void init() {
   2:     int pos = 0;
   3:     for (int i = 2; i < NUM; i++) {
   4:         if (position[i] == 0) {
   5:             pos++;
   6:             prime[pos] = i; // 加入素数表
   7:             position[i] = pos; // 素数对应的位置:序数
   8:             for (int j = 2 * i; j < NUM; j += i)
   9:                 position[j] = -1;
  10:         }
  11:     }
  12: }
  13:  
  14: public int resolve(int n) {
  15:     int pos = 0;
  16:     for (int i = 1; prime[i] * prime[i] <= n; i++) {
  17:         while (n % prime[i] == 0) {
  18:             pos = i;
  19:             n /= prime[i];
  20:         }
  21:     }
  22:     if (n != 1) // n本身就是素数
  23:     {
  24:         pos = position[n];
  25:     }
  26:     return pos;
  27: }

参考:

1、prime number theorem

2、Euler's totient function

3、埃拉托斯特尼筛法