1654 方程的解 - Wikioi

题目描述 Description
佳佳碰到了一个难题，请你来帮忙解决。对于不定方程a1+a2+… +ak-1 +ak=g(x)，其中k≥2且k ∈ N*，x是正整数，g(x) =xx mod 1000（即xx除以1000的余数），x，k是给定的数。我们要求的是这个不定方程的正整数解组数。举例来说，当k=3， x=2时，分别为(a1，a2,a3)＝（2，1，1），（1，2，1），（1，1，2)。
输入描述 Input Description
输人只有一行，为用空格隔开的两个正整数，依次为k，x。
输出描述 Output Description
输出只有一行，为方程的正整数解组数。
样例输入 Sample Input
3 2
样例输出 Sample Output
3
数据范围及提示 Data Size & Hint
【数据范围】
对于40%的数据，ans ≤ 1016;
对于100%的数据，k≤ 100,x≤231一1，k ≤g (x)。

傻逼dp题，只是要高精度（因为内存，我高精度压了十多位）

 const
     maxn=;
     h=;
 type
     big=array[..]of int64;
 var
     f,s:array[..maxn,..maxn]of big;
     n,k:int64;

 function q(x,y:int64):int64;
 begin
     if y= then exit();
     q:=q(x,y>>);
     q:=q*q mod ;
     if y and = then q:=q*x mod ;
 end;

 operator +(a,b:big)c:big;
 var
     i:longint;
 begin
     c:=a;
     for i:= to b[] do
       inc(c[i],b[i]);
     if c[]<b[] then c[]:=b[];
     for i:= to c[]- do
       begin
         inc(c[i+],c[i]div h);
         c[i]:=c[i]mod h;
       end;
     i:=c[];
     while c[i]>=h do
       begin
         c[i+]:=c[i]div h;
         c[i]:=c[i]mod h;
         inc(c[]);
       end;
 end;

 procedure print(a:big);
 var
     i:longint;
     k:int64;
 begin
     write(a[a[]]);
     for i:=a[]- downto  do
       begin
         k:=h;
         while k> do
           begin
             k:=k div ;
             if a[i]<k then write();
           end;
         write(a[i]);
       end;
 end;

 procedure main;
 var
     i,j:longint;
 begin
     read(n,k);
     k:=q(k,k);
     for i:= to k do
       begin
         s[,i][]:=;
         s[,i][]:=;
       end;
     for i:= to n do
       begin
         for j:=k downto i do
           f[i,j]:=s[i-,j-];
         for j:=i to k do
           s[i,j]:=f[i,j]+s[i,j-];
       end;
     print(f[n,k]);
 end;

 begin
     main;
 end.