bzoj2875

矩阵一开始写惨了。。

#include<bits/stdc++.h>
#define lowbit(a) ((a)&(-(a)))
#define clr(a,x) memset(a,x,sizeof(a))
#define rep(i,l,r) for(int i=l;i<(r);i++)
typedef long long ll;
using namespace std;
ll read()
{
    char c=getchar();
    ll ans=0,f=1;
    while(!isdigit(c)){
        if(c=='-') f=-1;
        c=getchar();
    }
    while(isdigit(c)){
        ans=ans*10+c-'0';
        c=getchar();
    }
    return ans*f;
}
ll mul(ll a,ll b,ll mod){
    ll c=0;
    while(b){
        if(b&1) c+=a;
        if(c>=mod) c-=mod;
        a<<=1;
        if(a>=mod) a-=mod;
        b>>=1;
    }
    return c;
}
ll m,a,c,x,n,g;
struct matrix{
    ll a[2][2];
    inline void clear(){
        clr(a,0);
    }
    inline matrix operator*(const matrix&A){
        matrix ans;
        ans.clear();
        rep(i,0,2){
            rep(j,0,2){
                rep(k,0,2){
                    ans.a[i][j]=(ans.a[i][j]+mul(a[i][k],A.a[k][j],m))%m;
                }
            }
        }
        return ans;
    }
    inline matrix operator=(const matrix&A){
        memcpy(a,A.a,sizeof(a));
        return *this;
    }
    inline matrix operator^(ll k){
        matrix ans,A=*this;
        ans.clear();
        ans.a[0][0]=ans.a[1][1]=1;
        while(k){
            if(k&1) ans=ans*A;
            A=A*A;
            k>>=1;
        }
        return ans;
    }
};
matrix A;
int main()
{    
    m=read(),a=read(),c=read(),x=read(),n=read(),g=read();
    A.a[0][0]=a;A.a[0][1]=0;A.a[1][0]=c;A.a[1][1]=1;
    matrix ans=A^n;
    cout<<(mul(x,ans.a[0][0],m)+ans.a[1][0])%m%g<<endl;
    return 0;
}

2875: [Noi2012]随机数生成器

Time Limit: 10 Sec  Memory Limit: 512 MB
Submit: 1290  Solved: 732
[Submit][Status][Discuss]

Description

 

Input

包含6个用空格分割的m,a,c,X0,n和g，其中a,c,X0是非负整数，m,n,g是正整数。

Output

输出一个数，即Xn mod g

Sample Input

11 8 7 1 5 3

Sample Output

2

HINT

 

Source

 

[Submit][Status][Discuss]