#pragma GCC target("avx") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include // #include // #include // #include // using namespace __gnu_pbds; // #include // namespace multiprecisioninteger = boost::multiprecision; // using cint=multiprecisioninteger::cpp_int; using namespace std; using ll=long long; // #define double long double using datas=pair; using ddatas=pair; using tdata=pair; using vec=vector; using mat=vector; using pvec=vector; using pmat=vector; // using llset=tree,rb_tree_tag,tree_order_statistics_node_update>; #define For(i,a,b) for(i=a;i<(ll)b;++i) #define bFor(i,b,a) for(i=b,--i;i>=(ll)a;--i) #define rep(i,N) For(i,0,N) #define rep1(i,N) For(i,1,N) #define brep(i,N) bFor(i,N,0) #define brep1(i,N) bFor(i,N,1) #define all(v) (v).begin(),(v).end() #define allr(v) (v).rbegin(),(v).rend() #define vsort(v) sort(all(v)) #define vrsort(v) sort(allr(v)) #define uniq(v) vsort(v);(v).erase(unique(all(v)),(v).end()) #define endl "\n" #define eb emplace_back #define print(x) cout< ostream& operator<<(ostream& os,const pair& p){return os<<"("< ostream& operator<<(ostream& os,const vector& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const set& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const map& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< inline bool chmax(T& a,T b){bool x=a inline bool chmin(T& a,T b){bool x=a>b;if(x)a=b;return x;} #ifdef DEBUG void debugg(){cout<void debugg(const T& x,const Args&... args){cout<<" "< eratosthenes(int N){ vector ch(N+1,1); for(int i=2;i<=N;++i){ if(ch[i]==1){ for(int j=1;i*j<=N;j++){ ch[i*j]=i; } } } return ch; } ll modpow(ll a,ll n,ll m=mod){ if(n<0)return 0; ll res=1; while(n>0){ if(n&1)res=res*a%m; a=a*a%m; n>>=1; } return res; } ll modinv(ll a,ll m=mod) { ll b=m,u=1,v=0,t; while(b){ t=a/b; a-=t*b; swap(a,b); u-=t*v; swap(u,v); } return (u+m)%m; } vector> ch[100001]; int main(){ startupcpp(); auto v=eratosthenes(1000000); int N; vector a(N); ll allprod=1,alllcm=1; cin>>N; map> fmp; map smp; for(int i=0;i>a[i]; (allprod*=a[i])%=mod; int x=a[i]; while(x!=1){ if(ch[i].empty()||ch[i].back().first!=v[x]){ ch[i].emplace_back(v[x],1); }else ch[i].back().second++; x/=v[x]; } for(auto [p,cnt]:ch[i]){ auto now=fmp[p]; if(chmax(fmp[p],pair(cnt,i)))smp[p]=now.first; else chmax(smp[p],cnt); } debug(ch[i]); } vector invlist(N,1); for(auto [p,x]:fmp){ (alllcm*=modpow(p,x.first))%=mod; (invlist[x.second]*=modpow(p,x.first-smp[p]))%=mod; } debug(fmp,smp); for(int i=0;i