HDU5391威爾遜定理

威爾遜定理 當且僅當p為素數,p | (p-1)!+1

若p為合數,則p=a*b;如果a!=b,那么p|(p-1)!,
如果a=b,如果p為4,那么p|(p-1)!=2,如果p大于4,那么sqrt§和sqrt(2q)肯定屬于(p-1)!中,可以整除

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<climits>
#include<cctype>
#include<queue>
#include<set>
#include<ctime>using namespace std;typedef long long ll;
const int INF=0x3f3f3f3f;
const int MAXN=1e5+5;ll mult(ll a, ll b, ll p) {a %= p;b %= p;ll res = 0, tmp = a;while(b) {if (b & 1) {res += tmp;if (res > p) res -= p;}tmp <<= 1;if (tmp > p) tmp -= p;b >>= 1;}return res;
}ll quick_pow(ll a, ll b, ll p) {ll res = 1;ll tmp = a % p;while(b) {if (b & 1) res = mult(res, tmp, p);tmp = mult(tmp, tmp, p);b >>= 1;}return res;
}bool check(ll a, ll n, ll x, ll t) {ll res = quick_pow(a, x, n);ll last = res;for (ll i = 1; i <= t; i++) {res = mult(res, res, n);if (res == 1 && last != 1 && last != n-1) return true;last = res;}return res != 1;
}bool Miller_Rabin(ll n) {if (n < 2) return false;if (n == 2) return true;if ((n & 1) == 0) return false;ll x = n-1;ll t = 0;while ((x & 1) == 0) {x >>= 1;t++;}srand(time(NULL));const ll tims = 8;for (ll i = 0; i < tims; i++) {ll a = rand() % (n-1) + 1;if (check(a, n, x, t)) return false;}return true;
}int main()
{int T,n;scanf("%d",&T);while(T--){scanf("%d",&n);if(n==4){printf("2\n");}else if(Miller_Rabin(n)){printf("%d\n",n-1);}else{printf("0\n");}}return 0;
}

本文來自互聯網用戶投稿,該文觀點僅代表作者本人,不代表本站立場。本站僅提供信息存儲空間服務,不擁有所有權,不承擔相關法律責任。
如若轉載,請注明出處:http://www.pswp.cn/news/383924.shtml
繁體地址,請注明出處:http://hk.pswp.cn/news/383924.shtml
英文地址,請注明出處:http://en.pswp.cn/news/383924.shtml

如若內容造成侵權/違法違規/事實不符,請聯系多彩編程網進行投訴反饋email:809451989@qq.com,一經查實,立即刪除!

相關文章

C++的基本認識

簡單介紹C 語言特點 支持數據封裝和數據隱藏 在C中&#xff0c;類是支持數據封裝的工具&#xff0c;對象則是數據封裝的實現。C通過建立用戶定義類支持數據封裝和數據隱藏。 在面向對象的程序設計中&#xff0c;將數據和對該數據進行合法操作的函數封裝在一起作為一個類的定…

OD 投籃大賽

