UVA - 10870 UVA - 10870

Problem A
Recurrences
Input: standard input
Output: standard output

Consider recurrent functions of the following form:

f(n) = a₁ f(n - 1) + a₂ f(n - 2) + a₃ f(n - 3) + ... + a_d f(n - d), for n > d.
a₁, a₂, ..., a_d - arbitrary constants.

A famous example is the Fibonacci sequence, defined as: f(1) = 1, f(2) = 1, f(n) = f(n - 1) + f(n - 2). Here d = 2, a₁ = 1, a₂ = 1.

Every such function is completely described by specifying d (which is called the order of recurrence), values of d coefficients: a₁, a₂, ..., a_d, and values of f(1), f(2), ..., f(d). You'll be given these numbers, and two integers n and m. Your program's job is to compute f(n) modulo m.

Input

Input file contains several test cases. Each test case begins with three integers: d, n, m, followed by two sets of d non-negative integers. The first set contains coefficients: a₁, a₂, ..., a_d. The second set gives values of f(1), f(2), ..., f(d).

You can assume that: 1 <= d <= 15, 1 <= n <= 2³¹ - 1, 1 <= m <= 46340. All numbers in the input will fit in signed 32-bit integer.

Input is terminated by line containing three zeroes instead of d, n, m. Two consecutive test cases are separated by a blank line.

Output

For each test case, print the value of f(n) (mod m) on a separate line. It must be a non-negative integer, less than m.

Sample Input                              Output for Sample Input

1 1 100 2 1   2 10 100 1 1 1 1   3 2147483647 12345 12345678 0 12345 1 2 3 0 0 0

1 55 423

sl: 比较裸的快速幂。 拿来练练手，主要是代码。

 1 // by caonima
 2 // hehe
 3 #include <cstdio>
 4 #include <cstring>
 5 #include <algorithm>
 6 #include <vector>
 7 using namespace std;
 8 const int MAX= 1e5+;
 9 typedef long long LL;
 typedef vector<LL> vec;
 typedef vector<vec> mat;
 LL a[MAX],f[MAX];
 LL d,n,m;
 
 mat mul(mat &A,mat &B) {
     mat C(A.size(),vec(B[].size(),));
     for(int i=;i<A.size();i++) {
         for(int k=;k<B.size();k++) {
             for(int j=;j<B[].size();j++) {
                 C[i][j]=(C[i][j]+A[i][k]*B[k][j])%m;
             }
         }
     }
     return C;
 }
 
 mat pow(mat A, LL n) {
     mat B(A.size(),vec(A.size(),));
     for(int i=;i<A.size();i++) B[i][i]=;
     while(n>) {
         if(n&) B=mul(A,B);
         A=mul(A,A);
         n>>=;
     }
     return B;
 }
 
 int main() {
     while(scanf("%lld %lld %lld",&d,&n,&m)==&&(d||n||m)) {
         for(int i=;i<d;i++) {
             scanf("%lld",&a[i]);
         }
         for(int i=;i<d;i++) {
             scanf("%lld",&f[i]);
         }
         mat A(d,vec(d,));
         vec B(d,);
         for(int i=;i<d;i++) A[][i]=a[i];
         for(int i=;i<d;i++) A[i][i-]=;
         for(int i=;i<d;i++) B[i]=f[i];
         A=pow(A,n-d);
         LL ans=;
         for(int i=;i<d;i++) {
             ans=(ans+(LL)A[][i]*B[d-i-])%m;
         }
         printf("%lld\n",ans);
     }
     return ;
 }