// Original C++ solution (slower):
#include<bits/stdc++.h>

using namespace std;

typedef long long ll;



inline ll gcd(ll a,ll b){ return b?gcd(b,a%b):a; }

inline ll lcm(ll a,ll b){ return a/gcd(a,b)*b; }



int n,q,l,r;

int p[100100];

int sum[100100], len[100100];

bool use[100100];

const ll mod = 1e9+7;

int par_sum[50][100100];



ll n_, c[1<<20];



void init(int n){

  n_ = 1;

  while(n_<n)n_*=2;



  for(int i=n_-1;i<2*n_-1;i++){

    if(i-n_+2<=n)c[i] = len[i-n_+2];

    else c[i] = 1;

  }



  for(int i=n_-2;i>=0;i--){

    c[i] = lcm(c[2*i+1],c[2*i+2]);

  }

}



ll query(int a,int b,int k,int l,int r){

  if(r<=a || b<=l)return 1;

  if(a<=l && r<=b)return c[k];

  ll vl = query(a,b,k*2+1,l,(l+r)/2);

  ll vr = query(a,b,k*2+2,(l+r)/2,r);

  return lcm(vl,vr);

}



int main(){

  scanf("%d%d",&n,&q);

  for(int i=1;i<=n;i++)scanf("%d",p+i);



  memset(use,0,sizeof(use));



  for(int i=1;i<=n;i++){

    if(!use[i]){

      int now = p[i];

      int tmplen = 1, tmpsum = i;

     

      while(now!=i){

	tmplen++; (tmpsum += now) %= mod;

	now = p[now];

      }



      now = i;

      while(!use[now]){

	use[now] = true;

	len[now] = tmplen;

	sum[now] = tmpsum;

      }

    }

  }



  for(int i=1;i<=40;i++){

    par_sum[i][0] = 0;

    for(int j=1;j<=n;j++){

      if(len[j]==i)par_sum[i][j] = (par_sum[i][j-1] + sum[j]) % mod;

      else par_sum[i][j] = par_sum[i][j-1];

    }

  }



  init(n);



  for(int i=0;i<q;i++){

    scanf("%d%d",&l,&r);

    ll cycle = query(l-1,r,0,0,n_);



    ll res = 0;

    for(int i=1;i<=40;i++){

      res += (((cycle/i)%mod)*(((par_sum[i][r]-par_sum[i][l-1])%mod+mod)%mod))%mod;

      res %= mod;

    }



    cout << res << endl;

  }

}

// Optimized solution to the same problem (faster):
#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

typedef pair<ll,ll> pll;

typedef vector<pll> vp;



const ll mod = 1e9+7;

inline ll lcm(ll a, ll b){return a/__gcd(a,b)*b;}



class SegmentTree{

  int n;

  vp node;

public:

  SegmentTree(vp a){

    int n_ = a.size();

    n=1;

    while(n<n_)n*=2;

    node.resize(2*n-1, pll(1, 0));

    for(int i=0;i<n_;i++)node[n-1+i] = a[i];

    for(int i=n-2;i>=0;i--){

      pll vl = node[2*i+1], vr = node[2*i+2];

      ll len = lcm(vl.first, vr.first);

      ll cost = ( (len/vl.first%mod) * vl.second) + ( (len/vr.first%mod) * vr.second);

      node[i] = pll(len,cost%mod);

    }

  }

  

  pll INTV_LCM(int a,int b,int k=0,int l=0,int r=0){

    if(l>=r)r = n;

    if(r<=a || b<=l)return pll(1,0);

    if(a<=l && r<=b)return node[k];

    pll vl = INTV_LCM(a,b,2*k+1,l,(l+r)/2);

    pll vr = INTV_LCM(a,b,2*k+2,(l+r)/2,r);

    ll len = lcm(vl.first, vr.first);

    ll cost = ( (len/vl.first%mod) * vl.second) + ( (len/vr.first%mod) * vr.second);

    return pll(len, cost%mod);

  }

};



int main(){

  int n,q;

  scanf("%d%d",&n,&q);

  vector<int> p(n+1);

  for(int i=1;i<=n;i++)scanf("%d",&p[i]);



  vector<pll> val(n);

  vector<int> use(n+1,0);

  for(int i=1;i<=n;i++){

    if(!use[i]){

      int cur = i;

      ll len = 0, sum = 0;

      

      do{

	len++; sum += p[cur];

	cur = p[cur];

      }while(cur != i);

      

      do{

	val[cur-1] = pll(len,sum%mod);

	use[cur] = 1;

	cur = p[cur];

      }while(cur != i);

    }

  }



  SegmentTree st(val);

  for(int i=0;i<q;i++){

    int l,r;

    scanf("%d%d",&l,&r);

    printf("%lld\n",(st.INTV_LCM(l-1,r)).second);

  }

}

// Fastest known solution for this problem: