经典算法（1）- 求1到N的质数

用Eratosthenes筛子法：

import java.util.ArrayList;
import java.util.List;

/**
 * 
 * @author ljs
 * Sieve of Eratosthenes: determine the prime numbers within the range 1....n (n>0)
 *
 */
public class Sieve {
	//intput: n is an positive integer (n>0)
	//output: a list of primes from 1 to n, as 2,3,5.... 
	public static List<Integer> solve(int n){
		int[] A = new int[n+1]; 
		//prepare
		for(int i=2;i<=n;i++){
			A[i]=i;
		}
		
		//sieve process
		int p = (int)(Math.sqrt(n));
		for(int i=2;i<=p;i++){
			//check if i is eliminated by previous steps
			if(A[i] != 0){
				int j = i*i;
				while(j<=n){
					//eliminated the sequence:  i, i*(i+1)....
					A[j]=0;
					j += i;
				}			
			}
		}
		
		//collect the result into list
		List<Integer> primes=new ArrayList<Integer>();
		for(int i=2;i<=n;i++){
			if(A[i] != 0)
			primes.add(A[i]);
		}
		return primes;
	}

	public static void main(String[] args) {
		int n = 1000;
		List<Integer> primes=Sieve.solve(n);
		System.out.format("total number of primes from 1 to %d is: %d%n", n, primes.size());
		int i = 0;
		for(int p:primes){
			if(i % 10 == 0)
				System.out.println();
			System.out.format(" %4d",p);
			i++;
		}
		System.out.println();
	}
}

测试：

total number of primes from 1 to 1000 is: 168

    2    3    5    7   11   13   17   19   23   29
   31   37   41   43   47   53   59   61   67   71
   73   79   83   89   97  101  103  107  109  113
  127  131  137  139  149  151  157  163  167  173
  179  181  191  193  197  199  211  223  227  229
  233  239  241  251  257  263  269  271  277  281
  283  293  307  311  313  317  331  337  347  349
  353  359  367  373  379  383  389  397  401  409
  419  421  431  433  439  443  449  457  461  463
  467  479  487  491  499  503  509  521  523  541
  547  557  563  569  571  577  587  593  599  601
  607  613  617  619  631  641  643  647  653  659
  661  673  677  683  691  701  709  719  727  733
  739  743  751  757  761  769  773  787  797  809
  811  821  823  827  829  839  853  857  859  863
  877  881  883  887  907  911  919  929  937  941
  947  953  967  971  977  983  991  997

转载于:https://www.cnblogs.com/ljsspace/archive/2011/06/05/2073359.html