棟棟最近迷上了隨機演算法，而隨機數是生成隨機演算法的基礎。棟棟準備使用線性同餘法（linear congruential method）來生成乙個隨機數列，這種方法需要設定四個非負整數引數m,a,c,x[0],按照下面的公式生成出一系列隨機數： x[
n+1]
=(ax
[n]+
c)mo
dm其中mod m表示前面的數除以m的餘數。從這個式子可以看出，這個序列的下乙個數總是由上乙個數生成的。
用這種方法生成的序列具有隨機序列的性質，因此這種方法被廣泛地使用，包括常用的c++和pascal的產生隨機數的庫函式使用的也是這種方法。
棟棟知道這樣產生的序列具有良好的隨機性，不過心急的他仍然想盡快知道x[n]是多少。由於棟棟需要的隨機數是0,1,…,g-1之間的，他需要將x[n]除以g取餘得到他想要的數，即x[n] mod g，你只需要告訴棟棟他想要的數x[n] mod g是多少就可以了。
輸入包含6個用空格分割的整數m,a,c,x[0],n和g，其中a,c,x[0]是非負整數，m,n,g是正整數。
輸出乙個數，即x[n] mod g
11 8 7 1 5 3
【樣例說明】
計算得x[n]=x[5]=8,故(x[n] mod g) = (8 mod 3) = 2
【資料規模】
40%的資料中m為質數
30%的資料中m與a-1互質
50%的資料中n<=10^6
100%的資料中n<=10^18
40%的資料m,a,c,x[0]<=10^4
85%的資料m,a,c,x[0]<=10^9
100%的資料中m,a,c,x[0]<=10^18
100%的資料中g<=10^8
對於所有資料,n>=1,m>=1,a>=0,c>=0,x[0]>=0,g>=1。
直接矩陣快速冪
乘法要用龜速乘
沒了。。。
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using
namespace
std;
#define ll long long
inline ll read()
struct dalao
void init()
};ll m,a,c,x0,g,n;
ll ppow(ll a,ll b,ll mod)
return s;
}dalao operator *(dalao a,dalao b)
dalao pow(dalao a,ll b)
return s;
}int main()
Bzoj2875 隨機數生成器
題目描述 棟棟最近迷上了隨機演算法，而隨機數是生成隨機演算法的基礎。棟棟準備使用線性同餘法 linear congruential method 來生成乙個隨機數列，這種方法需要設定四個非負整數引數m,a,c,x 0 按照下面的公式生成出一系列隨機數 x n 1 ax n c mod m用這種方法生...
隨機數生成器
標頭檔案 內容 rand,srand函式和rand max常量 rand max 在windows系統中為32767 在類unix系統中為2147483647 rand 函式返回乙個0 rand max的隨機整數 srand seed 函式 接受unsigned int 型別的引數seed，以see...
MATLAB隨機數生成器
1 rand 生成 0,1 區間上均勻分布的隨機數 基本語法 rand m,n,p 生成排列成m n p 多維向量的隨機數。如果只寫m，則生成m m矩陣 如果引數為 m,n 可以省略方括號。2 randn 生成服從標準正太分布 均值為0，方差為1 的隨機數 基本語法 randn m,n,p,解釋同1...