#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#define mo 1000000007
using namespace std;
int n,m,q,in1,in2,b[210],p[1510][210],res[210];
namespace mtt
{
    const int N=1000,mod=1e9+7;
    const double pi=acos(-1);
    struct com
    {
        double r,i;
        com(){}
        com(double _r,double _i):r(_r),i(_i){}
        com operator+(const com x){return {r+x.r,i+x.i};}
        com operator-(const com x){return {r-x.r,i-x.i};}
        com operator*(const com x){return {r*x.r-i*x.i,r*x.i+i*x.r};}
        com operator/(const double x){return {r/x,i/x};}
        com conj(){return {r,-i};}
    };
    int m,r[N],l;
    int a[N],b[N],x[N],y[N];
    com w[N];
    void fft(com *a)
    {
        for(int i=0;i<m;i++)
            if(i<r[i])
                swap(a[i],a[r[i]]);
        for(int i=1;i<m;i<<=1)
        {
            for(int j=0;j<m;j+=(i<<1))
            {
                for(int k=0;k<i;k++)
                {
                    com x=a[j+k],y=w[m/(i<<1)*k]*a[j+k+i];
                    a[j+k]=x+y;
                    a[j+k+i]=x-y;
                }    
            }
        }
    }
    void cal(int *x,int *y)
    {
        com a[N],b[N],da[N],db[N],dc[N],dd[N];
        for(int i=0;i<m;i++)
        {
            a[i]=com(x[i]&32767,x[i]>>15);
            b[i]=com(y[i]&32767,y[i]>>15);
        }
        fft(a),fft(b);
        for(int i=0;i<m;i++)
        {
            int j=(m-i)&(m-1);
            com va,vb,vc,vd;
            va=(a[i]+a[j].conj())*com(0.5,0);
            vb=(a[i]-a[j].conj())*com(0,-0.5);
            vc=(b[i]+b[j].conj())*com(0.5,0);
            vd=(b[i]-b[j].conj())*com(0,-0.5);
            da[j]=va*vc;
            db[j]=va*vd;
            dc[j]=vb*vc;
            dd[j]=vb*vd;
        }
        for(int i=0;i<m;i++)
        {
            a[i]=da[i]+db[i]*com(0,1);
            b[i]=dc[i]+dd[i]*com(0,1);
        }
        fft(a),fft(b);
        for(int i=0;i<m;i++)
        {
            long long va=a[i].r/m+0.5,vb=a[i].i/m+0.5,vc=b[i].r/m+0.5,vd=b[i].i/m+0.5;
            va%=mod,vb%=mod,vc%=mod,vd%=mod;
            x[i]=(va+((vb+vc)<<15)+((vd)<<30))%mod;
        }
    }
    vector<int>mul(vector<int>a,vector<int>b)
    {
        m=1,l=0;
        int ca=a.size(),cb=b.size();
        while(m<ca+cb-1)
            m<<=1,l++;
        for(int i=0;i<m;i++)
            r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));
        for(int i=0;i<m;i++)
            w[i]=com(cos(2*pi*i/m),sin(2*pi*i/m));
        for(int i=0;i<m;i++)
            x[i]=y[i]=0;
        for(int i=0;i<ca;i++)
            x[i]=a[i];
        for(int i=0;i<cb;i++)
            y[i]=b[i];
        cal(x,y);
        a.resize(ca+cb-1);
        for(int i=0;i<a.size();i++)
            a[i]=x[i];
        return a;
    }
}
inline int MO(int x)
{
	return x<mo?x:x-mo;
}
inline void wk(int* a,int* b,int* c)
{
	vector<int>bb,cc;
	for(int i=0;i<=m;i++)
	{
		bb.push_back(b[i]);
	}
	for(int i=0;i<=m;i++)
	{
		cc.push_back(c[i]);
	}
	vector<int>aa=mtt::mul(bb,cc);
	for(int i=0;i<=m;i++)
	{
		a[i]=aa[i];
	}
}
struct node
{
	int **dat;
	node *lch,*rch;
	node(int l,int r)
	{
		if(l==r)
		{
			dat=nullptr;
			return;
		}
		dat=new int*[r-l+1];
		for(int i=0;i<=r-l;i++)
		{
			dat[i]=new int[m+1];
		}
		int md=l+r>>1;
		memcpy(dat[md-l],p[md],sizeof(int)*(m+1));
		memcpy(dat[md+1-l],p[md+1],sizeof(int)*(m+1));
		for(int i=md-1;i>=l;i--)
		{
			wk(dat[i-l],dat[i-l+1],p[i]);
		}
		for(int i=md+2;i<=r;i++)
		{
			wk(dat[i-l],dat[i-l-1],p[i]);
		}
		lch=new node(l,md);
		rch=new node(md+1,r);
	}
	inline int wkot(int x,int y,int l,int r)
	{
		int md=l+r>>1;
		if(x==md+1)
		{
			int ans=0;
			for(int i=0;i<=m;i++)
			{
				ans=MO(ans+1ll*b[i]*dat[y-l][i]%mo);
			}
			return ans;
		}
		else if(y==md)
		{
			int ans=0;
			for(int i=0;i<=m;i++)
			{
				ans=MO(ans+1ll*b[i]*dat[x-l][i]%mo);
			}
			return ans;
		}
		else if(x<=md&&y>md)
		{
			int ans=0;
			wk(res,dat[x-l],dat[y-l]);
			for(int i=0;i<=m;i++)
			{
				ans=MO(ans+1ll*b[i]*res[i]%mo);
			}
			return ans;
		}
		else return y<md?lch->wkot(x,y,l,md):rch->wkot(x,y,md+1,r);
	}
};
node *no;
int main()
{
	// freopen("1.in","r",stdin);
	scanf("%d%d%d",&n,&m,&q);
	for(int i=0;i<=m;i++)
	{
		scanf("%d",&b[i]);
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<=m;j++)
		{
			scanf("%d",&p[i][j]);
			p[i][j]=MO(p[i][j]);
		}
	}
	no=new node(1,n);
	for(int i=1;i<=q;i++)
	{
		scanf("%d%d",&in1,&in2);
		if(n==1)
		{
			int ans=0;
			for(int i=0;i<=m;i++)
			{
				ans=MO(ans+1ll*p[1][i]*b[i]%mo);
			}
			printf("%d\n",ans);
		}
		else printf("%d\n",no->wkot(in1,in2,1,n));
	}
	return 0;
}