Cod sursa(job #927803)

Utilizator razvan.popaPopa Razvan razvan.popa Data 26 martie 2013 01:13:31
Problema Pod Scor 5
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.67 kb
#include<cstdio>
#include<algorithm>
#define INF (1 << 30)
#define FOR(i,a,b)\
   for(int i=a; i<=b; ++i)
#define FORR(i,a,b)\
   for(int i=a; i>=b; --i)
#define infile "pod.in"
#define outfile "pod.out"
#define mMax 1005
#define kMax 23
#define MOD 9901
using namespace std;

struct Matrix{
	int M[kMax][kMax];
	int L, C;

    #define A (*this)

	Matrix(int L, int C){
		A.L = L;
		A.C = C;

		FOR(i,0,L)
		   FOR(j,0,C)
		      M[i][j] = 0;
	}

	Matrix operator * (const Matrix &B)const{
		Matrix C(A.L, B.C);

		FOR(i, 1, A.L)
		   FOR(j, 1, B.C)
		      FOR(k, 1, A.C){
		         C.M[i][j] += A.M[i][k] * B.M[k][j];

		         if(C.M[i][j] > MOD)
		        	 C.M[i][j] %= MOD;
		}

		return C;
	}
} I(kMax, kMax), EXP(kMax, kMax);

int V[mMax];

int N, M, K, Sol;

void read(){
    freopen(infile, "r", stdin);
    
    scanf("%d %d %d", &N, &M, &K);
    
    FOR(i,1,M)
       scanf("%d", &V[i]);

    fclose(stdin);
}

Matrix pow(int K){
	if(!K)
		return I;

	Matrix r = pow(K>>1);
	r = r * r;

	if(K & 1)
		r = r * EXP;

	return r;
}

void solve(){
    sort(V+1, V+M+1);

    I.L = I.C = K;
    FOR(i,1,K)
       I.M[i][i] = 1;

    EXP.L = EXP.C = K;
    FOR(i,1,K)
        EXP.M[i+1][i] = 1, EXP.M[i][K] = 1;

    Matrix DP(1,K);
    DP.M[1][K] = 1;

    V[++ M] = N;

    FOR(i,1,M){
        DP = DP * pow(V[i] - V[i - 1]);

        if (i != M)
        	DP.M[1][K] = 0;
    }
    
    Sol = DP.M[1][K];
}

void print(){
    freopen(outfile, "w", stdout);
    
    printf("%d\n", Sol);
    
    fclose(stdout);
}

int main(){
    read();
    solve();
    print();

    return 0;
}