长沙邀请赛A&E

都是独立完成的，而且都比较顺利的敲了出来。。

A题，详见组合数学中关于用齐次方程法表达递推公式那一节，直接上代码了：

#include <iostream>
#include <cstdio>

#define cl(a) memset(a,0,sizeof(a))
#define ss(a) scanf("%I64d",&a)
#define ll __int64
#define N 35

using namespace std;

ll g[N][3][3],f[3],mod;

void init(ll a,ll b)
{
    cl(g);cl(f);
    g[0][1][1]=2*a;g[0][1][2]=mod+(b-a*a)%mod;
    g[0][2][1]=1;g[0][2][2]=0;
    f[0]=2;f[1]=2*a;
}

void mul(int x)
{
    int i,j,k;
    for (i=1;i<=2;i++)
        for (j=1;j<=2;j++)
        {
            g[x][i][j]=0;
            for (k=1;k<=2;k++) g[x][i][j]+=(g[x-1][i][k]*g[x-1][k][j])%mod;
            g[x][i][j]=(g[x][i][j])%mod;
        }
}

void zmul(int x)
{
    ll f1,f2;
    f1=f[1]*g[x][1][1]+f[0]*g[x][1][2];
    f2=f[1]*g[x][2][1]+f[0]*g[x][2][2];
    f[1]=f1%mod;
    f[0]=f2%mod;
}

int main()
{
    ll a,b,n;
    while (ss(a)!=EOF)
    {
        ss(b);ss(n);ss(mod);
        init(a,b);
        n--;
        int i=0;
        while (n>0)
        {
            if (i>0) mul(i);
            if (n&1) zmul(i);
            n>>=1;
            i++;
        }
        printf("%lld\n",f[1]);
    }
    return 0;
}

E题，抓住模prime^2这个细节，令前面的表达式为g(x),如果g(x)%(pri^2)=0,则g(x)%pri=0.

g(x)是不超过4次的多项式，g(x)=g(x+pri)(mod pri)，因此有循环节pri。循环节挺小的<10^4,枚举g(y)%pri=0,再通枚举y'=k*pri+y,g(y')%(pri*pri),得出最后结果。

根据多项式模循环，循环节为模数，可以圈定y'范围<pri*pri,则0<=k<pri

注意:在第一个小循环中，g(y)%pri=0,有多个解，每个解都要进行第二次枚举。附代码：

#include <iostream>
#include <cstdio>

#define N 10005
#define ss(a) scanf("%d",&a)
#define cl(a) memset(a,0,sizeof(a))
#define ll __int64

using namespace std;

int mod,flag,a[6];
ll p[6];

ll js(ll x,ll mm)
{
    ll res=0;
    p[0]=1;
    for (int i=0;i<=4;i++) 
    {
        if (i>0) p[i]=(p[i-1]*x)%mm;
        res=(res+(p[i]*a[i])%mm)%mm;
    }
    return res;
}

void deal(int k)
{
    int i;
    if (!js(k,mod))
    {
        for (i=0;i<mod;i++)
            if (!js(i*mod+k,mod*mod))
            {
                flag=0;
                printf("%d\n",i*mod+k);
                return;
            }
    }
    return;
}

int main()
{
    int T,n,i,cas;
    ss(T);
    for (cas=1;cas<=T;cas++)
    {
        cl(a);cl(p);
        ss(n);
        for (i=n;i>=0;i--) ss(a[i]);
        ss(mod);
        flag=1;
        printf("Case #%d: ",cas);
        for (i=0;i<mod&&flag;i++) deal(i);
        if (flag) printf("No solution!\n");
    }
    return 0;
}