import sys
input=lambda: sys.stdin.readline().rstrip()
n,m,k=map(int,input().split())
LR=[tuple(int(i) for i in input().split()) for _ in range(m)]
E=[[0]*n for _ in range(n)]
for l,r in LR:
  E[l-1][l-1]+=1
  if r<n:
    E[l-1][r]-=1
    E[r][l-1]-=1
    E[r][r]+=1
for i in range(n):
  for j in range(n-1):
    E[i][j+1]+=E[i][j]
for i in range(n):
  for j in range(n-1):
    E[j+1][i]+=E[j][i]

mod=10**9+7
def mul(a,b):
  c=[[0]*len(b[0]) for i in range(len(a))]
  for i in range(len(a)):
    for k in range(len(b)):
      for j in range(len(b[0])):
        c[i][j]=(c[i][j]+a[i][k]*b[k][j])%mod
  return c

def matrix_pow(a,x):
  b=[[0]*len(a) for i in range(len(a))]
  for i in range(len(a)):
    b[i][i] = 1
	
  while x>0:
    if x&1==1:
      b=mul(a,b)
    a=mul(a,a)
    x>>=1		
  return b
EK=matrix_pow(E,k)
print(EK[0][n-1])