#include "calc.h"
#include <assert.h>

#define MODULE 2005

int domul (int a, int b)
{
	int ret;
	a %= 2005; b %= 2005;
	ret = ((a*b )% 2005 + 2005) % 2005;
	return ret;
}

static int inverse[MODULE];

static int ggt (int a, int b)
{
	while (a && b) {
		if (a > b) {
			a %= b;
		} else {
			b %= a;
		}
	}
	return a+b;
}


static void init_inverse ()
{
	int i, j;
	static int inverse_done = 0;
	if (inverse_done)
		return;
	inverse_done = 1;
	for (i=0; i<MODULE; ++i) {
		if (inverse[i])
			continue;
		if (ggt (i, 2005) > 1)
			continue;
		for (j=0; j<MODULE; ++j) {
			if (domul (i, j) == 1) {
				inverse[i] = j;
				inverse[j] = i;
			}
		}
	}
}

int dodiv (int a, int b)
{
	init_inverse ();
	b %= MODULE;
	b += MODULE;
	b %= MODULE;
	assert (inverse[b]);
	return domul (a, inverse[b]);
}

int doadd (int a, int b)
{
	int ret;
	a %= MODULE; b %= MODULE;
	ret = a+b+2* MODULE;
	return ret % MODULE;
}

int dosub (int a, int b)
{
	int ret;
	a %= MODULE;
	b %= MODULE;
	ret = a - b + 2 * MODULE;
	return ret % MODULE;
}