Oplossingen COMA 2012

4183 days ago by Giles.Miclotte

7 is leuk, 23 niet!

Antwoord: 68323018645

We kunnen alle getallen van 0 tot 1000000 overlopen en kijken welke er aan de voorwaarden voldoen:

def coma1_0(l,nl,max): sum=0 for i in [0..max]: if l.divides(i): if not nl.divides(i): sum+=i return sum coma1_0(7,23,1000000) 
       
68323018645
68323018645

Alternatief. De gezochte getallen zijn wel deelbaar zijn door 7, maar niet door 23 (en telkens kleiner dan 1000000).  Dan is het duidelijk dat dit de getallen zijn die deelbaar zijn door 7, maar niet door $7\cdot 23$. De som van alle zevenvouden kunnen we berekenen met een rekenkundige rij, waarvan de eerste term 0 is en de laatste $i$, met $i$ het grootste zevenvoud kleiner dan 1000000. Analoog voor de $7\cdot 23$-vouden. Daaruit volgt:

def coma1_1(l,nl,max): m=floor(max/l) n=floor(m/nl) return l*(m*(m+1)/2-nl*n*(n+1)/2) coma1_1(7,23,1000000) 
       
68323018645
68323018645

Diophantus 2012

Antwoord: 5478676

We herschrijven de vergelijking als een kwadratische vergelijking in $x$:

$x^2 +x(2y-2012)+y^2+4024y=0$. Deze vergelijking heeft discriminant $D=4^2503^2\left(1-\frac{3y}{503}\right)$. $x$ en $y$ kunnen dus enkel geheel zijn als $y=-530k$, met $k\geq0$ geheel, en dan is $x=503(2+k\pm\sqrt{1+3k})$.

Definieer $s=\sqrt{1+3k}$, we vinden een bovengrens voor $s$ in functie van $n$ (=$300000$) door $503(2+k+2s)=n$ op te lossen, en hierna kunnen we de som berekenen.

def coma2_0(n, s=var('s')): bound=solve(503*(2+(s^2-1)/3+2*s)==n,s) r=max(abs(ceil(bound[0].rhs())),abs(ceil(bound[1].rhs()))) sum=0 for s in [1..r]: if mod(s^2,3)==1: k=(s^2-1)/3 y=-k*503 x1=503*(2+k+2*s) x2=503*(2+k-2*s) if abs(x1)<n: sum+=x1 if abs(x2)<n: sum+=x2 return sum coma2_0(300000) 
       
5478676
5478676

Alternatief:

Stel $a=x+y$ en $b=x-2y$. Dan geldt $x=\frac{2a+b}{3}$ en $y=\frac{a-b}{3}$ en wordt de vergelijking $a^2=2012b$. Merk op dat $(x,y)$ een gehele oplossing is als en slechts als $(a,b)$ een gehele oplossing is met $2a+b$ deelbaar door $3$. De gehele oplossingen van $a^2=2012b$ zijn de koppels $(a,\frac{a^2}{2012})$ met  $a$ en $\frac{a^2}{2012}$ geheel. Dit laatste is het geval als en slechts als a deelbaar is door $1006$, dus we kunnen stellen dat $a=1006c$, zodat we de oplossingen kunnen schrijven als $(1006c,503c^2)$. Samen met $x=\frac{2a+b}{3}$ vinden we dat de gezochte som de som is van alle $\frac{2012c+503c^2}{3}$, met $c$ geheel, $2012c+503c^2$ deelbaar door $3$, en $|\frac{2012c+503c^2}{3}|<300000$.

def coma2_1(n): antw=0 c=0 while(abs(2012*c+503*c^2)<3*n): if mod(2012*c+503*c^2,3)==0: antw+=(2012*c+503*c^2)/3 c+=1 c=-1 while(abs(2012*c+503*c^2)<3*n): if mod(2012*c+503*c^2,3)==0: antw+=(2012*c+503*c^2)/3 c-=1 print antw coma2_1(300000) 
       
5478676
5478676

Productje uitrekenen

Antwoord: 3896

Na wat proberen kunnen we vermoeden dat $f(n) \equiv n+1 \pmod{6037}$ voor $n \in \{1,...,3018\}$ en $f(3019) \equiv 1 \pmod{6037}$. We kunnen dit door sage laten controleren:

