]> www.ginac.de Git - cln.git/blob - examples/perfnum.cc
Update to currently largest known Mersenne prime.
[cln.git] / examples / perfnum.cc
1 // This program prints the largest now known even perfect number.
2
3 #include <cln/integer.h>
4 #include <cln/integer_io.h>
5
6 using namespace std;
7 using namespace cln;
8
9 int main ()
10 {
11         // previous ones were 1257787, 1398269, 2976221, 3021377, 6972593,
12         // 13466917, 20996011, 24036583, 25964951, 30402457, 32582657, 37156667,
13         // 42643801, 43112609, 57885161, 74207281, 77232917, 82589933
14         int p = 82589933;
15         cl_I x = (((cl_I)1 << p) - 1) << (p-1);
16         cout << x << endl;
17 }