#include<bits/stdc++.h>
using namespace std;
const int N=1500,M=200,Q=6e5;
const int mod=1e9+7;
int add(int x,int y){ return (x+=y)>=mod?x-mod:x; }
int sub(int x,int y){ return (x-=y)<0?x+mod:x; }
int n,m,qs;
int b[M+5];
int p[N+5][M+5];
struct Que{
    int l,r,id;
    Que(int l=0,int r=0,int id=0):l(l),r(r),id(id){}
};
int f[N+5][M+5];
int ans[Q+5];
void dp(int f[M+5],int g[M+5],int h[M+5]){
    fill(h,h+m+1,0);
    for(int i=0;i<=m;i++)
        for(int j=0;i+j<=m;j++)
            h[i+j]=(h[i+j]+1ll*f[i]*g[j])%mod;
}
void solve(int l,int r,vector<Que> &q){
    if(l==r){
        for(auto [l,r,i]:q){
            for(int j=0;j<=m;j++)
                ans[i]=(ans[i]+1ll*p[l][j]*b[j])%mod;
        }
        return;
    }
    int mid=l+r>>1;
    for(int i=0;i<=m;i++) f[mid][i]=p[mid][i];
    for(int i=0;i<=m;i++) f[mid+1][i]=p[mid+1][i];
    for(int i=mid-1;i>=l;i--)
        dp(f[i+1],p[i],f[i]);
    for(int i=mid+2;i<=r;i++)
        dp(f[i-1],p[i],f[i]);
    static int h[M+5];
    for(int i=1;i<=mid;i++){
        memset(h,0,sizeof h);
        for(int j=0;j<=m;j++)
            for(int k=j;k<=m;k++)
                h[k-j]=(h[k-j]+1ll*f[i][j]*b[k])%mod;
        for(int j=0;j<=m;j++) f[i][j]=h[j];
    }
    vector<Que> ql,qr;
    for(auto it:q){
        auto [l,r,i]=it;
        if(r<=mid) ql.push_back(it);
        else if(l>mid) qr.push_back(it);
        else{
            int *vl=f[l],*vr=f[r];
            for(int j=0;j<=m;j++)
                ans[i]=(ans[i]+1ll*vl[j]*vr[j])%mod;
        }
    }
    solve(l,mid,ql);
    solve(mid+1,r,qr);
}
int main(){
    scanf("%d %d %d",&n,&m,&qs);
    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]);
    vector<Que> q;
    for(int i=1;i<=qs;i++){
        int l,r;
        scanf("%d %d",&l,&r);
        q.push_back({l,r,i});
    }
    solve(1,n,q);
    for(int i=1;i<=qs;i++)
        printf("%d\n",ans[i]);
    return 0;
}