%hide f(x)=1326*x^3018 + 1675*x^3017 + 482*x^3016 + 1208*x^3015 + 5290*x^3014 + 4038*x^3013 + 4493*x^3012 + 710*x^3011 + 4701*x^3010 + 2050*x^3009 + 5158*x^3008 + 1645*x^3007 + 4334*x^3006 + 4080*x^3005 + 3718*x^3004 + 2104*x^3003 + 1759*x^3002 + 4014*x^3001 + 5385*x^3000 + 4242*x^2999 + 4121*x^2998 + 1826*x^2997 + 4065*x^2996 + 4861*x^2995 + 971*x^2994 + 2320*x^2993 + 4570*x^2992 + 589*x^2991 + 3088*x^2990 + 2943*x^2989 + 505*x^2988 + 5432*x^2987 + 394*x^2986 + 3420*x^2985 + 78*x^2984 + 3711*x^2983 + 1774*x^2982 + 2789*x^2981 + 456*x^2980 + 2821*x^2979 + 5053*x^2978 + 2603*x^2977 + 2259*x^2976 + 5653*x^2975 + 5166*x^2974 + 2467*x^2973 + 5045*x^2972 + 4802*x^2971 + 1631*x^2970 + 3814*x^2969 + 4128*x^2968 + 1407*x^2967 + 4412*x^2966 + 3614*x^2965 + 5691*x^2964 + 2076*x^2963 + 4273*x^2962 + 1719*x^2961 + 4774*x^2960 + 2259*x^2959 + 2775*x^2958 + 616*x^2957 + 433*x^2956 + 5295*x^2955 + 4265*x^2954 + 2350*x^2953 + 237*x^2952 + 2042*x^2951 + 4443*x^2950 + 1233*x^2949 + 1100*x^2948 + 2168*x^2947 + 853*x^2946 + 4051*x^2945 + 2378*x^2944 + 2012*x^2943 + 1663*x^2942 + 3105*x^2941 + 3*x^2940 + 633*x^2939 + 5297*x^2938 + 2434*x^2937 + 2011*x^2936 + 446*x^2935 + 1959*x^2934 + 144*x^2933 + 1634*x^2932 + 4504*x^2931 + 1933*x^2930 + 5478*x^2929 + 4640*x^2928 + 710*x^2927 + 2209*x^2926 + 4788*x^2925 + 2158*x^2924 + 4757*x^2923 + 2995*x^2922 + 4553*x^2921 + 1051*x^2920 + 1119*x^2919 + 5570*x^2918 + 2824*x^2917 + 6025*x^2916 + 4654*x^2915 + 2586*x^2914 + 52*x^2913 + 165*x^2912 + 5530*x^2911 + 2390*x^2910 + 685*x^2909 + 4131*x^2908 + 3664*x^2907 + 3144*x^2906 + 1159*x^2905 + 1666*x^2904 + 253*x^2903 + 3108*x^2902 + 988*x^2901 + 1559*x^2900 + 3534*x^2899 + 4983*x^2898 + 772*x^2897 + 1215*x^2896 + 4716*x^2895 + 5988*x^2894 + 3909*x^2893 + 3662*x^2892 + 2517*x^2891 + 886*x^2890 + 4261*x^2889 + 4488*x^2888 + 2198*x^2887 + 3393*x^2886 + 1888*x^2885 + 4530*x^2884 + 2440*x^2883 + 2096*x^2882 + 5226*x^2881 + 1960*x^2880 + 3035*x^2879 + 5838*x^2878 + 743*x^2877 + 1208*x^2876 + 4441*x^2875 + 4257*x^2874 + 3397*x^2873 + 4861*x^2872 + 3312*x^2871 + 56*x^2870 + 3305*x^2869 + 4041*x^2868 + 5220*x^2867 + 1711*x^2866 + 2286*x^2865 + 6000*x^2864 + 3757*x^2863 + 1859*x^2862 + 4754*x^2861 + 143*x^2860 + 5550*x^2859 + 2044*x^2858 + 732*x^2857 + 1215*x^2856 + 2210*x^2855 + 5891*x^2854 + 3313*x^2853 + 1426*x^2852 + 4943*x^2851 + 354*x^2850 + 4264*x^2849 + 962*x^2848 + 1565*x^2847 + 2577*x^2846 + 2494*x^2845 + 591*x^2844 + 1227*x^2843 + 2318*x^2842 + 4082*x^2841 + 4877*x^2840 + 4581*x^2839 + 5186*x^2838 + 4774*x^2837 + 1732*x^2836 + 330*x^2835 + 676*x^2834 + 1054*x^2833 + 3562*x^2832 + 2072*x^2831 + 1795*x^2830 + 2326*x^2829 + 62*x^2828 + 985*x^2827 + 1708*x^2826 + 1570*x^2825 + 764*x^2824 + 3738*x^2823 + 2772*x^2822 + 3000*x^2821 + 2586*x^2820 + 3394*x^2819 + 3322*x^2818 + 1791*x^2817 + 1763*x^2816 + 3310*x^2815 + 2610*x^2814 + 5859*x^2813 + 4391*x^2812 + 4927*x^2811 + 2882*x^2810 + 1370*x^2809 + 5556*x^2808 + 670*x^2807 + 5825*x^2806 + 1431*x^2805 + 2939*x^2804 + 3054*x^2803 + 4539*x^2802 + 5219*x^2801 + 657*x^2800 + 3298*x^2799 + 2352*x^2798 + 1311*x^2797 + 2782*x^2796 + 2640*x^2795 + 4921*x^2794 + 2293*x^2793 + 2144*x^2792 + 168*x^2791 + 5790*x^2790 + 824*x^2789 + 4570*x^2788 + 373*x^2787 + 679*x^2786 + 5552*x^2785 + 619*x^2784 + 2070*x^2783 + 1153*x^2782 + 166*x^2781 + 4459*x^2780 + 4318*x^2779 + 530*x^2778 + 2560*x^2777 + 743*x^2776 + 4942*x^2775 + 95*x^2774 + 5614*x^2773 + 1054*x^2772 + 5917*x^2771 + 5076*x^2770 + 1880*x^2769 + 2627*x^2768 + 5350*x^2767 + 3616*x^2766 + 2337*x^2765 + 5154*x^2764 + 2755*x^2763 + 4209*x^2762 + 545*x^2761 + 4523*x^2760 + 936*x^2759 + 3422*x^2758 + 1902*x^2757 + 109*x^2756 + 1102*x^2755 + 3612*x^2754 + 5097*x^2753 + 4527*x^2752 + 2698*x^2751 + 1123*x^2750 + 1406*x^2749 + 3599*x^2748 + 1155*x^2747 + 5519*x^2746 + 4509*x^2745 + 1489*x^2744 + 2110*x^2743 + 856*x^2742 + 5995*x^2741 + 5718*x^2740 + 4342*x^2739 + 2958*x^2738 + 2720*x^2737 + 4269*x^2736 + 189*x^2735 + 2312*x^2734 + 2591*x^2733 + 622*x^2732 + 5488*x^2731 + 2862*x^2730 + 5136*x^2729 + 5260*x^2728 + 4124*x^2727 + 3583*x^2726 + 3844*x^2725 + 2578*x^2724 + 2743*x^2723 + 4910*x^2722 + 4674*x^2721 + 637*x^2720 + 1406*x^2719 + 3267*x^2718 + 1733*x^2717 + 540*x^2716 + 1014*x^2715 + 3420*x^2714 + 3437*x^2713 + 3548*x^2712 + 925*x^2711 + 5666*x^2710 + 3268*x^2709 + 3386*x^2708 + 5889*x^2707 + 439*x^2706 + 1049*x^2705 + 2320*x^2704 + 4970*x^2703 + 966*x^2702 + 3076*x^2701 + 5180*x^2700 + 4799*x^2699 + 933*x^2698 + 2598*x^2697 + 2415*x^2696 + 5663*x^2695 + 2837*x^2694 + 2537*x^2693 + 712*x^2692 + 3627*x^2691 + 4396*x^2690 + 2705*x^2689 + 3514*x^2688 + 2639*x^2687 + 1905*x^2686 + 5301*x^2685 + 1064*x^2684 + 2164*x^2683 + 4946*x^2682 + 3163*x^2681 + 2711*x^2680 + 1383*x^2679 + 2447*x^2678 + 5416*x^2677 + 292*x^2676 + 3469*x^2675 + 4912*x^2674 + 1399*x^2673 + 335*x^2672 + 3329*x^2671 + 1831*x^2670 + 1921*x^2669 + 560*x^2668 + 345*x^2667 + 4170*x^2666 + 618*x^2665 + 5252*x^2664 + 2036*x^2663 + 1242*x^2662 + 3387*x^2661 + 2468*x^2660 + 3978*x^2659 + 579*x^2658 + 5287*x^2657 + 4025*x^2656 + 5977*x^2655 + 1589*x^2654 + 5792*x^2653 + 5672*x^2652 + 3075*x^2651 + 179*x^2650 + 3433*x^2649 + 559*x^2648 + 747*x^2647 + 1709*x^2646 + 1008*x^2645 + 956*x^2644 + 139*x^2643 + 5157*x^2642 + 1290*x^2641 + 5865*x^2640 + 4423*x^2639 + 4153*x^2638 + 126*x^2637 + 1004*x^2636 + 2112*x^2635 + 5266*x^2634 + 4042*x^2633 + 5824*x^2632 + 4798*x^2631 + 1727*x^2630 + 889*x^2629 + 1234*x^2628 + 4574*x^2627 + 1959*x^2626 + 1536*x^2625 + 5589*x^2624 + 5114*x^2623 + 5351*x^2622 + 277*x^2621 + 5784*x^2620 + 926*x^2619 + 4766*x^2618 + 2773*x^2617 + 532*x^2616 + 1276*x^2615 + 1080*x^2614 + 1744*x^2613 + 3646*x^2612 + 3394*x^2611 + 4102*x^2610 + 2260*x^2609 + 3214*x^2608 + 3143*x^2607 + 3666*x^2606 + 5955*x^2605 + 4782*x^2604 + 3060*x^2603 + 3581*x^2602 + 5517*x^2601 + 3454*x^2600 + 4389*x^2599 + 1787*x^2598 + 5534*x^2597 + 5301*x^2596 + 4473*x^2595 + 5985*x^2594 + 5583*x^2593 + 4467*x^2592 + 5228*x^2591 + 18*x^2590 + 1049*x^2589 + 386*x^2588 + 5960*x^2587 + 4863*x^2586 + 5283*x^2585 + 3833*x^2584 + 5418*x^2583 + 3017*x^2582 + 958*x^2581 + 4471*x^2580 + 2090*x^2579 + 5930*x^2578 + 4048*x^2577 + 2166*x^2576 + 3899*x^2575 + 2605*x^2574 + 147*x^2573 + 4632*x^2572 + 4376*x^2571 + 5602*x^2570 + 4728*x^2569 + 3414*x^2568 + 2114*x^2567 + 54*x^2566 + 2963*x^2565 + 5068*x^2564 + 2962*x^2563 + 3660*x^2562 + 2971*x^2561 + 3348*x^2560 + 1794*x^2559 + 1991*x^2558 + 5848*x^2557 + 2498*x^2556 + 3717*x^2555 + 5491*x^2554 + 1076*x^2553 + 5783*x^2552 + 5173*x^2551 + 613*x^2550 + 1513*x^2549 + 3632*x^2548 + 4213*x^2547 + 5023*x^2546 + 5213*x^2545 + 5434*x^2544 + 1892*x^2543 + 659*x^2542 + 761*x^2541 + 538*x^2540 + 2559*x^2539 + 1423*x^2538 + 2244*x^2537 + 3879*x^2536 + 1465*x^2535 + 1096*x^2534 + 2702*x^2533 + 3873*x^2532 + 4202*x^2531 + 846*x^2530 + 2496*x^2529 + 2395*x^2528 + 3689*x^2527 + 2242*x^2526 + 2493*x^2525 + 2370*x^2524 + 397*x^2523 + 3604*x^2522 + 5659*x^2521 + 6026*x^2520 + 5205*x^2519 + 108*x^2518 + 4010*x^2517 + 4699*x^2516 + 6008*x^2515 + 1931*x^2514 + 5306*x^2513 + 3904*x^2512 + 1347*x^2511 + 3281*x^2510 + 1220*x^2509 + 5240*x^2508 + 6009*x^2507 + 3483*x^2506 + 3729*x^2505 + 2176*x^2504 + 5756*x^2503 + 680*x^2502 + 1116*x^2501 + 5306*x^2500 + 1018*x^2499 + 1294*x^2498 + 2458*x^2497 + 5557*x^2496 + 608*x^2495 + 2606*x^2494 + 5130*x^2493 + 981*x^2492 + 2751*x^2491 + 3511*x^2490 + 5143*x^2489 + 1256*x^2488 + 1689*x^2487 + 3322*x^2486 + 4411*x^2485 + 4067*x^2484 + 1847*x^2483 + 4644*x^2482 + 3971*x^2481 + 1089*x^2480 + 4880*x^2479 + 2872*x^2478 + 3534*x^2477 + 3238*x^2476 + 4164*x^2475 + 2494*x^2474 + 39*x^2473 + 5741*x^2472 + 1231*x^2471 + 1109*x^2470 + 5315*x^2469 + 3657*x^2468 + 5201*x^2467 + 344*x^2466 + 4358*x^2465 + 1333*x^2464 + 5244*x^2463 + 198*x^2462 + 1712*x^2461 + 2164*x^2460 + 3008*x^2459 + 4841*x^2458 + 1976*x^2457 + 737*x^2456 + 5032*x^2455 + 2750*x^2454 + 5223*x^2453 + 3925*x^2452 + 2914*x^2451 + 4173*x^2450 + 1192*x^2449 + 3122*x^2448 + 4354*x^2447 + 2821*x^2446 + 601*x^2445 + 2212*x^2444 + 2036*x^2443 + 1836*x^2442 + 4242*x^2441 + 1968*x^2440 + 4659*x^2439 + 2201*x^2438 + 736*x^2437 + 2005*x^2436 + 1278*x^2435 + 2784*x^2434 + 4108*x^2433 + 1607*x^2432 + 1150*x^2431 + 56*x^2430 + 5452*x^2429 + 2219*x^2428 + 1864*x^2427 + 4609*x^2426 + 4003*x^2425 + 1951*x^2424 + 4407*x^2423 + 397*x^2422 + 4238*x^2421 + 5791*x^2420 + 5481*x^2419 + 2454*x^2418 + 4414*x^2417 + 1556*x^2416 + 5322*x^2415 + 2857*x^2414 + 1929*x^2413 + 5169*x^2412 + 2636*x^2411 + 2972*x^2410 + 3242*x^2409 + 3540*x^2408 + 4742*x^2407 + 4007*x^2406 + 1842*x^2405 + 5005*x^2404 + 5348*x^2403 + 1771*x^2402 + 1516*x^2401 + 5180*x^2400 + 5882*x^2399 + 5927*x^2398 + 4828*x^2397 + 4124*x^2396 + 3194*x^2395 + 2776*x^2394 + 1018*x^2393 + 5354*x^2392 + 1730*x^2391 + 374*x^2390 + 5844*x^2389 + 3384*x^2388 + 890*x^2387 + 1329*x^2386 + 5767*x^2385 + 4331*x^2384 + 3918*x^2383 + 268*x^2382 + 1271*x^2381 + 5707*x^2380 + 3124*x^2379 + 3376*x^2378 + 1568*x^2377 + 1662*x^2376 + 3732*x^2375 + 6018*x^2374 + 3916*x^2373 + 2570*x^2372 + 2437*x^2371 + 3358*x^2370 + 1652*x^2369 + 618*x^2368 + 1125*x^2367 + 5726*x^2366 + 1990*x^2365 + 2220*x^2364 + 3608*x^2363 + 5850*x^2362 + 1741*x^2361 + 5659*x^2360 + 1444*x^2359 + 3098*x^2358 + 122*x^2357 + 3347*x^2356 + 3215*x^2355 + 4149*x^2354 + 35*x^2353 + 1823*x^2352 + 2012*x^2351 + 329*x^2350 + 4226*x^2349 + 366*x^2348 + 2286*x^2347 + 158*x^2346 + 5706*x^2345 + 3516*x^2344 + 3124*x^2343 + 3464*x^2342 + 2984*x^2341 + 73*x^2340 + 5144*x^2339 + 2094*x^2338 + 72*x^2337 + 2821*x^2336 + 2782*x^2335 + 2194*x^2334 + 2532*x^2333 + 2436*x^2332 + 2678*x^2331 + 2470*x^2330 + 3675*x^2329 + 734*x^2328 + 2855*x^2327 + 3930*x^2326 + 4351*x^2325 + 2660*x^2324 + 2142*x^2323 + 1784*x^2322 + 2218*x^2321 + 3550*x^2320 + 1974*x^2319 + 5209*x^2318 + 4395*x^2317 + 5352*x^2316 + 1318*x^2315 + 3576*x^2314 + 2022*x^2313 + 3794*x^2312 + 680*x^2311 + 3499*x^2310 + 3826*x^2309 + 3021*x^2308 + 4888*x^2307 + 3121*x^2306 + 4442*x^2305 + 1134*x^2304 + 3555*x^2303 + 2096*x^2302 + 1054*x^2301 + 1162*x^2300 + 4593*x^2299 + 1839*x^2298 + 466*x^2297 + 1815*x^2296 + 48*x^2295 + 4328*x^2294 + 5520*x^2293 + 3155*x^2292 + 3725*x^2291 + 5733*x^2290 + 1828*x^2289 + 4878*x^2288 + 3862*x^2287 + 3586*x^2286 + 1064*x^2285 + 4872*x^2284 + 4550*x^2283 + 3731*x^2282 + 3899*x^2281 + 5219*x^2280 + 777*x^2279 + 5111*x^2278 + 3467*x^2277 + 2483*x^2276 + 5895*x^2275 + 532*x^2274 + 4089*x^2273 + 4135*x^2272 + 5034*x^2271 + 2511*x^2270 + 5366*x^2269 + 3954*x^2268 + 5301*x^2267 + 3368*x^2266 + 2070*x^2265 + 4320*x^2264 + 5742*x^2263 + 1223*x^2262 + 5604*x^2261 + 1272*x^2260 + 2839*x^2259 + 2256*x^2258 + 5229*x^2257 + 478*x^2256 + 2970*x^2255 + 2240*x^2254 + 4099*x^2253 + 2521*x^2252 + 2725*x^2251 + 3725*x^2250 + 910*x^2249 + 988*x^2248 + 5563*x^2247 + 3163*x^2246 + 4603*x^2245 + 4666*x^2244 + 2713*x^2243 + 4555*x^2242 + 2339*x^2241 + 2847*x^2240 + 1610*x^2239 + 5207*x^2238 + 1813*x^2237 + 3753*x^2236 + 2398*x^2235 + 1877*x^2234 + 2754*x^2233 + 1004*x^2232 + 3594*x^2231 + 2483*x^2230 + 247*x^2229 + 2150*x^2228 + 2071*x^2227 + 4619*x^2226 + 5708*x^2225 + 4706*x^2224 + 4018*x^2223 + 4179*x^2222 + 3808*x^2221 + 5918*x^2220 + 3689*x^2219 + 668*x^2218 + 2615*x^2217 + 3771*x^2216 + 4617*x^2215 + 658*x^2214 + 5187*x^2213 + 5897*x^2212 + 1876*x^2211 + 5243*x^2210 + 3381*x^2209 + 2481*x^2208 + 1768*x^2207 + 1389*x^2206 + 1245*x^2205 + 2814*x^2204 + 2298*x^2203 + 5377*x^2202 + 4755*x^2201 + 5550*x^2200 + 1902*x^2199 + 1986*x^2198 + 2096*x^2197 + 5174*x^2196 + 1844*x^2195 + 1787*x^2194 + 2838*x^2193 + 5812*x^2192 + 76*x^2191 + 5451*x^2190 + 5534*x^2189 + 3005*x^2188 + 2378*x^2187 + 5233*x^2186 + 3273*x^2185 + 3940*x^2184 + 2161*x^2183 + 356*x^2182 + 4480*x^2181 + 2543*x^2180 + 4732*x^2179 + 4818*x^2178 + 3537*x^2177 + 4869*x^2176 + 163*x^2175 + 1533*x^2174 + 1818*x^2173 + 3858*x^2172 + 4415*x^2171 + 3736*x^2170 + 5458*x^2169 + 4662*x^2168 + 1307*x^2167 + 536*x^2166 + 2798*x^2165 + 1201*x^2164 + 838*x^2163 + 402*x^2162 + 5501*x^2161 + 2034*x^2160 + 1355*x^2159 + 1702*x^2158 + 3144*x^2157 + 3802*x^2156 + 3448*x^2155 + 3007*x^2154 + 137*x^2153 + 2089*x^2152 + 3939*x^2151 + 4297*x^2150 + 3749*x^2149 + 2304*x^2148 + 208*x^2147 + 375*x^2146 + 1609*x^2145 + 3350*x^2144 + 3101*x^2143 + 5175*x^2142 + 567*x^2141 + 4237*x^2140 + 3708*x^2139 + 2463*x^2138 + 911*x^2137 + 2138*x^2136 + 2577*x^2135 + 4851*x^2134 + 3478*x^2133 + 5099*x^2132 + 1908*x^2131 + 1784*x^2130 + 5533*x^2129 + 2936*x^2128 + 5491*x^2127 + 3593*x^2126 + 707*x^2125 + 1578*x^2124 + 2173*x^2123 + 1429*x^2122 + 5217*x^2121 + 1140*x^2120 + 4849*x^2119 + 3769*x^2118 + 1660*x^2117 + 958*x^2116 + 2356*x^2115 + 3931*x^2114 + 500*x^2113 + 1260*x^2112 + 2225*x^2111 + 909*x^2110 + 520*x^2109 + 3340*x^2108 + 4497*x^2107 + 4660*x^2106 + 2750*x^2105 + 4819*x^2104 + 4662*x^2103 + 5038*x^2102 + 1356*x^2101 + 3190*x^2100 + 2792*x^2099 + 576*x^2098 + 2982*x^2097 + 2165*x^2096 + 4967*x^2095 + 3534*x^2094 + 2144*x^2093 + 1869*x^2092 + 4416*x^2091 + 5232*x^2090 + 1778*x^2089 + 639*x^2088 + 3362*x^2087 + 3243*x^2086 + 4518*x^2085 + 5876*x^2084 + 5268*x^2083 + 1388*x^2082 + 3833*x^2081 + 4183*x^2080 + 1975*x^2079 + 4141*x^2078 + 2567*x^2077 + 494*x^2076 + 5478*x^2075 + 3116*x^2074 + 784*x^2073 + 4750*x^2072 + 5578*x^2071 + 1681*x^2070 + 3162*x^2069 + 5329*x^2068 + 954*x^2067 + 1630*x^2066 + 4916*x^2065 + 187*x^2064 + 4952*x^2063 + 3446*x^2062 + 4958*x^2061 + 1144*x^2060 + 4277*x^2059 + 1864*x^2058 + 3386*x^2057 + 4181*x^2056 + 2996*x^2055 + 637*x^2054 + 310*x^2053 + 3629*x^2052 + 5310*x^2051 + 4836*x^2050 + 3677*x^2049 + 536*x^2048 + 2536*x^2047 + 4484*x^2046 + 3939*x^2045 + 5046*x^2044 + 5821*x^2043 + 5748*x^2042 + 3046*x^2041 + 1489*x^2040 + 3561*x^2039 + 4516*x^2038 + 1905*x^2037 + 3529*x^2036 + 1036*x^2035 + 3206*x^2034 + 4809*x^2033 + 218*x^2032 + 2093*x^2031 + 129*x^2030 + 2335*x^2029 + 3990*x^2028 + 492*x^2027 + 2929*x^2026 + 4713*x^2025 + 4652*x^2024 + 2707*x^2023 + 2445*x^2022 + 2618*x^2021 + 2129*x^2020 + 3079*x^2019 + 4272*x^2018 + 947*x^2017 + 1640*x^2016 + 1922*x^2015 + 3790*x^2014 + 763*x^2013 + 2932*x^2012 + 5549*x^2011 + 486*x^2010 + 2947*x^2009 + 1191*x^2008 + 4901*x^2007 + 181*x^2006 + 3010*x^2005 + 684*x^2004 + 2321*x^2003 + 431*x^2002 + 5630*x^2001 + 4568*x^2000 + 5695*x^1999 + 1989*x^1998 + 2347*x^1997 + 5762*x^1996 + 2169*x^1995 + 4664*x^1994 + 3300*x^1993 + 3151*x^1992 + 5285*x^1991 + 3148*x^1990 + 4777*x^1989 + 5161*x^1988 + 1116*x^1987 + 4802*x^1986 + 1754*x^1985 + 3725*x^1984 + 2989*x^1983 + 4675*x^1982 + 1797*x^1981 + 4798*x^1980 + 3616*x^1979 + 2221*x^1978 + 5406*x^1977 + 624*x^1976 + 4360*x^1975 + 4408*x^1974 + 5429*x^1973 + 3375*x^1972 + 4361*x^1971 + 823*x^1970 + 5102*x^1969 + 5232*x^1968 + 1782*x^1967 + 5186*x^1966 + 4106*x^1965 + 409*x^1964 + 3233*x^1963 + 4345*x^1962 + 3327*x^1961 + 4444*x^1960 + 4662*x^1959 + 803*x^1958 + 553*x^1957 + 1186*x^1956 + 496*x^1955 + 4542*x^1954 + 529*x^1953 + 474*x^1952 + 5888*x^1951 + 4216*x^1950 + 5920*x^1949 + 3912*x^1948 + 5752*x^1947 + 920*x^1946 + 595*x^1945 + 4551*x^1944 + 5560*x^1943 + 2259*x^1942 + 5865*x^1941 + 332*x^1940 + 629*x^1939 + 1601*x^1938 + 3729*x^1937 + 3282*x^1936 + 2519*x^1935 + 1947*x^1934 + 2138*x^1933 + 3739*x^1932 + 4867*x^1931 + 4541*x^1930 + 3013*x^1929 + 4555*x^1928 + 4654*x^1927 + 496*x^1926 + 1814*x^1925 + 4207*x^1924 + 2237*x^1923 + 5061*x^1922 + 2085*x^1921 + 5228*x^1920 + 740*x^1919 + 1879*x^1918 + 3817*x^1917 + 4581*x^1916 + 3057*x^1915 + 5620*x^1914 + 2135*x^1913 + 321*x^1912 + 3367*x^1911 + 1011*x^1910 + 1056*x^1909 + 549*x^1908 + 1467*x^1907 + 1461*x^1906 + 2681*x^1905 + 663*x^1904 + 1304*x^1903 + 830*x^1902 + 4526*x^1901 + 2214*x^1900 + 1489*x^1899 + 899*x^1898 + 2864*x^1897 + 233*x^1896 + 3189*x^1895 + 949*x^1894 + 3548*x^1893 + 2338*x^1892 + 17*x^1891 + 3821*x^1890 + 2551*x^1889 + 5896*x^1888 + 4600*x^1887 + 4477*x^1886 + 5300*x^1885 + 1113*x^1884 + 5465*x^1883 + 3165*x^1882 + 3367*x^1881 + 687*x^1880 + 5079*x^1879 + 1224*x^1878 + 3001*x^1877 + 1749*x^1876 + 4978*x^1875 + 639*x^1874 + 183*x^1873 + 1360*x^1872 + 5319*x^1871 + 5298*x^1870 + 4403*x^1869 + 239*x^1868 + 133*x^1867 + 4003*x^1866 + 2221*x^1865 + 2576*x^1864 + 5777*x^1863 + 3858*x^1862 + 2956*x^1861 + 4262*x^1860 + 4050*x^1859 + 5582*x^1858 + 5683*x^1857 + 1072*x^1856 + 2480*x^1855 + 506*x^1854 + 4289*x^1853 + 1697*x^1852 + 2468*x^1851 + 5471*x^1850 + 2346*x^1849 + 1765*x^1848 + 5627*x^1847 + 767*x^1846 + 1516*x^1845 + 4446*x^1844 + 5626*x^1843 + 5226*x^1842 + 1623*x^1841 + 2204*x^1840 + 612*x^1839 + 1131*x^1838 + 947*x^1837 + 3319*x^1836 + 3884*x^1835 + 3212*x^1834 + 3777*x^1833 + 918*x^1832 + 1224*x^1831 + 1873*x^1830 + 5781*x^1829 + 2466*x^1828 + 2321*x^1827 + 4856*x^1826 + 5444*x^1825 + 5499*x^1824 + 4988*x^1823 + 872*x^1822 + 2105*x^1821 + 3897*x^1820 + 1770*x^1819 + 1267*x^1818 + 2367*x^1817 + 4480*x^1816 + 3560*x^1815 + 202*x^1814 + 104*x^1813 + 5641*x^1812 + 5727*x^1811 + 2984*x^1810 + 1482*x^1809 + 1557*x^1808 + 2327*x^1807 + 4392*x^1806 + 4576*x^1805 + 5747*x^1804 + 376*x^1803 + 1580*x^1802 + 963*x^1801 + 5371*x^1800 + 3734*x^1799 + 3948*x^1798 + 5279*x^1797 + 3621*x^1796 + 342*x^1795 + 4984*x^1794 + 5424*x^1793 + 143*x^1792 + 3891*x^1791 + 2206*x^1790 + 40*x^1789 + 3614*x^1788 + 3581*x^1787 + 4030*x^1786 + 1086*x^1785 + 1973*x^1784 + 666*x^1783 + 1389*x^1782 + 4188*x^1781 + 398*x^1780 + 5430*x^1779 + 396*x^1778 + 4839*x^1777 + 326*x^1776 + 4853*x^1775 + 4081*x^1774 + 2408*x^1773 + 247*x^1772 + 2726*x^1771 + 3856*x^1770 + 3023*x^1769 + 5244*x^1768 + 5303*x^1767 + 482*x^1766 + 4223*x^1765 + 3411*x^1764 + 1578*x^1763 + 1743*x^1762 + 139*x^1761 + 1409*x^1760 + 2633*x^1759 + 3262*x^1758 + 2011*x^1757 + 3383*x^1756 + 1189*x^1755 + 4236*x^1754 + 2140*x^1753 + 860*x^1752 + 239*x^1751 + 2294*x^1750 + 3299*x^1749 + 3463*x^1748 + 1581*x^1747 + 4128*x^1746 + 785*x^1745 + 1933*x^1744 + 3089*x^1743 + 4706*x^1742 + 5925*x^1741 + 5561*x^1740 + 3470*x^1739 + 1792*x^1738 + 4941*x^1737 + 4903*x^1736 + 310*x^1735 + 3859*x^1734 + 1790*x^1733 + 3052*x^1732 + 5808*x^1731 + 5961*x^1730 + 1440*x^1729 + 3035*x^1728 + 2656*x^1727 + 1178*x^1726 + 1776*x^1725 + 1898*x^1724 + 2908*x^1723 + 439*x^1722 + 3108*x^1721 + 5163*x^1720 + 4359*x^1719 + 3726*x^1718 + 59*x^1717 + 121*x^1716 + 988*x^1715 + 2489*x^1714 + 4551*x^1713 + 4564*x^1712 + 959*x^1711 + 1083*x^1710 + 3286*x^1709 + 5465*x^1708 + 829*x^1707 + 517*x^1706 + 2834*x^1705 + 4244*x^1704 + 936*x^1703 + 1588*x^1702 + 914*x^1701 + 1312*x^1700 + 3103*x^1699 + 1434*x^1698 + 4362*x^1697 + 5769*x^1696 + 5948*x^1695 + 5558*x^1694 + 179*x^1693 + 5730*x^1692 + 5781*x^1691 + 392*x^1690 + 234*x^1689 + 4123*x^1688 + 5327*x^1687 + 3581*x^1686 + 4998*x^1685 + 1518*x^1684 + 4906*x^1683 + 3184*x^1682 + 3883*x^1681 + 1195*x^1680 + 4003*x^1679 + 5885*x^1678 + 3712*x^1677 + 1474*x^1676 + 4897*x^1675 + 3476*x^1674 + 1387*x^1673 + 39*x^1672 + 4565*x^1671 + 5219*x^1670 + 5235*x^1669 + 4412*x^1668 + 4277*x^1667 + 4992*x^1666 + 4812*x^1665 + 5260*x^1664 + 1693*x^1663 + 1751*x^1662 + 509*x^1661 + 3366*x^1660 + 3144*x^1659 + 102*x^1658 + 3898*x^1657 + 2626*x^1656 + 3853*x^1655 + 3758*x^1654 + 1662*x^1653 + 1645*x^1652 + 4143*x^1651 + 3628*x^1650 + 3596*x^1649 + 5740*x^1648 + 399*x^1647 + 4876*x^1646 + 4074*x^1645 + 41*x^1644 + 4841*x^1643 + 2490*x^1642 + 2774*x^1641 + 309*x^1640 + 3727*x^1639 + 4823*x^1638 + 5068*x^1637 + 3198*x^1636 + 5458*x^1635 + 598*x^1634 + 2723*x^1633 + 530*x^1632 + 507*x^1631 + 720*x^1630 + 2861*x^1629 + 5611*x^1628 + 1090*x^1627 + 3706*x^1626 + 5985*x^1625 + 675*x^1624 + 2158*x^1623 + 451*x^1622 + 500*x^1621 + 6032*x^1620 + 5109*x^1619 + 981*x^1618 + 4246*x^1617 + 2794*x^1616 + 1043*x^1615 + 3556*x^1614 + 1079*x^1613 + 5077*x^1612 + 5875*x^1611 + 900*x^1610 + 887*x^1609 + 5677*x^1608 + 3792*x^1607 + 4985*x^1606 + 2727*x^1605 + 4031*x^1604 + 1684*x^1603 + 2358*x^1602 + 4184*x^1601 + 84*x^1600 + 606*x^1599 + 3441*x^1598 + 4866*x^1597 + 3676*x^1596 + 1739*x^1595 + 4059*x^1594 + 5581*x^1593 + 5767*x^1592 + 134*x^1591 + 2725*x^1590 + 4013*x^1589 + 3802*x^1588 + 5105*x^1587 + 4343*x^1586 + 156*x^1585 + 460*x^1584 + 2825*x^1583 + 2214*x^1582 + 1589*x^1581 + 5336*x^1580 + 134*x^1579 + 3144*x^1578 + 3702*x^1577 + 5930*x^1576 + 5659*x^1575 + 23*x^1574 + 2527*x^1573 + 3986*x^1572 + 3849*x^1571 + 5555*x^1570 + 1705*x^1569 + 981*x^1568 + 4753*x^1567 + 4804*x^1566 + 3975*x^1565 + 4753*x^1564 + 3806*x^1563 + 5508*x^1562 + 4741*x^1561 + 2053*x^1560 + 3132*x^1559 + 458*x^1558 + 5420*x^1557 + 5185*x^1556 + 2375*x^1555 + 4389*x^1554 + 3111*x^1553 + 2693*x^1552 + 2995*x^1551 + 5591*x^1550 + 5458*x^1549 + 966*x^1548 + 2034*x^1547 + 2173*x^1546 + 5002*x^1545 + 4296*x^1544 + 616*x^1543 + 1317*x^1542 + 4667*x^1541 + 1681*x^1540 + 3333*x^1539 + 2672*x^1538 + 3040*x^1537 + 1740*x^1536 + 5079*x^1535 + 3281*x^1534 + 4190*x^1533 + 2777*x^1532 + 2320*x^1531 + 4045*x^1530 + 4200*x^1529 + 2372*x^1528 + 262*x^1527 + 4061*x^1526 + 5930*x^1525 + 3475*x^1524 + 4833*x^1523 + 3155*x^1522 + 1574*x^1521 + 698*x^1520 + 1534*x^1519 + 1904*x^1518 + 5565*x^1517 + 5043*x^1516 + 4600*x^1515 + 1787*x^1514 + 5565*x^1513 + 4600*x^1512 + 4098*x^1511 + 5122*x^1510 + 5823*x^1509 + 16*x^1508 + 1243*x^1507 + 146*x^1506 + 2177*x^1505 + 1275*x^1504 + 3601*x^1503 + 3892*x^1502 + 5149*x^1501 + 2634*x^1500 + 1208*x^1499 + 1604*x^1498 + 3496*x^1497 + 2158*x^1496 + 1786*x^1495 + 1731*x^1494 + 3219*x^1493 + 5120*x^1492 + 4527*x^1491 + 4076*x^1490 + 5914*x^1489 + 5421*x^1488 + 5034*x^1487 + 892*x^1486 + 5809*x^1485 + 5541*x^1484 + 860*x^1483 + 3763*x^1482 + 2480*x^1481 + 4217*x^1480 + 3489*x^1479 + 5904*x^1478 + 5790*x^1477 + 5445*x^1476 + 1769*x^1475 + 2132*x^1474 + 1103*x^1473 + 1109*x^1472 + 3236*x^1471 + 554*x^1470 + 5750*x^1469 + 4095*x^1468 + 146*x^1467 + 5450*x^1466 + 211*x^1465 + 4392*x^1464 + 1331*x^1463 + 839*x^1462 + 2937*x^1461 + 2469*x^1460 + 2856*x^1459 + 4466*x^1458 + 4944*x^1457 + 4165*x^1456 + 2610*x^1455 + 3114*x^1454 + 3630*x^1453 + 2191*x^1452 + 5157*x^1451 + 3014*x^1450 + 4470*x^1449 + 2430*x^1448 + 2772*x^1447 + 2051*x^1446 + 4258*x^1445 + 5612*x^1444 + 3792*x^1443 + 2446*x^1442 + 2670*x^1441 + 3185*x^1440 + 5433*x^1439 + 5053*x^1438 + 2679*x^1437 + 6023*x^1436 + 4507*x^1435 + 2294*x^1434 + 1177*x^1433 + 3759*x^1432 + 1218*x^1431 + 1742*x^1430 + 3028*x^1429 + 2930*x^1428 + 4013*x^1427 + 3071*x^1426 + 1703*x^1425 + 2816*x^1424 + 3518*x^1423 + 130*x^1422 + 2254*x^1421 + 308*x^1420 + 4341*x^1419 + 3393*x^1418 + 4633*x^1417 + 2048*x^1416 + 3434*x^1415 + 948*x^1414 + 3075*x^1413 + 1874*x^1412 + 668*x^1411 + 655*x^1410 + 2385*x^1409 + 5449*x^1408 + 977*x^1407 + 3685*x^1406 + 2458*x^1405 + 3454*x^1404 + 3248*x^1403 + 4447*x^1402 + 4394*x^1401 + 1973*x^1400 + 4200*x^1399 + 5985*x^1398 + 1439*x^1397 + 334*x^1396 + 4690*x^1395 + 742*x^1394 + 2392*x^1393 + 1959*x^1392 + 6032*x^1391 + 3223*x^1390 + 630*x^1389 + 2934*x^1388 + 90*x^1387 + 4819*x^1386 + 2470*x^1385 + 5224*x^1384 + 91*x^1383 + 5631*x^1382 + 4826*x^1381 + 2817*x^1380 + 2659*x^1379 + 5231*x^1378 + 4445*x^1377 + 3447*x^1376 + 1067*x^1375 + 3953*x^1374 + 5727*x^1373 + 2777*x^1372 + 5862*x^1371 + 2127*x^1370 + 458*x^1369 + 155*x^1368 + 2705*x^1367 + 5539*x^1366 + 1364*x^1365 + 2472*x^1364 + 5012*x^1363 + 5943*x^1362 + 3352*x^1361 + 3709*x^1360 + 2856*x^1359 + 2631*x^1358 + 3523*x^1357 + 680*x^1356 + 1353*x^1355 + 5072*x^1354 + 5909*x^1353 + 1738*x^1352 + 140*x^1351 + 899*x^1350 + 5082*x^1349 + 4593*x^1348 + 3293*x^1347 + 2599*x^1346 + 5224*x^1345 + 2179*x^1344 + 5092*x^1343 + 3911*x^1342 + 3187*x^1341 + 4059*x^1340 + 1652*x^1339 + 3681*x^1338 + 5379*x^1337 + 5724*x^1336 + 4116*x^1335 + 1132*x^1334 + 4683*x^1333 + 3895*x^1332 + 5682*x^1331 + 5189*x^1330 + 4438*x^1329 + 4558*x^1328 + 3913*x^1327 + 3627*x^1326 + 131*x^1325 + 6010*x^1324 + 4226*x^1323 + 5811*x^1322 + 207*x^1321 + 3959*x^1320 + 2109*x^1319 + 1821*x^1318 + 587*x^1317 + 1798*x^1316 + 2228*x^1315 + 259*x^1314 + 2384*x^1313 + 2125*x^1312 + 206*x^1311 + 1699*x^1310 + 1539*x^1309 + 3754*x^1308 + 291*x^1307 + 732*x^1306 + 2901*x^1305 + 4896*x^1304 + 2986*x^1303 + 3818*x^1302 + 5115*x^1301 + 3019*x^1300 + 4497*x^1299 + 1965*x^1298 + 4382*x^1297 + 4964*x^1296 + 5054*x^1295 + 1662*x^1294 + 1677*x^1293 + 5099*x^1292 + 3924*x^1291 + 329*x^1290 + 2885*x^1289 + 4749*x^1288 + 1093*x^1287 + 5537*x^1286 + 1462*x^1285 + 5455*x^1284 + 2034*x^1283 + 5222*x^1282 + 4488*x^1281 + 3936*x^1280 + 3453*x^1279 + 4174*x^1278 + 2137*x^1277 + 5813*x^1276 + 2845*x^1275 + 4262*x^1274 + 5703*x^1273 + 1421*x^1272 + 1194*x^1271 + 298*x^1270 + 3542*x^1269 + 1653*x^1268 + 2913*x^1267 + 3221*x^1266 + 4661*x^1265 + 2894*x^1264 + 3876*x^1263 + 2490*x^1262 + 2365*x^1261 + 3663*x^1260 + 860*x^1259 + 5335*x^1258 + 5192*x^1257 + 1050*x^1256 + 3638*x^1255 + 3850*x^1254 + 1140*x^1253 + 1634*x^1252 + 4765*x^1251 + 4247*x^1250 + 2401*x^1249 + 4764*x^1248 + 2675*x^1247 + 5748*x^1246 + 1635*x^1245 + 3457*x^1244 + 1621*x^1243 + 83*x^1242 + 3970*x^1241 + 4963*x^1240 + 3833*x^1239 + 3263*x^1238 + 4616*x^1237 + 5035*x^1236 + 3926*x^1235 + 2904*x^1234 + 5276*x^1233 + 501*x^1232 + 5337*x^1231 + 3076*x^1230 + 2139*x^1229 + 1179*x^1228 + 3545*x^1227 + 5217*x^1226 + 393*x^1225 + 5791*x^1224 + 5683*x^1223 + 5592*x^1222 + 5202*x^1221 + 5593*x^1220 + 2298*x^1219 + 5564*x^1218 + 1450*x^1217 + 4605*x^1216 + 820*x^1215 + 5589*x^1214 + 2029*x^1213 + 3133*x^1212 + 666*x^1211 + 1958*x^1210 + 973*x^1209 + 4867*x^1208 + 5111*x^1207 + 2963*x^1206 + 3954*x^1205 + 1681*x^1204 + 5215*x^1203 + 3132*x^1202 + 3514*x^1201 + 3159*x^1200 + 1255*x^1199 + 4042*x^1198 + 2651*x^1197 + 3416*x^1196 + 4718*x^1195 + 349*x^1194 + 1022*x^1193 + 4959*x^1192 + 1413*x^1191 + 4570*x^1190 + 672*x^1189 + 4970*x^1188 + 237*x^1187 + 3963*x^1186 + 2164*x^1185 + 5570*x^1184 + 1576*x^1183 + 983*x^1182 + 2843*x^1181 + 134*x^1180 + 906*x^1179 + 1917*x^1177 + 1092*x^1176 + 4962*x^1175 + 3352*x^1174 + 3102*x^1173 + 2991*x^1172 + 4737*x^1171 + 2985*x^1170 + 157*x^1169 + 3736*x^1168 + 2412*x^1167 + 1337*x^1166 + 2577*x^1165 + 5207*x^1164 + 2912*x^1163 + 237*x^1162 + 990*x^1161 + 4416*x^1160 + 2835*x^1159 + 3953*x^1158 + 4460*x^1157 + 3509*x^1156 + 5396*x^1155 + 5507*x^1154 + 463*x^1153 + 6*x^1152 + 1247*x^1151 + 2060*x^1150 + 5684*x^1149 + 4211*x^1148 + 1510*x^1147 + 3658*x^1146 + 1791*x^1145 + 3919*x^1144 + 1011*x^1143 + 2786*x^1142 + 1277*x^1141 + 1253*x^1140 + 3001*x^1139 + 355*x^1138 + 3277*x^1137 + 913*x^1136 + 4102*x^1135 + 2397*x^1134 + 1017*x^1133 + 556*x^1132 + 1032*x^1131 + 2996*x^1130 + 5577*x^1129 + 3450*x^1128 + 1280*x^1127 + 174*x^1126 + 2806*x^1125 + 18*x^1124 + 5393*x^1123 + 1358*x^1122 + 4051*x^1121 + 2712*x^1120 + 3132*x^1119 + 5334*x^1118 + 1234*x^1117 + 1484*x^1116 + 3210*x^1115 + 404*x^1114 + 1463*x^1113 + 5156*x^1112 + 3181*x^1111 + 5122*x^1110 + 5723*x^1109 + 2218*x^1108 + 5311*x^1107 + 1261*x^1106 + 531*x^1105 + 5761*x^1104 + 5777*x^1103 + 1750*x^1102 + 2800*x^1101 + 956*x^1100 + 1776*x^1099 + 4705*x^1098 + 6011*x^1097 + 2677*x^1096 + 2588*x^1095 + 3399*x^1094 + 6021*x^1093 + 963*x^1092 + 3639*x^1091 + 1251*x^1090 + 3864*x^1089 + 2552*x^1088 + 3765*x^1087 + 697*x^1086 + 4484*x^1085 + 4231*x^1084 + 3151*x^1083 + 217*x^1082 + 1141*x^1081 + 1252*x^1080 + 1135*x^1079 + 2160*x^1078 + 3886*x^1077 + 2130*x^1076 + 4161*x^1075 + 175*x^1074 + 462*x^1073 + 3261*x^1072 + 5374*x^1071 + 3735*x^1070 + 4909*x^1069 + 553*x^1068 + 4088*x^1067 + 1601*x^1066 + 50*x^1065 + 5927*x^1064 + 5658*x^1063 + 2919*x^1062 + 95*x^1061 + 5083*x^1060 + 3223*x^1059 + 5259*x^1058 + 4177*x^1057 + 5564*x^1056 + 3429*x^1055 + 2633*x^1054 + 1914*x^1053 + 5381*x^1052 + 4031*x^1051 + 3949*x^1050 + 2083*x^1049 + 2955*x^1048 + 757*x^1047 + 1371*x^1046 + 356*x^1045 + 1809*x^1044 + 2411*x^1043 + 2498*x^1042 + 3179*x^1041 + 5837*x^1040 + 3250*x^1039 + 1211*x^1038 + 1864*x^1037 + 2413*x^1036 + 2384*x^1035 + 4422*x^1034 + 1508*x^1033 + 1943*x^1032 + 3037*x^1031 + 1562*x^1030 + 1131*x^1029 + 5636*x^1028 + 4162*x^1027 + 2748*x^1026 + 2254*x^1025 + 3837*x^1024 + 3317*x^1023 + 5364*x^1022 + 655*x^1021 + 1323*x^1020 + 5520*x^1019 + 762*x^1018 + 4123*x^1017 + 609*x^1016 + 2882*x^1015 + 3180*x^1014 + 1749*x^1013 + 4266*x^1012 + 1494*x^1011 + 5859*x^1010 + 4782*x^1009 + 5665*x^1008 + 3958*x^1007 + 3265*x^1006 + 4872*x^1005 + 4166*x^1004 + 8*x^1003 + 2655*x^1002 + 968*x^1001 + 2021*x^1000 + 1955*x^999 + 5965*x^998 + 290*x^997 + 709*x^996 + 3985*x^995 + 2013*x^994 + 469*x^993 + 2381*x^992 + 288*x^991 + 1351*x^990 + 2825*x^989 + 3621*x^988 + 341*x^987 + 1541*x^986 + 1167*x^985 + 2211*x^984 + 817*x^983 + 3707*x^982 + 4505*x^981 + 477*x^980 + 4625*x^979 + 299*x^978 + 5105*x^977 + 5231*x^976 + 1610*x^975 + 245*x^974 + 2062*x^973 + 3819*x^972 + 3735*x^971 + 880*x^970 + 621*x^969 + 3069*x^968 + 2255*x^967 + 3461*x^966 + 3447*x^965 + 1516*x^964 + 560*x^963 + 4560*x^962 + 2891*x^961 + 6016*x^960 + 4810*x^959 + 3516*x^958 + 1474*x^957 + 2248*x^956 + 2676*x^955 + 5286*x^954 + 3267*x^953 + 5*x^952 + 3366*x^951 + 98*x^950 + 3957*x^949 + 1669*x^948 + 3557*x^947 + 3669*x^946 + 3889*x^945 + 358*x^944 + 5780*x^943 + 1913*x^942 + 3996*x^941 + 4151*x^940 + 2995*x^939 + 3246*x^938 + 497*x^937 + 1738*x^936 + 3638*x^935 + 3461*x^934 + 1225*x^933 + 3920*x^932 + 2768*x^931 + 3034*x^930 + 3239*x^929 + 3465*x^928 + 2996*x^927 + 1819*x^926 + 3253*x^925 + 2496*x^924 + 590*x^923 + 1948*x^922 + 2714*x^921 + 1020*x^920 + 4805*x^919 + 1466*x^918 + 5610*x^917 + 2448*x^916 + 198*x^915 + 61*x^914 + 4987*x^913 + 3995*x^912 + 1582*x^911 + 5791*x^910 + 3629*x^909 + 5773*x^908 + 5190*x^907 + 2912*x^906 + 5779*x^905 + 1839*x^904 + 3037*x^903 + 4700*x^902 + 4151*x^901 + 4143*x^900 + 1208*x^899 + 1897*x^898 + 2497*x^897 + 51*x^896 + 442*x^895 + 5492*x^894 + 4755*x^893 + 5795*x^892 + 5525*x^891 + 3420*x^890 + 3206*x^889 + 1631*x^888 + 3705*x^887 + 4148*x^886 + 1859*x^885 + 2488*x^884 + 4611*x^883 + 3473*x^882 + 5466*x^881 + 4254*x^880 + 3123*x^879 + 3976*x^878 + 4567*x^877 + 1969*x^876 + 4076*x^875 + 1379*x^874 + 2843*x^873 + 194*x^872 + 4258*x^871 + 3243*x^870 + 5022*x^869 + 1672*x^868 + 2454*x^867 + 4065*x^866 + 2810*x^865 + 5214*x^864 + 2296*x^863 + 1453*x^862 + 1207*x^861 + 629*x^860 + 2014*x^859 + 5941*x^858 + 2805*x^857 + 2778*x^856 + 1991*x^855 + 3043*x^854 + 5530*x^853 + 5003*x^852 + 5677*x^851 + 3745*x^850 + 3204*x^849 + 593*x^848 + 3231*x^847 + 4860*x^846 + 4143*x^845 + 5347*x^844 + 291*x^843 + 4377*x^842 + 407*x^841 + 5167*x^840 + 5406*x^839 + 2561*x^838 + 322*x^837 + 540*x^836 + 4046*x^835 + 897*x^834 + 4176*x^833 + 3530*x^832 + 3649*x^831 + 1042*x^830 + 1775*x^829 + 3655*x^828 + 4928*x^827 + 1819*x^826 + 340*x^825 + 2561*x^824 + 4554*x^823 + 3386*x^822 + 4974*x^821 + 5642*x^820 + 2469*x^819 + 375*x^818 + 3653*x^817 + 1242*x^816 + 5745*x^815 + 4809*x^814 + 3789*x^813 + 1073*x^812 + 3514*x^811 + 5443*x^810 + 562*x^809 + 653*x^808 + 4162*x^807 + 5582*x^806 + 1727*x^805 + 2192*x^804 + 5534*x^803 + 2237*x^802 + 2979*x^801 + 1101*x^800 + 16*x^799 + 1470*x^798 + 3755*x^797 + 1423*x^796 + 3714*x^795 + 2623*x^794 + 1662*x^793 + 5424*x^792 + 3004*x^791 + 2826*x^790 + 654*x^789 + 3206*x^788 + 631*x^787 + 1835*x^786 + 4030*x^785 + 1767*x^784 + 2839*x^783 + 3481*x^782 + 6021*x^781 + 1162*x^780 + 2324*x^779 + 1212*x^778 + 5357*x^777 + 5418*x^776 + 5483*x^775 + 2035*x^774 + 4388*x^773 + 5911*x^772 + 2175*x^771 + 5108*x^770 + 4178*x^769 + 5058*x^768 + 3370*x^767 + 5888*x^766 + 46*x^765 + 5798*x^764 + 96*x^763 + 5844*x^762 + 2470*x^761 + 2683*x^760 + 3748*x^759 + 4657*x^758 + 4117*x^757 + 1268*x^756 + 3639*x^755 + 1266*x^754 + 3843*x^753 + 5686*x^752 + 4774*x^751 + 5648*x^750 + 5238*x^749 + 4118*x^748 + 2928*x^747 + 2136*x^746 + 2096*x^745 + 4053*x^744 + 99*x^743 + 1511*x^742 + 2811*x^741 + 2810*x^740 + 4537*x^739 + 4973*x^738 + 4304*x^737 + 4781*x^736 + 5668*x^735 + 4400*x^734 + 5508*x^733 + 5628*x^732 + 4509*x^731 + 2765*x^730 + 3769*x^729 + 4640*x^728 + 597*x^727 + 5404*x^726 + 2658*x^725 + 1416*x^724 + 5399*x^723 + 3204*x^722 + 140*x^721 + 1367*x^720 + 274*x^719 + 1116*x^718 + 5693*x^717 + 3433*x^716 + 5092*x^715 + 1355*x^714 + 400*x^713 + 4773*x^712 + 419*x^711 + 1118*x^710 + 2171*x^709 + 3830*x^708 + 5185*x^707 + 107*x^706 + 5520*x^705 + 1946*x^704 + 5416*x^703 + 3948*x^702 + 4198*x^701 + 5388*x^700 + 4388*x^699 + 4756*x^698 + 1214*x^697 + 1102*x^696 + 2492*x^695 + 3200*x^694 + 44*x^693 + 4515*x^692 + 2903*x^691 + 1413*x^690 + 5243*x^689 + 4304*x^688 + 360*x^687 + 4823*x^686 + 3671*x^685 + 1061*x^684 + 3058*x^683 + 2107*x^682 + 838*x^681 + 221*x^680 + 2534*x^679 + 2345*x^678 + 156*x^677 + 57*x^676 + 3752*x^675 + 4121*x^674 + 3597*x^673 + 2245*x^672 + 4372*x^671 + 2533*x^670 + 2355*x^669 + 2093*x^668 + 2736*x^667 + 5821*x^666 + 4386*x^665 + 1275*x^664 + 148*x^663 + 5478*x^662 + 2353*x^661 + 4406*x^660 + 4686*x^659 + 3724*x^658 + 5517*x^657 + 3197*x^656 + 4888*x^655 + 3048*x^654 + 2150*x^653 + 2715*x^652 + 5077*x^651 + 5474*x^650 + 1872*x^649 + 2143*x^648 + 956*x^647 + 3702*x^646 + 5892*x^645 + 4133*x^644 + 4547*x^643 + 158*x^642 + 2340*x^641 + 2020*x^640 + 5162*x^639 + 4771*x^638 + 4548*x^637 + 806*x^636 + 2895*x^635 + 3019*x^634 + 3191*x^633 + 2745*x^632 + 881*x^631 + 2144*x^630 + 2504*x^629 + 4499*x^628 + 916*x^627 + 2069*x^626 + 669*x^625 + 3154*x^624 + 2311*x^623 + 5823*x^622 + 1138*x^621 + 4246*x^620 + 3151*x^619 + 2648*x^618 + 789*x^617 + 5414*x^616 + 1225*x^615 + 1745*x^614 + 846*x^613 + 311*x^612 + 4157*x^611 + 707*x^610 + 1890*x^609 + 3233*x^608 + 2955*x^607 + 1168*x^606 + 3063*x^605 + 5579*x^604 + 5842*x^603 + 4567*x^602 + 1951*x^601 + 2084*x^600 + 2312*x^599 + 72*x^598 + 2973*x^597 + 295*x^596 + 4018*x^595 + 5755*x^594 + 1205*x^593 + 438*x^592 + 666*x^591 + 586*x^590 + 5537*x^589 + 814*x^588 + 3279*x^587 + 4154*x^586 + 3127*x^585 + 926*x^584 + 596*x^583 + 5926*x^582 + 2786*x^581 + 707*x^580 + 154*x^579 + 2428*x^578 + 4104*x^577 + 1705*x^576 + 5973*x^575 + 1106*x^574 + 3235*x^572 + 2702*x^571 + 5054*x^570 + 5305*x^569 + 5774*x^568 + 2912*x^567 + 2432*x^566 + 2177*x^565 + 1862*x^564 + 4930*x^563 + 5580*x^562 + 5613*x^561 + 4494*x^560 + 2773*x^559 + 396*x^558 + 5146*x^557 + 4318*x^556 + 1994*x^555 + 1038*x^554 + 2137*x^553 + 406*x^552 + 3461*x^551 + 4880*x^550 + 3704*x^549 + 2007*x^548 + 4181*x^547 + 4284*x^546 + 4504*x^545 + 3560*x^544 + 3259*x^543 + 5561*x^542 + 2396*x^541 + 614*x^540 + 5322*x^539 + 4134*x^538 + 851*x^537 + 4979*x^536 + 803*x^535 + 1509*x^534 + 2075*x^533 + 3897*x^532 + 3752*x^531 + 4202*x^530 + 2007*x^529 + 5586*x^528 + 1099*x^527 + 3179*x^526 + 3323*x^525 + 811*x^524 + 1735*x^523 + 3318*x^522 + 1087*x^521 + 3094*x^520 + 2611*x^519 + 5390*x^518 + 5770*x^517 + 5918*x^516 + 1785*x^515 + 4310*x^514 + 4459*x^513 + 5181*x^512 + 4524*x^511 + 5034*x^510 + 3973*x^509 + 5548*x^508 + 636*x^507 + 509*x^506 + 884*x^505 + 775*x^504 + 6014*x^503 + 322*x^502 + 2517*x^501 + 469*x^500 + 1541*x^499 + 5452*x^498 + 3419*x^497 + 3002*x^496 + 3417*x^495 + 4799*x^494 + 5560*x^493 + 5768*x^492 + 4105*x^491 + 4109*x^490 + 5941*x^489 + 4179*x^488 + 1706*x^487 + 5435*x^486 + 855*x^485 + 853*x^484 + 579*x^483 + 4590*x^482 + 4370*x^481 + 5546*x^480 + 4641*x^479 + 750*x^478 + 518*x^477 + 1751*x^476 + 1591*x^475 + 1870*x^474 + 2400*x^473 + 5933*x^472 + 531*x^471 + 55*x^470 + 1034*x^469 + 389*x^468 + 5518*x^467 + 2879*x^466 + 1762*x^465 + 322*x^464 + 2565*x^463 + 1610*x^462 + 2430*x^461 + 4081*x^460 + 2281*x^459 + 1005*x^458 + 4417*x^457 + 4685*x^456 + 2658*x^455 + 3658*x^454 + 2039*x^453 + 3710*x^452 + 1038*x^451 + 5906*x^450 + 925*x^449 + 5137*x^448 + 2777*x^447 + 166*x^446 + 5303*x^445 + 5827*x^444 + 2966*x^443 + 1879*x^442 + 56*x^441 + 4378*x^440 + 5051*x^439 + 5064*x^438 + 932*x^437 + 2079*x^436 + 2177*x^435 + 1357*x^434 + 4665*x^433 + 5563*x^432 + 5004*x^431 + 5064*x^430 + 3308*x^429 + 5522*x^428 + 3172*x^427 + 1638*x^426 + 5625*x^425 + 2514*x^424 + 3790*x^423 + 3147*x^422 + 2840*x^421 + 79*x^420 + 5647*x^419 + 3068*x^418 + 3685*x^417 + 1759*x^416 + 1790*x^415 + 4933*x^414 + 3564*x^413 + 1587*x^412 + 5067*x^411 + 188*x^410 + 213*x^409 + 1395*x^408 + 3192*x^407 + 4555*x^406 + 1413*x^405 + 2673*x^404 + 1586*x^403 + 1286*x^402 + 2724*x^401 + 2723*x^400 + 5792*x^399 + 1773*x^398 + 3745*x^397 + 3055*x^396 + 5323*x^395 + 3069*x^394 + 237*x^393 + 1263*x^392 + 2163*x^391 + 1646*x^390 + 3091*x^389 + 4114*x^388 + 5731*x^387 + 503*x^386 + 4060*x^385 + 3898*x^384 + 3098*x^383 + 1751*x^382 + 1840*x^381 + 1503*x^380 + 870*x^379 + 2970*x^378 + 915*x^377 + 2238*x^376 + 5730*x^375 + 3933*x^374 + 907*x^373 + 1478*x^372 + 3114*x^371 + 3787*x^370 + 474*x^369 + 5863*x^368 + 200*x^367 + 2830*x^366 + 682*x^365 + 1926*x^364 + 2*x^363 + 5832*x^362 + 3485*x^361 + 5069*x^360 + 2455*x^359 + 886*x^358 + 3690*x^357 + 5182*x^356 + 2379*x^355 + 317*x^354 + 1002*x^353 + 1507*x^352 + 4854*x^351 + 409*x^350 + 362*x^349 + 3283*x^348 + 3303*x^347 + 3336*x^346 + 1410*x^345 + 619*x^344 + 4910*x^343 + 4969*x^342 + 3095*x^341 + 3984*x^340 + 3520*x^339 + 942*x^338 + 2813*x^337 + 275*x^336 + 3951*x^335 + 4222*x^334 + 4761*x^333 + 5063*x^332 + 859*x^331 + 1321*x^330 + 4070*x^329 + 1602*x^328 + 1705*x^327 + 200*x^326 + 5885*x^325 + 5328*x^324 + 5945*x^323 + 1226*x^322 + 5824*x^321 + 2261*x^320 + 4241*x^319 + 232*x^318 + 560*x^317 + 5505*x^316 + 4141*x^315 + 2929*x^314 + 2769*x^313 + 2072*x^312 + 1755*x^311 + 5011*x^310 + 4958*x^309 + 3654*x^308 + 6005*x^307 + 3172*x^306 + 2709*x^305 + 1956*x^304 + 5284*x^303 + 4792*x^302 + 5408*x^301 + 3924*x^300 + 530*x^299 + 1570*x^298 + 549*x^297 + 2128*x^296 + 513*x^295 + 3684*x^294 + 2802*x^293 + 484*x^292 + 3349*x^291 + 357*x^290 + 1217*x^289 + 2106*x^288 + 4360*x^287 + 5137*x^286 + 1547*x^285 + 132*x^284 + 882*x^283 + 4926*x^282 + 5067*x^281 + 199*x^280 + 4020*x^279 + 5942*x^278 + 863*x^277 + 831*x^276 + 5603*x^275 + 4562*x^274 + 4200*x^273 + 2319*x^272 + 1795*x^271 + 2951*x^270 + 1343*x^269 + 4856*x^268 + 5415*x^267 + 3295*x^266 + 4609*x^265 + 5479*x^264 + 3579*x^263 + 3627*x^262 + 2388*x^261 + 4576*x^260 + 479*x^259 + 5171*x^258 + 986*x^257 + 2303*x^256 + 44*x^255 + 4674*x^254 + 4496*x^253 + 1870*x^252 + 5448*x^251 + 708*x^250 + 900*x^249 + 378*x^248 + 5332*x^247 + 3829*x^246 + 2434*x^245 + 1418*x^244 + 5656*x^243 + 4605*x^242 + 142*x^241 + 5129*x^240 + 1273*x^239 + 2534*x^238 + 2653*x^237 + 5726*x^236 + 2937*x^235 + 646*x^234 + 1763*x^233 + 4921*x^232 + 1367*x^231 + 2406*x^230 + 2641*x^229 + 2765*x^228 + 1432*x^227 + 1704*x^226 + 3155*x^225 + 5070*x^224 + 4174*x^223 + 5455*x^222 + 298*x^221 + 2494*x^220 + 1714*x^219 + 4868*x^218 + 5479*x^217 + 4636*x^216 + 5396*x^215 + 5046*x^214 + 2310*x^213 + 1265*x^212 + 2413*x^211 + 4424*x^210 + 3846*x^209 + 5016*x^208 + 793*x^207 + 3596*x^206 + 3820*x^205 + 4999*x^204 + 4088*x^203 + 1844*x^202 + 1597*x^201 + 5829*x^200 + 2520*x^199 + 1868*x^198 + 1910*x^197 + 5583*x^196 + 1931*x^195 + 957*x^194 + 5385*x^193 + 781*x^192 + 5847*x^191 + 4656*x^190 + 5717*x^189 + 532*x^188 + 2617*x^187 + 1686*x^186 + 812*x^185 + 5013*x^184 + 1206*x^183 + 281*x^182 + 4900*x^181 + 1440*x^180 + 957*x^179 + 2325*x^178 + 3318*x^177 + 2762*x^176 + 1995*x^175 + 2694*x^174 + 5048*x^173 + 143*x^172 + 2299*x^171 + 4100*x^170 + 3127*x^169 + 5807*x^168 + 5232*x^167 + 606*x^166 + 4627*x^165 + 1866*x^164 + 2212*x^163 + 3519*x^162 + 5478*x^161 + 3894*x^160 + 2609*x^159 + 498*x^158 + 4900*x^157 + 2441*x^156 + 4290*x^155 + 5793*x^154 + 2988*x^153 + 904*x^152 + 1329*x^151 + 3099*x^150 + 2932*x^149 + 3114*x^148 + 5691*x^147 + 5859*x^146 + 3637*x^145 + 4635*x^144 + 4613*x^143 + 1663*x^142 + 3276*x^141 + 5341*x^140 + 4052*x^139 + 159*x^138 + 974*x^137 + 2186*x^136 + 4269*x^135 + 5561*x^134 + 1738*x^133 + 257*x^132 + 3284*x^131 + 532*x^130 + 4893*x^129 + 1920*x^128 + 3930*x^127 + 1558*x^126 + 77*x^125 + 3693*x^124 + 860*x^123 + 4349*x^122 + 1825*x^121 + 1576*x^120 + 1526*x^119 + 2010*x^118 + 37*x^117 + 1256*x^116 + 2083*x^115 + 296*x^114 + 5211*x^113 + 927*x^112 + 4196*x^111 + 2439*x^110 + 3782*x^109 + 5873*x^108 + 5699*x^107 + 2724*x^106 + 4900*x^105 + 2305*x^104 + 277*x^103 + 105*x^102 + 3904*x^101 + 2025*x^100 + 2620*x^99 + 744*x^98 + 1806*x^97 + 4883*x^96 + 3419*x^95 + 5181*x^94 + 5306*x^93 + 6012*x^92 + 3473*x^91 + 3893*x^90 + 3403*x^89 + 1145*x^88 + 1802*x^87 + 3146*x^86 + 1804*x^85 + 5649*x^84 + 697*x^83 + 5437*x^82 + 5115*x^81 + 3678*x^80 + 4323*x^79 + 3328*x^78 + 1759*x^77 + 5254*x^76 + 4714*x^75 + 3597*x^74 + 5805*x^73 + 762*x^72 + 5682*x^71 + 3852*x^70 + 1284*x^69 + 1183*x^68 + 1174*x^67 + 3308*x^66 + 4052*x^65 + 1993*x^64 + 5807*x^63 + 388*x^62 + 5304*x^61 + 2326*x^60 + 4930*x^59 + 641*x^58 + 4361*x^57 + 3251*x^56 + 3676*x^55 + 2255*x^54 + 4979*x^53 + 4350*x^52 + 1024*x^51 + 4379*x^50 + 3819*x^49 + 5294*x^48 + 4262*x^47 + 464*x^46 + 4351*x^45 + 3337*x^44 + 5795*x^43 + 5363*x^42 + 5145*x^41 + 283*x^40 + 3783*x^39 + 5340*x^38 + 3945*x^37 + 2674*x^36 + 5033*x^35 + 6004*x^34 + 5237*x^33 + 217*x^32 + 5844*x^31 + 424*x^30 + 4067*x^29 + 3614*x^28 + 398*x^27 + 3359*x^26 + 3618*x^25 + 1842*x^24 + 6004*x^23 + 2622*x^22 + 759*x^21 + 3090*x^20 + 192*x^19 + 387*x^18 + 1394*x^17 + 1388*x^16 + 1518*x^15 + 3288*x^14 + 1417*x^13 + 2457*x^12 + 2977*x^11 + 2046*x^10 + 2290*x^9 + 1086*x^8 + 2480*x^7 + 4933*x^6 + 55*x^5 + 633*x^4 + 2253*x^3 + 3490*x^2 + 5283*x + 3019 
       
