#include <fstream>
using namespace std;

const long long N=21,K=31,S=15001,mod=3210121;
long long fact[S],invfact[S],s,k,n,rez;

ifstream in("cowfood.in");
ofstream out("cowfood.out");

inline long long calc(long long a,long long b)
{
    return (int)((long long)a*b%mod);
}

inline long long comb(long long n,long long k)
{
    return calc(fact[n],calc(invfact[k],invfact[n-k]));
}

struct Try
{
    long long a[K];
    long long val()
    {
        long long sum=0;
        for (long long i=1;i<=k;i++)
            sum+=a[i];
        return comb(s-sum+k,k);
    }
    void merge(Try x)
    {
        for (long long i=1;i<=k;i++)
            a[i]=max(a[i],x.a[i]);
    }
    void init()
    {
        for (long long i=0;i<=k;i++)
            a[i]=0;
    }
    void get()
    {
        for (long long i=1;i<=k;i++)
            in>>a[i];
    }
} a[N],zero;

void euclid(long long a,long long b,long long &x,long long &y)
{
    if (!b)
    {
        x=1;
        y=0;
        return;
    }
    long long x1,y1;
    euclid(b,a%b,x1,y1);
    x=y1;
    y=x1-a/b*y1;
}

long long inv(long long x)
{
    long long a,b;
    euclid(x,mod,a,b);
    a=(a%mod+mod)%mod;
    return a;
}

void calc_comb(long long n)
{
    fact[0]=invfact[0]=1;
    for (long long i=1;i<=n;i++)
        fact[i]=calc(i,fact[i-1]);
    invfact[n]=inv(fact[n]);
    for (long long i=n-1;i;i--)
        invfact[i]=calc(invfact[i+1],i+1);
}

void pie(long long p,long long semn,Try lim)
{
    if (p==n+1)
    {
        rez=(rez+semn*lim.val())%mod;
        return;
    }
    pie(p+1,semn,lim);
    lim.merge(a[p]);
    pie(p+1,-semn,lim);
}

int main()
{
    calc_comb(S-1);
    in>>k>>s>>n;
    for (long long i=1;i<=n;i++)
        a[i].get();
    rez-=s*k+1;
    pie(1,1,zero);
    while (rez<0)
        rez+=mod;
    out<<rez<<"\n";
    return 0;
}