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=;
 long long m;
 long long num[];
 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=;i<=n;i++)
     {
         ans=(ans+((num[i]*base)%mod))%mod;
         base=(base*x)%mod;
     }
     ans=ans%mod;
     while(ans<)
         ans+=mod;
     return ans;
 }
 long long ans;
 int is()
 {
     for(long long i=0LL;i<=m;i++)
     {
         //if(i==9716LL)
             //printf("YES,%lld\n",fac(i,m*m));
         if(!fac(i,m))
         {
             //printf("fac(i,m*m):%lld\n",fac(i,m*m));
             for(long long j=0LL;j<=m;j++)
             {
                 if(!fac(i+j*m,m*m))
                 {
                     ans=i+j*m;
                     return ;
                 }
             }
             return ;
         }
     }
     return ;
 }
 int main()
 {
     #ifdef test
     freopen("test.in","r",stdin);
     #endif
     int t;
     //long long a=935124339326LL;
     //a=-a;
     //printf("%lld\n",a%(9811LL);
     scanf("%d",&t);
     for(int ttt=;ttt<=t;ttt++)
     {
         scanf("%lld",&n);
         for(int i=n;i>=;i--)
             scanf("%lld",&num[i]);
         scanf("%lld",&m);
         printf("Case #%d: ",ttt);
         if(is())
             printf("%lld\n",ans);
         else
             printf("No solution!\n");
     }
     return ;
 }