%time print mod(f(3019),6037) for i in [1..3018]: if not f(mod(i,6037)) == mod(i+1,6037): print "ons vermoeden klopte niet" 
       
1
CPU time: 94.23 s,  Wall time: 94.39 s
1
CPU time: 94.23 s,  Wall time: 94.39 s

Nu volgt dat $a_i$, opgevat als element van $\operatorname{GF}(6037)$, gelijk is aan $i\,\mathrm{mod}\,3019$. Schrijven we $100010000001000000000=3019q+r$, vinden we dat het gezochte getal (nog steeds opgevat als element van $\operatorname{GF}(6037)$) gelijk is aan $(3019!)^q\cdot r!$. Rekening houdend met $(3019!)^{6036}\equiv 1$ kunnen we nu het antwoord berekenen:

def coma3(m): G = GF(6037) q,r = 100010000001000000000.quo_rem(3019) q = q % 6036 return G(factorial(3019))^q * G(factorial(r)) coma3(100010000001000000000) 
       
3896
3896

Wiskunde is te moeilijk,
dan maar programmeren

Antwoord: 17

We zoeken eerst de maximale oneven lengte: loop alle mogelijke palindromen van lengte 3 af, en zoek deze in de lijst met cijfers van $(2012^{2012})^{2012}$. Als je er één vindt, kijk dan of de twee cijfers er links en rechts van ook gelijk zijn, indien ja vinden we een palindroom van lengte 5 en kijken we ook naar de cijfers daar rond, enzovoort. De maximale even lengte vinden we analoog (begin met alle mogelijke palindromen met lengte 2).

