Hdu 4569 Special equations 2013长沙邀请赛

Special equations

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 20    Accepted Submission(s): 12
Special Judge

Problem Description

　　Let f(x) = a _nx ⁿ  +...+ a ₁x +a ₀, in which a _i  (0 <= i <= n) are all known integers. We call f(x) 0 (mod m) congruence equation. If m is a composite, we can factor m into powers of primes and solve every such single equation after which we merge them using the Chinese Reminder Theorem. In this problem, you are asked to solve a much simpler version of such equations, with m to be prime's square.

 

Input

　　The first line is the number of equations T, T<=50.
　　Then comes T lines, each line starts with an integer deg (1<=deg<=4), meaning that f(x)'s degree is deg. Then follows deg integers, representing a _n  to a ₀  (0 < abs(a _n) <= 100; abs(a _i) <= 10000 when deg >= 3, otherwise abs(a _i) <= 100000000, i<n). The last integer is prime pri (pri<=10000).  
　　Remember, your task is to solve f(x) 0 (mod pri*pri)

 

Output

　　For each equation f(x) 0 (mod pri*pri), first output the case number, then output anyone of x if there are many x fitting the equation, else output "No solution!"

 

Sample Input

  
  

   
   4
2 1 1 -5 7
1 5 -2995 9929
2 1 -96255532 8930 9811
4 14 5458 7754 4946 -2210 9601
  
  

 

Sample Output

  
  

   
   Case #1: No solution!
Case #2: 599
Case #3: 96255626
Case #4: No solution!
  
  

 

Source

2013 ACM-ICPC长沙赛区全国邀请赛——题目重现

 

Recommend

zhoujiaqi2010

 

 

容易看出，若x₀是同余方程

f(x) º 0 (mod p^a)                      (1)

的解，则它必是方程

f(x) º 0 (mod p^{a- 1})                     (2)

的解，因此，必有与x₀相应的方程(2)的某个解x₁，使

x₀º x₁ (mod p^{a- 1})，x₀ = x₁+p^{a- 1}t₀，

此处，t₀是某个适当的整数。

这提示我们：可以从方程(2)的解中去求方程(1)的解。于是，现在的问题是，对于方程(2)的每个解x₁，是否必有方程(1)的解x₀与之对应？若有，如何去确定它？

定理  设p是素数，a³ 2是整数，f(x) = a_nxⁿ+L+ a₁x + a₀是整系数多项式，又设x₁是同余方程(2)的一个解。以f¢(x)表示f(x)的导函数。

(ⅰ) 若f¢(x₁)0 (mod p)，则存在整数t，使得

x = x₁+p^{a- 1}t                       (3)

是同余方程(1)的解。

(ⅱ) 若f¢(x₁) º 0 (mod p)，并且f(x₁)º 0 (mod p^a)，则对于t = 0，1, 2, L, p - 1，式(3)中的x都是方程(1)的解。

证明  我们来说明，如何确定式(3)中的t，使x₁+ p^{a- 1}t满足同余方程(1)，即要使

a_n(x₁+p^{a- 1}t)ⁿ+a_{n- 1}(x₁+p^{a- 1}t)^{n- 1}+L+a₁(x₁+p^{a- 1}t)+a₀º 0 (mod p^a)， (4)

即

f(x₁) +p^{a- 1}tf¢(x₁)º 0 (mod p^a)，

tf¢(x₁)º-(mod p)。                  (5)

下面考虑两种情形。

(ⅰ)  若f¢(x)0 (mod p)，则关于t的同余方程(5)有唯一解t º t₀ (mod p)，即t = t₀+ pk（kÎZ），于是

x º x₁+ p^{a- 1}t₀ (mod p^a)

是同余方程(1)的解。

(ⅱ) 若f¢(x₁) º 0 (mod p)，并且f(x₁)º 0 (mod p^a)，则式(5)对于任意的整数t成立，即同余方程(5)有p个解

t º i(mod p)，0 £ i £ p - 1。

于是x º x₁+ p^{a- 1}i (mod p^a)，0 £ i £ p - 1，都是同余方程(1)的解。证毕。

在定理中，没有对f¢(x₁)º 0 (mod p)并且 f(x₁)0 (mod p^a)的情形进行讨论。事实上，此时，同余方程(5)无解。即，我们无法从同余方程(2)的解x₁出发去求得同余方程(1)的解。

由定理，可以把解同余方程(1)，转化为解同余方程

f(x) º 0 (mod p)。                     (6)

事实上，由方程(6)的解，利用定理，可以求出方程f(x) º 0 (mod p²)的解，再利用定理，又可以求出方程f(x)º 0 (mod p³)的解，LL，直到求出方程(1)的解。

 

 

 

#include<stdio.h>
#include<stdlib.h>
#include<vector>
using namespace std;
typedef __int64 II;
int gcd(int a,int b){
    return b==0? a:gcd(b,a%b);
}
typedef struct {
    int deg,a[15];
}function;
function f,diff;
II op(function f, II i){
    II ret=0;
    for(int loop=f.deg;loop>=0;loop--){
        ret=ret*i+(II)f.a[loop];
    }
    return ret;
}
II op(function f, II i, II mod){
    II ret=0;
    for(int loop=f.deg;loop>=0;loop--){
        
        ret=ret*i+f.a[loop];
        ret=ret%mod;
    }
    return ret;
}
II mm(II a,II b,II c){
    if(b==0) return 1%c;
    II tmp=mm(a,b/2,c);
    return tmp*tmp%c*(b%2? a:1)%c;
}
II ext_gcd(II a,II &x,II b,II &y){//a*x+b*y=(x,y)=d
    if(b==0){
        x=1;y=0;
        return a;
    }
    II xx,yy,d=ext_gcd(b,xx,a%b,yy);
    x=yy;
    y=xx-a/b*yy;
    return d;
}

vector<II> linear(II a,II b,II m){//ax=b(mod m)
    a%=m,b%=m;
    vector<II> ret;
    ret.clear();
    II x,y,d=ext_gcd(a,x,m,y);
    //if(d<0) d=-d;
    if(b%d==0){
        II e=(x*(b/d)%m+m)%m;//基础解
        //for(int i=0;i<(d<0? -d:d);i++) ret.push_back((e+i*(m/d))%m);
        ret.push_back(e);
    }
    return ret;
}
int main(){
    //freopen("in_std.txt","r",stdin);
    //freopen("out_std.txt","w",stdout);
    int T;
    II pri,p2;
    scanf("%d",&T);
    for(int cas=1;cas<=T;cas++){
        scanf("%d",&f.deg);
        for(int i=f.deg;i>=0;i--) scanf("%d",&f.a[i]);
        scanf("%I64d",&pri);
        p2=pri*pri;
        II i;
        for(i=0;i<pri;i++)
            if(op(f,i,pri)==0) break;
        if(i>=pri) {printf("Case #%d: No solution!\n",cas);continue;}
        diff.deg=f.deg-1;
        for(int loop=f.deg-1;loop>=0;loop--){
            diff.a[loop]=f.a[loop+1]*(loop+1);//?
        }
        II fx1=op(diff,i),t;
        if(fx1%pri!=0){//fx1与pri互质 
            II fx=op(f,i);
            vector<II> fuck=linear(fx1,-fx/pri,pri);
            t=fuck[0];
            printf("Case #%d: %I64d\n",cas,((i+t*pri)%p2+p2)%p2);
        }else{
            II fx=op(f,i);
            if(fx%p2==0) printf("Case #%d: %I64d\n",cas,i);
            else printf("Case #%d: No solution!\n",cas);
        }        
    }
    //while(1);
    return 0;   
}