/****************************************************************************/
/* Busca numeros primos gemelos a partir de un número dado.Usa LIP.         */
/*                                                                          */
/* Jaime Suarez <mcripto@bigfoot.com> 2003                                  */
/* en http://elparaiso.mat.uned.es                                            */
/****************************************************************************/

#include "lip.h"
#include <stdio.h>

main(int argc, char *argv[])
{
	verylong a=0,b=0,tmp=0;

	if (argc!=2) {
		printf("%s <n> \n",argv[0]);
		printf("busca un par de primos gemelos mayores que n.\n");
		return 1;
	}
	zsread(argv[1],&a);
	
	/* Si no es impar le sumamos uno para empezar */
	if (!zodd(a)) { zsadd(a,1,&tmp); zcopy(tmp,&a);}
	
	while(1) {
		zsadd(a,2,&b);                                     /* b=a+2 */
		if (zprobprime(a,100) && zprobprime(b,100)) break;
		zcopy(b,&a);                                       /* a=a+2 */
	}
		
	zwrite(a);printf(" y ");zwrite(b);
	printf(" son primos gemelos.\n");
	return 0;
}