定积分旋转体体积公式用编程语言实现椭圆曲线群上的倍点公式

2020年10月11日发(作者：谢介子)
《信息安全概论》习题
学号： 姓名： 班级：

一、所使用的工具软件及环境：windows7+eclipse+jdk1.6
二、题目要求：
用编程语言实现椭圆曲线群上的倍点公式

三、实现代码：
import .*;
public class Test
{
public static int x,y,n,a,p;
public static int getX(int x)求逆
{
int i=1;
while(true)
{
if(modP(x*i)==1)
{
break;
}
i++;
}
return i;
}
public static int modP(int x)求模
{
if(x<0)
{
while(x<0)
{
x=x+p;
}
}
else
x=x%p;
return x;
}
public static void getR()计算倍点，得结果
{
int X=x,Y=y;
for(int i=1;i {
if(X==x&Y==y)
{
int k=getK(a,x,y);
X=modP(k*k-x-X);
Y=modP(k*(x-X)-y);
}
else
{
int k=getk(x,y,X,Y);
X=modP(k*k-x-X);
Y=modP(k*(x-X)-y);
}
}
n(+X++Y+);
}
public static int getK(int a,int x,int y)计算倍点K
{
int k=(3*x*x+a)*getX(2*y);
return modP(k);
}
public static int getk(int x1,int y1,int x2,int y2)计算点加K
{
int k=(y2-y1)*getX(x2-x1);
return modP(k);
}
public static void main(String args[]) throws IOException
{
BufferedReader f=new BufferedReader(new
InputStreamReader());
n(请按如下格式输入要计算的倍点：);
n(椭圆曲线方程：y^2=aX^3+bx+c,请输入：a);
String s=ne();
a=nt(s);
n(所求点的坐标：Q（x，y）,请输入x，y);
s=ne();
x=nt(s);
s=ne();
y=nt(s);
n(请输入素数p：);
s=ne();
p=nt(s);
n(所求的倍点:nQ,请输入n);
s=ne();
n=nt(s);
(点Q（+x++y+的+n+倍点结果+n+是);
getR();
}

}

四、程序运行结果：