def coma4_0(): x=str((2012^2012)^2012); l=len(x) y=['0','1','2','3','4','5','6','7','8','9']; number=len(y) max=''; index=0 for i in [0..number-1]: a=y[i] b=y[i] pal=a+b; index=x.find(pal,0) while(index>-1): k=1 while(a==b and index-k>-1 and index+1+k<l): a=''+x[index-k] b=''+x[index+k] if(a==b): pal=a+pal+b k+=1 if(len(pal)>len(max)): max=pal a=y[i] b=y[i] pal=a+b index=x.find(a+b, index+1); tempmax=max; max=''; for i in [0..number-1]: for j in [0..number-1]: a=y[j]; m=y[i]; b=y[j]; pal=a+m+b; index=x.find(pal,0); while(index>-1): k=2 while(a==b and index+1-k>-1 and index+1+k<l): a=''+x[index+1-k]; b=''+x[index+1+k]; if(a==b): pal=a+pal+b k+=1 if(len(pal)>len(max)): max=pal a=y[j] m=y[i] b=y[j] pal=a+m+b index=x.find(pal,index+1) if(len(tempmax)>len(max)): max=tempmax return len(max) coma4_0() 
       
17
17

Alternatief kunnen we de hele lijst van cijfers van $(2012^{2012})^{2012}$ overlopen, en (analoog als hierboven) rond elk van deze cijfers een zolang mogelijk palindroom proberen bouwen. Voor de even palindromen overlopen we dan alle koppels van opeenvolgende getallen, en indien ze gelijk zijn proberen we er een langer palindroom rond te bouwen. Deze aanpak kost wel een pak meer rekentijd (een kleine 3 minuten in totaal).

