#包括#包括#包括int层次=3/*腿模基础测试*/int Grenze=1000/*Legt den Bereich节日,fuer den a(k)berechnet wird*/mpz_t ErgFeld[1000];文件*Ausgabe;int测试(mpz_t k){int碑文;mpz_t基础;mpz_init_set_ui(基础,层次);mpz_t Wert;mpz _ init(Wert);mpz_powm(Wert,基础,k,k);如果(mpz_cmp_si(Wert,Grenze)<0){Stelle=mpz_get_ui(Wert);if(mpz_cmp_si(ErgFeld[Stelle],-1)==0){mpz_set(ErgFeld[Stelle],k);mpz_clear(基础);mpz_clear(Wert);gmp_fprintf(Ausgabe,“a(%d)=%Zd\n”,钨铬钴合金,k);fflush(奥斯加贝);返回(1);}}mpz_clear(基础);mpz_clear(Wert);返回(0);}整型main(){整数i;mpz_t扎勒;int Fehler=格伦泽;mpz_init_set_ui(Zaehler,2);Ausgabe=fopen(“3Erg.txt”,“w”);if(Ausgabe==空){printf(“Dateifehler!\n”);}其他的{对于(i=0;i2){费勒-=试验(扎勒);mpz_add_ui(扎勒,扎勒,1);if(mpz_fdiv_ui(Zaehler,1000000)==0){gmp_printf(“Jetzt:%Zd\n”,扎勒);printf(“Es fehlen noch%d.\n”,Fehler);}}}mpz_clear(扎勒);fclose(奥斯加贝);返回(0);}