結果

問題 No.1653 Squarefree
ユーザー aut_jul_dur
提出日時 2021-08-20 23:00:04
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,209 bytes
コンパイル時間 3,508 ms
コンパイル使用メモリ 231,688 KB
実行使用メモリ 20,104 KB
最終ジャッジ日時 2024-10-14 06:11:01
合計ジャッジ時間 23,311 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27 WA * 11
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define all(x) (x).begin(),(x).end()
const long long MOD = 1000000007;
const long long INF = 9*1e18;
using ll = long long;
using mint = modint1000000007;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }
long double PI = 3.14159265358979;
vector<ll> primelist;
bool Prime(ll num)
{
    if (num < 2) return false;
    else if (num == 2) return true;
    else if (num % 2 == 0) return false;
 
    double sqrtNum = sqrt(num);
    for (ll i = 3; i <= sqrtNum; i += 2)
    {
        if (num % i == 0)
        {
            return false;
        }
    }
    return true;
}
void primeinput(){
  for(ll i = 2;i<1002000;i++){
    if(Prime(i)){
      primelist.push_back(i);
    }
  }
}
int main(){
  ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  ll L,R;
  cin>>L>>R;
  primeinput();
  vector<ll> A,B;
  for(ll i = L;i<=R;i++){
    A.push_back(i);
    B.push_back(0);
  }
  rep(i,primelist.size()){
    ll temp = 0;
    if(primelist[i] >= L){
      temp = primelist[i]-L;
    }
    else{
      temp = primelist[i]-(L%primelist[i]);
    }
   // cout << primelist[i] << " " << temp << endl;
    while(1){
      if(temp >= A.size()){
        break;
      }
      ll cnt = 0;
      while(A[temp]%primelist[i] == 0){
        if(primelist[i] == 7){
          //cout << A[temp] << " " << primelist[i] << endl;
        }
        A[temp]/=primelist[i];
        cnt++;
      }
      if(cnt >= 2){
        B[temp] = 1;
      }
      temp += primelist[i];
    }
  }
  /*rep(i,B.size()){
    cout << B[i] << " ";
  }
  cout << endl;*/
  ll ans = 0;
  rep(i,A.size()){
    if(B[i] == 1){
      continue;
    }
    if(A[i] == 1 && B[i] != 1){
      ans++;
      continue;
    }
    ll temp = sqrt(A[i]);
    //cout << temp << endl;
    bool d = 0;
    for(int k = -1;k<=1;k++){
      if((temp+k)*(temp+k) == A[i]){
        d = 1;
      }
    }
    if(!d){
      //cout << L+i << endl;
      ans++;
    }
  }
  cout << ans << endl;
  return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0