def coma4_1(): A=str((2012^2012)^2012) L=len(A) onevenmax=1 for n in [0..L-1]: i=0 fout=false while n-i>0 and n+i<L-1 and fout==false: i+=1 if A[n-i]!=A[n+i]: fout=true if 2*i-1>onevenmax: onevenmax=2*i-1 evenmax=1 for n in [0..L-2]: i=0 fout=false while n-i+1>0 and n+i<L-1 and fout==false: i+=1 if A[n-i+1]!=A[n+i]: fout=true if 2*i-2>evenmax: evenmax=2*i-2 return max(onevenmax,evenmax) coma4_1() 
       
17
17

Nullen en eentjes

Antwoord: 3033576902

Zij $l_n$ de lengte van $a_n$. Men kan vinden dat $l_n$ recursief gedefiniëerd is als volgt:

\[ \begin{aligned} l_1 &= 2 \\ l_2 &= 4 \\ l_3 &= 7 \\ l_4 &= 11 \\ l_5 &= 14 \\ l_k &= \begin{cases} l_{k-1} + l_{k-3} + 1 & \text{als $k>5$ even is,}\\ l_{k-1} + l_{k-3} & \text{ als $k>5$ oneven is}\end{cases}\end{aligned} \]

%time def coma5_0(n): a=[2,4,7,11,14] while len(a)<n: a.append(a[len(a)-1]+a[len(a)-3]+1) a.append(a[len(a)-1]+a[len(a)-3]) return int(str(a[n-1])[:10]) # return floor(a[n-1]/(10^(floor(log(a[n-1],10))-9))) coma5_0(221112) 
       