/*** 題目描述* 你現在是一場采用特殊賽制投籃大賽的記錄員。這場比賽由若干回合組成&#xff0c;過去幾回合的得分可能會影響以后幾回合的得分。* 比賽開始時&#xff0c;記錄時空白的。你會得到一個記錄操作的字符串列表aops&#xff0c;其中ops[i]是你需要記錄的第i項操作&a…

IO多路復用之epoll總結

http://www.cnblogs.com/Anker/p/3263780.html 1、基本知識 epoll是在2.6內核中提出的&#xff0c;是之前的select和poll的增強版本。相對于select和poll來說&#xff0c;epoll更加靈活&#xff0c;沒有描述符限制。epoll使用一個文件描述符管理多個描述符&#xff0c;將用戶關…

2018南京區域賽 J-Prime Game

完全沒有頭緒 聽完隊友講的我還是楞了好半天菜慢慢理解.我好菜啊 首先要弄懂題目的意思,轉換一下題意就是求每個素因子出現區間的次數.區間長度最短為1.我們分析,第一個數的因子會影響1* n個區間(暫時不考慮重復),第二個數的因子會影響2 * (n-1)個區間,以此類推.因此我們只需要…

3_V1-類和對象 -- 默認成員函數

定義一個日期類 #include <iostream> #include <assert.h> using namespace std;class Date { public:void Display(); private:int _year;int _month;int _day; }; 注意: 在定義一個類的時候往往會將其成員變量定義為私有,成員函數定義為公有.這是為了達到軟件…

C++ 類模板二(類模版與友元函數)

http://www.cnblogs.com/zhanggaofeng/p/5661829.html //類模版與友元函數 #include<iostream> using namespace std;template<typename T> class Complex{ public:Complex(T a,T b);void Print() const//const修飾的是this指針{cout << this->Real <&…

HDU - 2973威爾遜定理

核心問題就是那個等式 我們觀察到等式可以寫成(n-1)!-1/n-[(n-1)!/n]的形式&#xff0c;這樣就應該聯想到威爾遜定理了。 回顧一下威爾遜定理的內容&#xff1a;當且僅當n為素數的時候n|(n-1)!-1&#xff0c;n為合數且大于4的時候n|(n-1)!【參見威爾遜定理的證明】 對于這個…

linux網絡編程之posix 線程(四):posix 條件變量與互斥鎖 示例生產者--消費者問題

http://blog.csdn.net/jnu_simba/article/details/9129939 一、posix 條件變量 一種線程間同步的情形&#xff1a;線程A需要等某個條件成立才能繼續往下執行&#xff0c;現在這個條件不成立&#xff0c;線程A就阻塞等待&#xff0c;而線程B在執行過程中使這個條件成立了&#x…

3-V2-類和對象 -- const內聯 靜態成員 友元

const修飾成員函數 在成員函數后面加一個const, const修飾this指針指向的對象, 保證調用這個const成員函數的對象在函數內不會被改變 注意:成員函數如果不修改成員變量,則加上const,成員函數如果要修改成員變量,此時就不能給其加上const修飾了 1.const對象不能調用非const…

C語言 二級指針內存模型混合實戰

http://www.cnblogs.com/zhanggaofeng/p/5485833.html //二級指針內存模型混合實戰 #define _CRT_SECURE_NO_WARNINGS #include <stdio.h> #include <stdlib.h> #include <string.h>//將內存模型①和內存模型②的數據拷貝到內存模型③ char ** threemodel(ch…

擴展歐幾里得算法

對于a?xb?yca*xb*yca?xb?yc,這樣一個二元一次方程組&#xff0c;我們想要得到他的一組解可以用擴展歐幾里得算法&#xff0c;參數列表的a,b,x,y就是方程中的a,b,x,y&#xff0c;d計算出來是gcd(a,b)。 算法求出來的是a?xb?ygcd(a,ba*xb*ygcd(a,ba?xb?ygcd(a,b的一組解…

Linux 網絡編程八(epoll應用--大并發處理)

http://www.cnblogs.com/zhanggaofeng/p/5901316.html //頭文件 pub.h #ifndef _vsucess#define _vsucess#ifdef __cplusplus extern "C" {#endif //服務器創建socket int server_socket(int port);//設置非阻塞 int setnonblock(int st);//接收客戶端socket int ser…

約瑟夫問題

n個人編號為0…n-1圍成一個圈,從0開始報數,每經過k個人那個人就退出這個圈不再報數,問最后留下來的人的編號. 樸素的做法當然是模擬,但是n,k的值一旦變得比較大的時候就難以解決問題. 我們考慮歸納的解決問題 當只有一個人的時候答案顯然為0, 假設我們已知n-1個人的時候答案為…

【數據結構與算法】內部排序之三:堆排序(含完整源碼)

轉載請注明出處&#xff1a;http://blog.csdn.net/ns_code/article/details/20227303 前言 堆排序、快速排序、歸并排序&#xff08;下篇會寫這兩種排序算法&#xff09;的平均時間復雜度都為O&#xff08;n*logn&#xff09;。要弄清楚堆排序&#xff0c;就要先了解下二叉堆這…

模線性方程(中國剩余定理+擴展中國剩余定理)

已知一系列除數和模數,求最小的滿足條件的數 我們先考慮一般的情況&#xff0c;即模數不互質。&#xff08;擴展中國剩余定理&#xff09; 我們考慮兩個方程的情況 x%MR xk1?MRxk1 * MRxk1?MR x%mr xk2?mrxk2 * mrxk2?mr 所以k1?MRk2?mrk1 * MRk2 * mrk1?MRk2?mr 即…

C++:Vector和List的實現

Vector的實現 //test.h #pragma once#include <iostream> #include <cstdio> #include <string.h> #include <assert.h>using namespace std;typedef int DataType;#define TESTHEADER printf("\n%s\n", __FUNCTION__)class Vector { publi…

【數據結構】(面試題)使用兩個棧實現一個隊列(詳細介紹)

http://blog.csdn.net/hanjing_1995/article/details/51539578 使用兩個棧實現一個隊列 思路一&#xff1a; 我們設定s1是入棧的&#xff0c;s2是出棧的。 入隊列&#xff0c;直接壓到s1即可 出隊列&#xff0c;先把s1中的元素倒入到s2中&#xff0c;彈出s2中的棧頂元素&#x…

POJ 1006 Biorhythms

中國剩余定理的模板題 只是有一個問題就是求出來Xk*MR中的R比給定的日期還大&#xff0c;但是如果負數的整除就不是向下取整了&#xff0c;為了解決這個問題&#xff0c;我們將R減小M&#xff0c;這樣總是正的&#xff0c;求出來的就沒有什么問題。 #include <iostream>…

POJ 3696 歐拉函數+快速冪

題目的意思大概就是問是否存在一串全是8的數字是L的倍數 直接想沒有什么想法&#xff0c;要想到用簡潔的形式將這個數字表示出來&#xff0c;對于每一位都是8的數字我們可以用 X8*(10k-1)/9的形式表示出來&#xff0c;那么題目的意思就是求X使L|X&#xff0c;我們先處理一下8和…

兩個棧實現一個隊列,兩個隊列實現一個棧

http://blog.csdn.net/zw_1510/article/details/51927554 問題1&#xff1a;用兩個棧實現一個隊列&#xff0c;實現隊列的push和delete操作 棧的特性是先進后出&#xff08;FILO&#xff09;,隊列的特性是先進先出&#xff08;FIFO&#xff09;,在實現delete時&#xff0c;我們…