Subversion Repositories mod_sim_exp

/*
 ============================================================================
 Name        : ModExp.c
 Author      : 
 Version     :
 Copyright   : Your copyright notice
 Description : Hello World in C, Ansi-style
 ============================================================================
 */
 
#include <stdio.h>
#include <stdlib.h>
#include <gmp.h>
 
int main(void) {
	unsigned int base_width, exp_width;
	mpz_t m, g0, g1, e0, e1, R, R2, result, tmp, gt0, gt1, gt01;
	gmp_randstate_t state;
 
	gmp_randinit_default(state);
	mpz_init(m);
	mpz_init(g0);
	mpz_init(g1);
	mpz_init(e0);
	mpz_init(e1);
	mpz_init(R);
	mpz_init(R2);
	mpz_init(result);
	mpz_init(tmp);
	mpz_init(gt0);
	mpz_init(gt1);
	mpz_init(gt01);
 
	printf("-- input generator program\n");
	printf("--  generates test values per bit input pair\n");
 
	while (1){
 
		//read in base_width
		scanf("%d", &base_width);
		if (base_width == 0) break;
		scanf("%d", &exp_width);
 
		//generate modulus (must be uneven)
		mpz_urandomb(m, state, base_width);
		mpz_setbit(m, 0); //uneven
 
		//generate g0
		mpz_urandomb(g0, state, base_width);
 
		//generate g1
		mpz_urandomb(g1, state, base_width);
 
		//generate e0
		mpz_urandomb(e0, state, exp_width);
 
		//generate e1
		mpz_urandomb(e1, state, exp_width);
 
		//calculate R
		mpz_set_ui(R, 2);
		mpz_powm_ui(R, R, base_width, m); //R = 2^n mod m
 
		//calculate R2
		mpz_set_ui(R2, 2);
		mpz_powm(R2, R, R2, m); //R2 = Rē mod m = 2^2n mod m
 
		//calc precompute values
		mpz_mul(gt0, g0, R);
		mpz_powm_ui(gt0, gt0, 1, m);
		mpz_mul(gt1, g1, R);
		mpz_powm_ui(gt1, gt1, 1, m);
		mpz_mul(gt01, g0, g1);
		mpz_mul(gt01, gt01, R);
		mpz_powm_ui(gt01, gt01, 1, m);
 
		//calculate result
		mpz_powm(result, g0, e0, m);
		mpz_powm(tmp, g1, e1, m);
		mpz_mul(result, result, tmp);
		mpz_powm_ui(result, result, 1, m);
 
		printf("-- base_width, exp_width, g0, g1, e0, e1, m, R^2, result\n");
		printf("%d\n", base_width);
		printf("%d\n", exp_width);
		gmp_printf("%Zx\n", g0);
		gmp_printf("%Zx\n", g1);
		gmp_printf("%Zx\n", e0);
		gmp_printf("%Zx\n", e1);
		gmp_printf("%Zx\n", m);
		gmp_printf("%Zx\n", R2);
		gmp_printf("%Zx\n", R);
		gmp_printf("%Zx\n", gt0);
		gmp_printf("%Zx\n", gt1);
		gmp_printf("%Zx\n", gt01);
		gmp_printf("%Zx\n", result);
	}
 
	mpz_clear(g0);
	mpz_clear(g1);
	mpz_clear(e0);
	mpz_clear(e1);
	mpz_clear(m);
	mpz_clear(result);
	mpz_clear(tmp);
	mpz_clear(R);
	mpz_clear(R2);
	mpz_clear(gt0);
	mpz_clear(gt1);
	mpz_clear(gt01);
 
	return EXIT_SUCCESS;
}