3033576902
CPU time: 1.48 s,  Wall time: 1.49 s
3033576902
CPU time: 1.48 s,  Wall time: 1.49 s

Het gaat iets sneller als we niet alle tussenresultaten bijhouden:

%time # werkt als n >= 5 def coma5_1(n): (drie,vier,vijf) = (7,11,14) for i in range(int((n-5)/2)): temp = vijf+drie+1 (drie,vier,vijf) = (vijf,temp,temp+vier) if n%2 == 0: (drie,vier,vijf) = (vier,vijf,vijf+drie+1) return int(str(vijf)[:10]) coma5_1(221112) 
       
3033576902
CPU time: 0.46 s,  Wall time: 0.46 s
3033576902
CPU time: 0.46 s,  Wall time: 0.46 s
Alternatief:

Door in een $a_i$ tussen elke "$01$" een streep te zetten, verdelen we de $a_n$ in stukjes die we atomen zullen noemen. Een atoom wordt steeds omgezet in één of meerdere atomen, en waarin een atoom wordt omgezet hangt niet af wat er voor of na komt. We kunnen de atomen oplijsten; er zijn er 10:

Nr. Atoom Lengte

Wordt omgezet in nrs.

1 1 1 2
2 11 2 3,1
3 10 2 5
4 110 3 3,4
5 1110 4 6
6 11110 5 7,4
7 100 3 9
8 1100 4 3,8
9 11100 5 10
10 111100 6 7,8

