密码学—椭圆曲线-求P的周期

#include<stdio.h>
int Ni(int m,int p)//模p下m的乘法逆元
{
 int n,i;
 if(m==1)
 {
  n=1;
  return n;
 }
 while(m<0) m+=p; 
 for( i=1;i<p;i++)
 {
  if((m*i)%p==1)
  {
   n=i;
   return n;
  }
 } 
}

int Lmn(int x1,int y1,int x2,int y2,int p)//模p下
{
 int result;
 if(x1==x2&&y1==y2)
  result=(3*x1*x1+1)*Ni(2*y1,p);
 else
  result=(y2-y1)*Ni(x2-x1,p);
 while(result<0)
  result+=p;