HDU4569 Special equations

/*
 HDU4569 Special equations
 http://acm.hdu.edu.cn/showproblem.php?pid=4569
 数论 
 题意：f(x)为一n次方程求是否存在x, s.t. f(x)=0 (mod p^2)
 其中p为质数
 首先，我们只需考虑0-p中的x即可，因为其他的x总可以先取模
 到这个区间，因此可以在1e4内找到f(x)=0(mod p)的x
 考虑到 x=0(mod p^2) => x=0(mod p) 
 因此只需在这些x中找是否满足x=0(mod p^2)即可
 当我们找到一个满足x0=0(mod p)的x0时，
 就的到了通解x=x0+k*p;在通解中找满足题意的x即可
 */
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
#include <queue>
#include <iostream>
#include <map>
#include <set>
//#define test
using namespace std;
const int Nmax=1005;
long long m;
long long num[15];
long long n;
long long fac(long long x,long long mod)
{
    long long ans=0LL;
    long long base=1LL;
    x=x%mod;
    for(int i=0;i<=n;i++)
    {
        ans=(ans+((num[i]*base)%mod))%mod;
        base=(base*x)%mod;
    }
    ans=ans%mod;
    while(ans<0)
        ans+=mod;
    return ans;
}
long long ans;
int is()
{
    for(long long i=0LL;i<=m;i++)
    {
        //if(i==9716LL)
            //printf("YES,%lld
",fac(i,m*m));
        if(!fac(i,m))
        {
            //printf("fac(i,m*m):%lld
",fac(i,m*m));
            for(long long j=0LL;j<=m;j++)
            {
                if(!fac(i+j*m,m*m))
                {
                    ans=i+j*m;
                    return 1;
                }
            }
            return 0;
        }
    }
    return 0;
}
int main()
{
    #ifdef test
    freopen("test.in","r",stdin);
    #endif
    int t; 
    //long long a=935124339326LL;
    //a=-a;
    //printf("%lld
",a%(9811LL);
    scanf("%d",&t);
    for(int ttt=1;ttt<=t;ttt++)
    {
        scanf("%lld",&n);
        for(int i=n;i>=0;i--)
            scanf("%lld",&num[i]);
        scanf("%lld",&m);
        printf("Case #%d: ",ttt);
        if(is())
            printf("%lld
",ans);
        else
            printf("No solution!
");
    }
    return 0;
}