| 2 |
- |
1 |
#include "os.h"
|
|
|
2 |
#include <mp.h>
|
|
|
3 |
#include <libsec.h>
|
|
|
4 |
|
|
|
5 |
void
|
|
|
6 |
main(void)
|
|
|
7 |
{
|
|
|
8 |
mpint *z = mpnew(0);
|
|
|
9 |
mpint *p = mpnew(0);
|
|
|
10 |
mpint *q = mpnew(0);
|
|
|
11 |
mpint *nine = mpnew(0);
|
|
|
12 |
|
| 33 |
7u83 |
13 |
fmtinstall('B', mpfmt);
|
| 2 |
- |
14 |
strtomp("2492491", nil, 16, z); // 38347921 = x*y = (2**28-9)/7,
|
|
|
15 |
// an example of 3**(n-1)=1 mod n
|
|
|
16 |
strtomp("15662C00E811", nil, 16, p);// 23528569104401, a prime
|
|
|
17 |
uitomp(9, nine);
|
|
|
18 |
|
|
|
19 |
if(probably_prime(z, 5) == 1)
|
|
|
20 |
fprint(2, "tricked primality test\n");
|
|
|
21 |
if(probably_prime(nine, 5) == 1)
|
|
|
22 |
fprint(2, "9 passed primality test!\n");
|
|
|
23 |
if(probably_prime(p, 25) == 1)
|
|
|
24 |
fprint(2, "ok\n");
|
|
|
25 |
|
|
|
26 |
DSAprimes(q, p, nil);
|
|
|
27 |
print("q=%B\np=%B\n", q, p);
|
|
|
28 |
|
|
|
29 |
exits(0);
|
|
|
30 |
}
|
|
|
31 |
|
|
|
32 |
// example output, checked with Maple:
|
|
|
33 |
// seed EB7B6E35F7CD37B511D96C67D6688CC4DD440E1E
|
|
|
34 |
// q=E0F0EF284E10796C5A2A511E94748BA03C795C13
|
|
|
35 |
// = 1284186945063585093695748280224501481698995297299
|
|
|
36 |
// p=C41CFBE4D4846F67A3DF7DE9921A49D3B42DC33728427AB159CEC8CBBDB12B5F0C244F1A734AEB9840804EA3C25036AD1B61AFF3ABBC247CD4B384224567A863A6F020E7EE9795554BCD08ABAD7321AF27E1E92E3DB1C6E7E94FAAE590AE9C48F96D93D178E809401ABE8A534A1EC44359733475A36A70C7B425125062B1142D
|
|
|
37 |
// = 137715385439333164327584575331308277462546592976152006175830654712456008630139443747529133857837818585400418619916530061955288983751958831927807888408309879880101870216437711393638413509484569804814373511469405934988856674935304074081350525593807908358867354528898618574659752879015380013845760006721861915693
|
|
|
38 |
// r=DF310F4E54A5FEC5D86D3E14863921E834113E060F90052AD332B3241CEF2497EFA0303D6344F7C819691A0F9C4A773815AF8EAECFB7EC1D98F039F17A32A7E887D97251A927D093F44A55577F4D70444AEBD06B9B45695EC23962B175F266895C67D21C4656848614D888A4
|
|
|
39 |
// = 107239359478548771267308764204625458348785444483302647285245969203446101233421655396874997253111222983406676955642093641709149748793954493558324738441197139556917622937892491175016280660608595599724194374948056515856812347094848443460715881455884639869144172708
|
|
|
40 |
// g=2F1C308DC46B9A44B52DF7DACCE1208CCEF72F69C743ADD4D2327173444ED6E65E074694246E07F9FD4AE26E0FDDD9F54F813C40CB9BCD4338EA6F242AB94CD410E676C290368A16B1A3594877437E516C53A6EEE5493A038A017E955E218E7819734E3E2A6E0BAE08B14258F8C03CC1B30E0DDADFCF7CEDF0727684D3D255F1
|
|
|
41 |
// = 33081848392740465806285326014906437543653045153885419334085917570615301913274531387168723847139029827598735376746057461417880810924280288611116213062512408829164220104555543445909528701551198146080221790002337033997295756585193926863581671466708482411159477816144226847280417522524922667065714073338662508017
|