题目大意：就是给出a,b,n,m；让你求s(n);

解题思路：因为n很可能很大，所以一步一步的乘肯定会超时，我建议看代码之前，先看一下快速幂和矩阵快速幂，这样看起来就比较容易，这里我直接贴别人的推导，应该很容易懂。

看到这里你应该明白了大概吧！好吧现在继续看我的代码吧！！

AC代码：

#include
long long c[2][2],d[2];
int main()
{
    long long a,b,n,m,x,y,p,q;
    while(scanf("%I64d%I64d%I64d%I64d",&a,&b,&n,&m)!=EOF)
    {
        c[0][0]=2*a;
        c[0][1]=b-a*a;
        c[1][0]=1;
        c[1][1]=0;
        d[0]=2*a;
        d[1]=2;
        while(n)//矩阵快速幂的实现过程
        {
            if(n&1)
            {
                x=(c[0][0]*d[0]+c[0][1]*d[1])%m;
                y=(c[1][0]*d[0]+c[1][1]*d[1])%m;
                d[0]=x;
                d[1]=y;
            }
            n=n/2;
            x=(c[0][0]*c[0][0]+c[0][1]*c[1][0])%m;
            y=(c[0][0]*c[0][1]+c[0][1]*c[1][1])%m;
            p=(c[1][0]*c[0][0]+c[1][1]*c[1][0])%m;
            q=(c[1][0]*c[0][1]+c[1][1]*c[1][1])%m;
            c[0][0]=x;
            c[0][1]=y;
            c[1][0]=p;
            c[1][1]=q;
        }
        printf("%I64d\n",(d[1]+m)%m);
    }
    return 0;
}