We kunnen het aantal atomen van elke soort in een $a_n$ voorstellen door een vector $v_n$ met lengte 10. Dan is $v_{n+1}=Tv_{n}$, met $T$ de matrix met een $1$ op $(i,j)$ als en slechts als atoom $i$ wordt omgezet in atoom $j$. Samen met $v_2=(0,1,1,0,0,0,0,0,0,0)$ laat dit ons toe om $v_{221112}$ te berekenen, waaruit we eenvoudig de lengte van $a_{221112}$ kunnen halen.

%time def coma5_2(n): T=matrix(10,10) one=[(1,2),(2,3),(2,1),(3,5),(4,3),(4,4),(5,6),(6,7),(6,4),(7,9),(8,3),(8,8),(9,10),(10,7),(10,8)] for k in one: T[k[1]-1,k[0]-1]=1 l=(matrix([1,2,2,3,4,5,3,4,5,6])*T^(n-2)*matrix(10,1,[0,1,1,0,0,0,0,0,0,0]))[0,0] return int(str(l)[:10]) # return floor(l/(10^(floor(log(l,10))-9))) coma5_2(221112) 
       
3033576902
CPU time: 0.41 s,  Wall time: 0.41 s
3033576902
CPU time: 0.41 s,  Wall time: 0.41 s

