c语言简单实现elgamal椭圆曲线算法

#include<cstdio>
#include<cmath>
#include<iostream>
#include<time.h>




using namespace std;
struct zuobiao
{
int x;
int y;
};


int  EX_gcd(int a, int b, int &s,int &t)//无问题
{
if (a == 0 && b == 0) {
return -1;
}
if (b == 0) {
s = 1;
t = 0;
return a;
}
int  d = EX_gcd(b, a%b, t,s);
t -= a / b*s;
return d;
}
int  f1(zuobiao P,zuobiao q,int a,int p) //求λ的值  这里的除法a/b=a*b^-1（modp）
{ int  k,c,b,e,t;
if ((P.x==q.x)&&(P.y==q.y)) { //此处只有加法规则的第三步，无（1）（2）
c = ((3 * q.x*q.x + a)%p+p)%p;
b= (2 *P.y);
EX_gcd(b, p, e,t);
while (e < 0) {
e = (e + p)% p;
}
k = (c*e+p)%p;
//cout << k << endl;
}
else {
c = ((q.y - P.y) % p + p ) % p;
b = ((q.x - P.x) % p + p) % p;