bc 80

可貴打一場bc，FST了兩道，顯然我還不在最佳狀態，可是我已經感受到提升了。

這一場比賽開始一個小時後纔看的題，當時只剩40分鐘，不過因爲我急於漲分，因此決定看題，先看C題，很簡單的矩陣快速冪+費馬小定理降冪，可是調了很久。。。確實不在狀態。。。不過仍是寫完了，，，而後剩下十分鐘順手過了B和A題。最後掛了兩道，A題是本身沒考慮周到，C題是以前用費馬小定理的時候一直沒注意不互質的時候，也就是a%p==0的時候。而後本身確實也沒認真去對待比賽，最後的hack沒提早造數據，隨手hack，固然hack失敗，提交代碼以前也沒考慮極端狀況，順手一交，也難怪FST。

D題也是能夠作的，就是個拓展中國剩餘定理。

C題：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#define REP(i,a,b) for(int i=a;i<=b;i++)
#define MS0(a) memset(a,0,sizeof(a))

using namespace std;

typedef long long ll;
const int maxn=1000100;
const int INF=1e9+10;

ll n,a,b,c,p;
ll p1;

struct Matrix
{
    ll a[3][3];
    friend Matrix operator*(Matrix A,Matrix B)
    {
        Matrix C;MS0(C.a);
        REP(i,0,2) REP(j,0,2) REP(k,0,2) C.a[i][j]=(C.a[i][j]+p1*3+A.a[i][k]*B.a[k][j])%p1;
        return C;
    }
};

Matrix qpow(Matrix n,ll k)
{
    Matrix res;MS0(res.a);
    REP(i,0,2) res.a[i][i]=1;
    while(k){
        if(k&1) res=res*n;
        n=n*n;
        k>>=1;
    }
    return res;
}

ll g(ll n)
{
    if(n==1) return 0;
    if(n==2) return b;
    Matrix tmp;MS0(tmp.a);
    tmp.a[0][0]=(c+1)%p1;
    tmp.a[0][1]=(p1*3+1-c)%p1;
    tmp.a[0][2]=(p1*3-1)%p1;
    tmp.a[1][0]=1;
    tmp.a[2][1]=1;
    tmp=qpow(tmp,n-3);
    ll a[3]={(b*c+b)%p1,b,0};
    ll res=0;
    REP(i,0,2) res=(res+3*p1+tmp.a[0][i]*a[i])%p1;
    return res;
}

ll qpow2(ll n,ll k,ll p)
{
    ll res=1;
    while(k){
        if(k&1) res=(res*n)%p;
        n=(n*n)%p;
        k>>=1;
    }
    return res;
}

ll solve()
{
    if(a%p==0) return 0;
    ll x=g(n);
    ll res=qpow2(a,x,p);
    return res;
}

int main()
{
    //freopen("in.txt","r",stdin);
    int T;cin>>T;
    while(T--){
        scanf("%I64d%I64d%I64d%I64d%I64d",&n,&a,&b,&c,&p);
        p1=p-1;
        printf("%I64d\n",solve());
    }
    return 0;
}

View Code

D題：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<vector>
#define REP(i,a,b) for(int i=a;i<=b;i++)
#define MS0(a) memset(a,0,sizeof(a))

using namespace std;

typedef  long long ll;
const int maxn=1000100;
const int INF=1e9+10;

int n;
int a[maxn],b[maxn];
int pos[maxn];
ll Y[maxn],Z[maxn];

void exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
    if(b==0){
        x=1;y=0;d=a;
        return;
    }
    exgcd(b,a%b,d,y,x);
    y-=a/b*x;
}

ll MLE(ll a,ll b,ll n)
{
    ll x,y,d;
    exgcd(a,n,d,x,y);
    if(b%d) return -1;
    x*=b/d;
    return (x%n+n)%n;
}

ll LCM(ll a,ll b)
{
    return a/__gcd(a,b)*b;
}

ll ex_china(ll *a,ll *m,int n)
{
    ll A=0,M=1;
    REP(i,1,n){
        ll k=MLE(m[i],A-a[i],M);
        if(k==-1) return -1;
        A=k*m[i]+a[i];
        M=LCM(M,m[i]);
    }
    if(A%M==0) return M;
    return (A+M)%M;
}

int main()
{
    freopen("in.txt","r",stdin);
    int T;cin>>T;
    while(T--){
        scanf("%d",&n);
        REP(i,1,n) scanf("%d",&a[i]),b[a[i]]=i;
        int s=0;
        int cnt=n;
        vector<int> v;
        REP(i,1,n) v.push_back(i);
        REP(i,1,n){
            for(int j=0;j<v.size();j++){
                if(v[j]==b[i]){
                    Y[i]=j+1-s;
                    Z[i]=(int)v.size();
                    vector<int>::iterator it=(v.begin()+j);
                    v.erase(it);
                    s=j;
                    break;
                }
            }
        }
        ll ans=ex_china(Y,Z,n);
        if(ans==-1) puts("Creation August is a SB!");
        else printf("%I64d\n",ans);
    }
    return 0;
}

View Code