//testuj odwrotnosc modulo void Tester::autotest_inverse(int len){ Large one = Large::Set("1", base, in_base); Large zero = Large::Set("0", base, in_base); TesterRational tr = TesterRational(); int i; Large a; Large mod; Large inv; srand( (unsigned)time(NULL) ); for (i = 0; i < TIMES; i++) { a = Large::Set(tr.randomString(len), base, in_base); mod = Large::Set(tr.randomString(len), base, in_base); //jeœli a jest wiêksza ni¿ mod to nie ma sensu szukaæ if(a > mod){ i--; continue; } inv = a.inverseMod(mod); Large result = Large(); (a * inv).divide(mod, result); //wynik jest b³êdny gdy reszta nie jest równa jeden i inv nie jest zerem(znaczy to ze nie ma odwrotnosci) if (result.compareAbsolute(one)!=0 && inv.compareAbsolute(zero)!=0) { cout << " auto_Inv:Error " << endl; cout << " a " << a.toString() << " inv " << inv.toString() << " mod " << mod.toString() << endl;; } cout << "." ; } cout << "k" << endl; }
bool process_number(int * pan, Large a) { for (int i = 0; i < a.size(); i++) { if (pan[a.get(i)] == 1) { return false; } else { pan[a.get(i)] = 1; } } return true; }
bool is_pandigital(Large a, Large b, Large c) { if ((a.toString()+b.toString()+c.toString()).length() != 9) { return false; } int pan [10]; int i; for (i = 0; i < 10; i++) { pan[i] = 0; } return process_number(pan,a) && process_number(pan,b) && process_number(pan,c) && pan[0] == 0; }
//testuj odwrotnosc modulo void Tester::test_inverse(){ //test1 Large one = Large::Set("1", base, in_base); Large zero = Large::Set("0", base, in_base); Large test1 = Large::Set("12", base, in_base); Large mod = Large::Set("2f", base, in_base); Large inv = test1.inverseMod(mod); Large result; //jeœli istnieje liczba odwrotna(nie jest zerem) i odwrotnoœæ razy liczba mod to jeden, // to wynik jest poprawny (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse1: " << test1.toString() << " " << mod.toString() << endl; } //test2 test1 = Large::Set("3", base, in_base); mod = Large::Set("7", base, in_base); inv = test1.inverseMod(mod); (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse2: " << test1.toString() << " " << mod.toString() << endl; } //test3 test1 = Large::Set("c", base, in_base); mod = Large::Set("11", base, in_base); inv = test1.inverseMod(mod); (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse3: " << test1.toString() << " " << mod.toString() << endl; } //test4 test1 = Large::Set("21", base, in_base); mod = Large::Set("61", base, in_base); inv = test1.inverseMod(mod); (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse4: " << test1.toString() << " " << mod.toString() << endl; } //test5 test1 = Large::Set("b", base, in_base); mod = Large::Set("17", base, in_base); inv = test1.inverseMod(mod); (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse5: " << test1.toString() << " " << mod.toString() << endl; } //test6 test1 = Large::Set("5", base, in_base); mod = Large::Set("13", base, in_base); inv = test1.inverseMod(mod); (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse6: " << test1.toString() << " " << mod.toString() << endl; } //test7 test1 = Large::Set("80722", base, in_base); mod = Large::Set("41413", base, in_base); inv = test1.inverseMod(mod); (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse7: " << test1.toString() << " " << mod.toString() << endl; } //test8 test1 = Large::Set("39782", base, in_base); mod = Large::Set("43d35", base, in_base); inv = test1.inverseMod(mod); Large test = test1.inverseMod(mod); (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse8: " << test1.toString() << " " << mod.toString() << endl; } //test9 test1 = Large::Set("45881", base, in_base); mod = Large::Set("152", base, in_base); inv = test1.inverseMod(mod); test = test1.inverseMod(mod); (test1 * inv).divide(mod, result); if (inv.compareAbsolute(zero)!=0 && result.compareAbsolute(one)!=0){ cout << "Test_Inverse9: " << test1.toString() << " " << mod.toString() << endl; } }
void Tester::autotest_crt(int len){ TesterRational tr = TesterRational(); Large a; Large a_m; Large b; Large b_m; Large c; Large c_m; vector<Large> vect = vector<Large>(); srand( (unsigned)time(NULL) ); Large result(base); for(int i=0; i<TIMES; i++){ a = Large::Set(tr.randomString(len/2), base, in_base); a_m = Large::Set(tr.randomString(len), base, in_base); a.divide(a_m, a); b = Large::Set(tr.randomString(len/2), base, in_base); b_m = Large::Set(tr.randomString(len), base, in_base); while(!Large::coPrime(a_m, b_m)){ b_m = Large::Set(tr.randomString(len), base, in_base); } b.divide(b_m, b); c = Large::Set(tr.randomString(len/2), base, in_base); c_m = Large::Set(tr.randomString(len), base, in_base); while(!Large::coPrime(b_m, c_m) || !Large::coPrime(a_m, c_m)){ c_m = Large::Set(tr.randomString(len), base, in_base); } c.divide(c_m, c); vect.clear(); vect.push_back(a); vect.push_back(a_m); vect.push_back(b); vect.push_back(b_m); vect.push_back(c); vect.push_back(c_m); result = Large::crt(vect); if(result.mod(a_m) != a){ cout << "AutoTest CRT: Error przy a= " << a.toString() << " mod " << a_m.toString() << " wynik " << result.toString() << endl; cout << "AutoTest CRT: Error przy b= " << b.toString() << " mod " << b_m.toString() << " wynik " << result.toString() << endl; cout << "AutoTest CRT: Error przy c= " << c.toString() << " mod " << c_m.toString() << " wynik " << result.toString() << endl; break; } if(result.mod(b_m) != b){ cout << "AutoTest CRT: Error przy a= " << a.toString() << " mod " << a_m.toString() << " wynik " << result.toString() << endl; cout << "AutoTest CRT: Error przy b= " << b.toString() << " mod " << b_m.toString() << " wynik " << result.toString() << endl; cout << "AutoTest CRT: Error przy c= " << c.toString() << " mod " << c_m.toString() << " wynik " << result.toString() << endl; break; } if(result.mod(c_m) != c){ cout << "AutoTest CRT: Error przy a= " << a.toString() << " mod " << a_m.toString() << " wynik " << result.toString() << endl; cout << "AutoTest CRT: Error przy b= " << b.toString() << " mod " << b_m.toString() << " wynik " << result.toString() << endl; cout << "AutoTest CRT: Error przy c= " << c.toString() << " mod " << c_m.toString() << " wynik " << result.toString() << endl; break; } cout << "." ; } }