《素数测试》


Solovag-StrassonRobert Solovag和Volker Strasson开发了一种概率的基本测试算法。
这个算法使用了雅可比函数来测试p是否为素数：
（1） 选择一个小于p的随机数a。
（2） 如果GCD（a,p）<>1，那么p通不过测试，它是合数。
（3） 计算j=a^(p-1)/2 mod p。

（4） 计算雅可比符号J（a,p）。
（5） 如果j<>J（a,p），那么p肯定不是素数。
（6） 如果j=J(a,p),那麽p不是素数的可能性值多是50%数a被称为一个证据，
如果a不

import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStreamReader;
import java.math.BigInteger;

public class bigTest {

public static void main(String[] args) throws IOException {
// TODO Auto-generated method stub
BufferedReader br = new BufferedReader(new InputStreamReader(System.in));
String num1 = "";// 第一个数
String num2 = "";// 第二个数

num1 = br.readLine();
while (

11G 王旁青头五夫一，
12F 土干十寸未甘雨，
13D 大犬戊其古石厂，
14S 木丁西甫一四里，
15A 工戈草头右框七。

21H 目上卜止虎头具，
22J 日早两竖与虫依，
23K 口中两川三个竖，
24L 田甲方框四车里，
25M 山由贝骨下框集。

31T 禾竹反文双人立，
32R 白斤气丘叉手提，
33E 月用力豸毛衣臼，
34W 人八登头单人几，
35Q 金夕鸟儿犭边鱼。

41Y 言文方点谁人去，
42U 立辛六羊病门里，
43I 水族三点鳖头小，
44O 火业广鹿四点米，
45

#include<stdio.h>
#include"miracl.h"
#define LIMIT 15000
int main()
{
int n,p;
big x,y;
miracl *mip=mirsys(5000,10);
x=mirvar(0);
y=mirvar(0);
gprime(LIMIT);
printf("input number to be factored\n");
cinnum(x,stdin);
if(isprime(x))
{
   printf("this number is prime!\n");
   return 0;
}
n=0;
p=mip->PRIMES[0];
printf("factors are \n

f[n_] := Sum[1/k, {k, 1, n}]
r[n_] := f[n] + Exp[f[n]]*Log[f[n]]
g[n_] := DivisorSigma[1, n]
h[n_] := GreaterEqual[N[f[n]], g[n]]

Select[Range[1, 10000], h]

在公式中插入文本

可以通过\mbox{text}在公式中添加text，
\documentclass{article}
\usepackage{CJK}
\begin{CJK*}{GBK}{song}
\begin{document}
$$\mbox{对任意的$x>0$}, \mbox{有 }f(x)>0. $$
\end{CJK*}
\end{document}