/*DavitMarg*/
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <vector>
#include <string>
#include <cstring>
#include <list>
#include <map>
#include <unordered_map>
#include <set>
#include <unordered_set>
#include <queue>
#include <iomanip>
#include <bitset>
#include <stack>
#include <cassert>
#include <iterator>
#include <fstream>
#define mod 1000000007ll
#define LL long long
#define LD long double
#define MP make_pair
#define PB push_back
#define all(v) v.begin(), v.end()
using namespace std;
#ifndef death
#include "horses.h"
#endif
const int N=500005;
LL gcd(LL a,LL b)
{
if(!a || !b)
return a+b;
return gcd(b,a%b);
}
LL binpow(LL a,LL n)
{
if(n==0)
return 1;
if(n%2==0)
return binpow((a*a)%mod,n/2);
return (a*binpow(a,n-1))%mod;
}
LL inv(LL a)
{
return binpow(a,mod-2);
}
int n,x[N],y[N];
LL t[4*N],m[4*N],d[4*N];
void build(int v,int l,int r)
{
t[v]=d[v]=m[v]=1;
if(l==r)
return;
int mid=(l+r)/2;
build(v*2,l,mid);
build(v*2+1,mid+1,r);
}
void push(int v,int l,int r,bool f=1)
{
d[v]%=mod;
m[v]%=mod;
//LL g=gcd(d[v],m[v]);
//d[v]/=g;
//m[v]/=g;
if(d[v]==m[v])
return;
if(l!=r)
{
d[v*2]*=d[v];
d[v*2+1]*=d[v];
m[v*2]*=m[v];
m[v*2+1]*=m[v];
if(f)
{
int mid=(l+r)/2;
push(v*2,l,mid,0);
push(v*2+1,mid+1,r,0);
}
}
t[v]*=d[v];
t[v]%=mod;
t[v]*=inv(m[v]);
d[v]=m[v]=1;
}
void upd(int v,int l,int r,int i,int j,int dv,int mv)
{
push(v,l,r);
if(i>j)
return;
if(l==i && r==j)
{
d[v]*=dv;
m[v]*=mv;
push(v,l,r);
return;
}
int mid=(l+r)/2;
upd(v*2,l,mid,i,min(j,mid),dv,mv);
upd(v*2+1,mid+1,r,max(i,mid+1),j,dv,mv);
t[v]=max(t[v*2],t[v*2+1]);
}
LL get(int v,int l,int r,int i,int j)
{
push(v,l,r);
if(i>j)
return 0;
if(l==i && r==j)
return t[v];
int mid=(l+r)/2;
return max(
get(v*2,l,mid,i,min(j,mid)),
get(v*2+1,mid+1,r,max(i,mid+1),j));
}
int updateX(int pos,int val)
{
upd(1,0,n-1,pos,n-1,val,x[pos]);
x[pos]=val;
return get(1,0,n-1,0,n-1)%mod;
}
int updateY(int pos,int val)
{
upd(1,0,n-1,pos,pos,val,y[pos]);
y[pos]=val;
return get(1,0,n-1,0,n-1)%mod;
}
int init(int N,int X[],int Y[])
{
n=N;
for(int i=0;i<n;i++)
{
x[i]=X[i];
y[i]=Y[i];
}
build(1,0,n-1);
for(int i=0;i<n;i++)
{
upd(1,0,n-1,i,n-1,x[i],1);
upd(1,0,n-1,i,i,y[i],1);
}
return get(1,0,n-1,0,n-1)%mod;
}
#ifdef death
int main()
{
int N,X[102],Y[102],M;
cin>>N;
for(int i=0;i<N;i++)
cin>>X[i];
for(int i=0;i<N;i++)
cin>>Y[i];
cout<<init(N,X,Y)<<endl;
cin>>M;
while(M--)
{
int TYPE,POS,VAL;
cin>>TYPE>>POS>>VAL;
if(TYPE==1)
cout<<updateX(POS,VAL)<<endl;
else
cout<<updateY(POS,VAL)<<endl;
}
return 0;
}
#endif
/*
2
2 1
3 4
0
*/
