結果
問題 | No.1567 Integer Coefficient Equation |
ユーザー |
![]() |
提出日時 | 2021-06-26 14:31:18 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 212 ms / 2,000 ms |
コード長 | 1,514 bytes |
コンパイル時間 | 2,337 ms |
コンパイル使用メモリ | 195,116 KB |
最終ジャッジ日時 | 2025-01-22 13:23:56 |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 38 |
ソースコード
#define rep(i, n) for (int i = 0; i < (int)(n); i++) #define ALL(v) v.begin(), v.end() typedef long long ll; #include <bits/stdc++.h> using namespace std; template<typename T> vector<T> smallest_prime_factors(T n){ vector<T> spf(n+1); for(int i=0;i<=n;i++) spf[i]=i; for(T i=2;i*i<=n;i++){ if(spf[i]==i){ //素数だったら for(T j=i*i;j<=n;j+=i){ if(spf[j]==j) spf[j]=i; //iを持つ整数かつまだ素数が決まっていないなら } } } return spf; } //高速因数分解(クエリ<=10^5,整数<=10^6程度) template<typename T> bool factolization(T x,vector<T> &spf){ int cnt=0; while(x!=1) { x/=spf[x]; cnt++; if(cnt>=3){ break; } } if(cnt>=3) return true; return false; } int main(){ ios::sync_with_stdio(false); std::cin.tie(nullptr); constexpr int MAX=400040; auto spf=smallest_prime_factors(MAX); vector<int> A(MAX+1); for(int i=1;i<MAX;i++){ if(factolization(i,spf)==1) A[i]=1; } //素数の3乗 for(int i=1;i*i*i<MAX;i++){ if(spf[i]==i) A[i*i*i]=0; } vector<int> S(MAX+2); for(int i=1;i<MAX;i++) S[i]=A[i]+S[i-1]; int t; cin>>t; while(t--){ ll p,l,r; cin>>p>>l>>r; l-=p,r-=p; if(l>0 && r>0) cout<<S[r]-S[l-1]<<endl; else if(l<0 && r<0) cout<<S[-l]-S[-r-1]<<endl; else if(l==0 && r==0) cout<<1<<endl; else if(l==0) cout<<1+S[r]<<endl; else if(r==0) cout<<1+S[-l]<<endl; else cout<<1+S[r]+S[-l]<<endl; } return 0; }