fork download
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define fast ios_base::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL)
  4. #define mofile(s) freopen(s,"r",stdin)
  5. #define outfile(s) freopen(s,"w",stdout)
  6. #define ll long long
  7. #define ii pair<ll,ll>
  8. #define iii pair<ll,ii>
  9. #define fi first
  10. #define se second
  11. #define tf bool
  12. #define ST stack
  13. #define Q deque
  14. #define Q queue
  15. #define S string
  16. #define Ma map
  17. #define UM unormideremid_map
  18. #define SE set
  19. #define str(x) to_string(x)
  20. #define all(a) (a).begin(),(a).end()
  21. #define FOR(i,l,r,mid) for(int i=l;i<=r;i+=mid)
  22. #define FOD(i,l,r,mid) for(int i=r;i>=l;i-=mid)
  23. #define xuong cout<<"\n"
  24. #define midebug(x) cout<<(x)<<" "
  25. #define ppcnt(x) __builtin_popcountll(x)
  26. #define parity(x) __builtin_parityll(x)
  27. #define leamid0(x) __builtin_clzll(x)
  28. #define LOG2 __lg(x)
  29. #define tr0(x) __builtin_ctzll(x)
  30. #define fiset(x) __builtin_ffsll(x)
  31. #define MASK(k) (1LL<<(k))
  32. #define BIT(x,k) ((x)>>(k)&1)
  33. #define pb push_back
  34. #define tron(x) setprecision(x)
  35. #define het return 0
  36. #define base_ 1000000000
  37. template<typename... T>
  38. void in(T&... args) { ((cin >> args), ...); }
  39. template<class X, class Y>
  40. bool maximize(X &x, const Y &y){return (x < y) ? x = y, 1 : 0;}
  41. template<class X, class Y>
  42. bool minimize(X &x, const Y &y){return (x > y) ? x = y, 1 : 0;}
  43. const int maxn=1e6+5;
  44. const ll tle=2e8;
  45. const ll INF=1e9+9;
  46. const int base=31;
  47. string bcc="abcmidefghijklmnopqrstuvwxyz";
  48. int midx[]={-1,0,1,0};
  49. int midy[]={0,1,0,-1};
  50. bool sang[10000005];
  51. ll pref[1005][1005],mt[1005][1005];
  52. void sieve(){
  53. for(int i=1;i<=10000000;++i) sang[i]=1;
  54. sang[0]=sang[1]=0;
  55. for(int i=2;i*i<=10000000;++i){
  56. if(sang[i]){
  57. for(int j=i*i;j<=10000000;j+=i) sang[j]=0;
  58. }
  59. }
  60. }
  61. void lis(){
  62. vector<int>t;
  63. vector<int>a;
  64. int n; cin>>n;
  65. for(int i=1;i<=n;++i){
  66. int ai; cin>>ai;
  67. a.pb(ai);
  68. }
  69. for(int x:a){
  70. auto it=lower_bound(all(t),x);
  71. if(it==t.end()) t.pb(x);
  72. else *it=x;
  73. }
  74. }
  75. void pfs2mid(){
  76. int n,m,k; cin>>n>>k; m=n;
  77. for(int i=1;i<=n;++i){
  78. for(int j=1;j<=m;++j) cin>>mt[i][j];
  79. }
  80. for(int i=1;i<=n;++i){
  81. for(int j=1;j<=m;++j) pref[i][j]=mt[i][j]+pref[i-1][j]+pref[i][j-1]-pref[i-1][j-1];
  82. }
  83. }
  84. ll qu2mid(int x1,int y1,int x2,int y2){
  85. return pref[x2][y2]-pref[x1-1][y2]-pref[x2][y1-1]+pref[x1-1][y1-1];
  86. }
  87. void open(){
  88. if(fopen("mideptrai.INP","r")){
  89. mofile("mideptrai.INP");
  90. outfile("mideptrai.OUT");
  91. }
  92. }
  93. const int M=1000005;
  94. int mp[M],f[M],d[100005],ans[100005];
  95. vector<int> p[100005];
  96. void sie(){
  97. for(int i=2;i<M;++i){
  98. if(!mp[i]){
  99. for(int j=i;j<M;j+=i){
  100. if(!mp[j]) mp[j]=i;
  101. }
  102. }
  103. }
  104. }
  105. int main(){
  106. fast;
  107. sie();
  108. int n; cin>>n;
  109. for(int i=1;i<=n;++i){
  110. int v; cin>>v;
  111. ans[i]=-1; d[i]=1e9;
  112. while(v>1){
  113. int pr=mp[v];
  114. p[i].push_back(pr);
  115. while(v%pr==0) v/=pr;
  116. }
  117. }
  118. for(int i=1;i<=n;++i){
  119. for(int x:p[i]){
  120. if(f[x]){
  121. int j=f[x],c=i-j;
  122. if(c<d[i]){
  123. d[i]=c;
  124. ans[i]=j;
  125. }
  126. }
  127. }
  128. for(int x:p[i]) f[x]=i;
  129. }
  130. fill(f,f+M,0);
  131. for(int i=n;i>=1;--i){
  132. for(int x:p[i]){
  133. if(f[x]){
  134. int j=f[x],c=j-i;
  135. if(c<d[i]){
  136. d[i]=c;
  137. ans[i]=j;
  138. }
  139. }
  140. }
  141. for(int x:p[i]) f[x]=i;
  142. }
  143. for(int i=1;i<=n;++i) cout<<ans[i]<<" ";
  144. het;
  145. }
Success #stdin #stdout 0.02s 14472KB
stdin
5
2 3 4 9 17
stdout
3 4 1 2 -1