!

Antwoord: 1032069286

Als $m=\prod_{i=1}^l{p_i^{e_i}}$, dan is $\phi(m)=n \iff \prod_{i=1}^l{(p_i-1)p_i^{e_i-1}}=n$. De mogelijke priemfactoren van $m$ zijn dus de priemgetallen $p$ met de eigenschap dat $n$ deelbaar is door $p-1$, en deze priemgetallen kunnen we gemakkelijk oplijsten (we noteren $p_1>p_2>\ldots >p_l$). We zullen nu een lijst opstellen met alle mogelijke waarden van m. Dit kunnen we recursief doen: als $p_i$ de grootste priemfactor van $m$ is, zijn de mogelijke waarden van $\frac{m}{p_i}$ bepaald door de oplossingen van $p_i^{e_i-1}\prod_{j=i+1}^l{(p_j-1)p_j^{e_j-1}}=\frac{n}{p_i-1}$. Dit aantal kunnen we dan recursief bepalen, zoals hieronder uitgewerkt. $\phi^{-1}_k{(n)}$ vinden we dan door de gevonden lijst met $m$-waarden te sorteren en het $k$-de element te nemen.

def primes6(n): f=divisors(n) primes=[] for p in f: if is_prime(int(p+1)): primes.append(p+1) return sorted(primes,reverse=true) def coma6(n): result=[] if n==0: return result if n==1: return [1,2] p=primes6(n) for i in range(len(p)): for k in rec6(p[i], n / (p[i] - 1), i,p): result.append(k) return sorted(result) def rec6(q,phi,i,p): result=[] if phi == 1 : result.append(q) if (phi % p[i]) == 0: for s in rec6(q * p[i], phi / p[i], i, p): result.append(s) for j in range(i+1,len(p)): if( phi % (p[j] - 1) == 0 ): for s in rec6(q * p[j], phi / (p[j] - 1), j, p): result.append(s) return result 
       
%time nr=[1,12,123,1234,12345] sum([coma6(factorial(8+i))[nr[i]-1] for i in [0..4]]) 
       
1032069286
CPU time: 14.04 s,  Wall time: 14.07 s
1032069286
CPU time: 14.04 s,  Wall time: 14.07 s