libbrandt/gp-scripts/smc.gp

/* search for a random prime with the specified amount of bits */
\\ randomprime(bits:small=127)=
\\ {
\\ 	local(r:int=4);
\\ 	while(!isprime(r),
\\ 		r = bitor(2^(bits-1) + random(2^(bits-1)-1), 1);
\\ 	);
\\ 	r;
\\ }

smc_hextodec(s:str) =
{
	local(v:vecsmall = Vecsmall(s), ret:int = 0);
	for(i = 1, #v,
		ret = (ret<<4) + if(v[i]<=57 && v[i]>=48, v[i]-48, v[i]<=70 && v[i]>=65, v[i]-55, v[i]<=102 && v[i]>=97, v[i]-87, error("invalid input format"))
	);
	ret;
}

smc_genbid(k:small, bid:small, g)=
{
	vector(k,j,g^(bid==j));
}

smc_genalpha(k:small, b:vec, r:vec, y)=
{
	vector(k, j, b[j]*y^r[j]);
}

smc_genbeta(k:small, r:vec, g)=
{
	vector(k, j, g^r[j]);
}