#614909 (C++17) No.1383 Numbers of Product

提出ソース

結果

問題	No.1383 Numbers of Product
ユーザー	👑 tute7627
提出日時	2021-02-07 21:47:27
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	567 ms / 2,000 ms
コード長	4,728 bytes
コンパイル時間	2,411 ms
コンパイル使用メモリ	211,508 KB
実行使用メモリ	68,352 KB
最終ジャッジ日時	2024-11-22 23:18:08
合計ジャッジ時間	17,834 ms
ジャッジサーバーID （参考情報）	judge5 / judge2

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,816 KB
testcase_01	AC	2 ms 6,816 KB
testcase_02	AC	2 ms 6,816 KB
testcase_03	AC	2 ms 6,820 KB
testcase_04	AC	2 ms 6,816 KB
testcase_05	AC	2 ms 6,816 KB
testcase_06	AC	2 ms 6,820 KB
testcase_07	AC	342 ms 43,264 KB
testcase_08	AC	331 ms 42,624 KB
testcase_09	AC	407 ms 52,096 KB
testcase_10	AC	546 ms 67,584 KB
testcase_11	AC	535 ms 66,688 KB
testcase_12	AC	545 ms 67,456 KB
testcase_13	AC	519 ms 64,768 KB
testcase_14	AC	445 ms 55,424 KB
testcase_15	AC	356 ms 45,568 KB
testcase_16	AC	545 ms 66,816 KB
testcase_17	AC	2 ms 6,816 KB
testcase_18	AC	2 ms 6,816 KB
testcase_19	AC	2 ms 6,816 KB
testcase_20	AC	2 ms 6,816 KB
testcase_21	AC	2 ms 6,820 KB
testcase_22	AC	2 ms 6,820 KB
testcase_23	AC	2 ms 6,820 KB
testcase_24	AC	2 ms 6,816 KB
testcase_25	AC	2 ms 6,820 KB
testcase_26	AC	2 ms 6,816 KB
testcase_27	AC	2 ms 6,816 KB
testcase_28	AC	10 ms 6,816 KB
testcase_29	AC	422 ms 53,248 KB
testcase_30	AC	556 ms 68,224 KB
testcase_31	AC	2 ms 6,820 KB
testcase_32	AC	133 ms 16,640 KB
testcase_33	AC	2 ms 6,820 KB
testcase_34	AC	2 ms 6,820 KB
testcase_35	AC	503 ms 62,464 KB
testcase_36	AC	485 ms 60,032 KB
testcase_37	AC	553 ms 68,352 KB
testcase_38	AC	555 ms 68,224 KB
testcase_39	AC	559 ms 68,224 KB
testcase_40	AC	558 ms 68,224 KB
testcase_41	AC	567 ms 68,352 KB
testcase_42	AC	549 ms 68,352 KB
testcase_43	AC	556 ms 68,224 KB
testcase_44	AC	554 ms 68,224 KB
testcase_45	AC	548 ms 68,224 KB
testcase_46	AC	2 ms 6,816 KB
testcase_47	AC	548 ms 66,816 KB
testcase_48	AC	546 ms 67,328 KB
testcase_49	AC	550 ms 67,968 KB
testcase_50	AC	545 ms 67,968 KB
testcase_51	AC	2 ms 6,816 KB
testcase_52	AC	2 ms 6,816 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

//#define _GLIBCXX_DEBUG

#include<bits/stdc++.h>
using namespace std;

#define endl '\n'
#define lfs cout<<fixed<<setprecision(10)
#define ALL(a)  (a).begin(),(a).end()
#define ALLR(a)  (a).rbegin(),(a).rend()
#define spa << " " <<
#define fi first
#define se second
#define MP make_pair
#define MT make_tuple
#define PB push_back
#define EB emplace_back
#define rep(i,n,m) for(ll i = (n); i < (ll)(m); i++)
#define rrep(i,n,m) for(ll i = (ll)(m) - 1; i >= (ll)(n); i--)
using ll = long long;
using ld = long double;
const ll MOD1 = 1e9+7;
const ll MOD9 = 998244353;
const ll INF = 1e18;
using P = pair<ll, ll>;
template<typename T1, typename T2>bool chmin(T1 &a,T2 b){if(a>b){a=b;return true;}else return false;}
template<typename T1, typename T2>bool chmax(T1 &a,T2 b){if(a<b){a=b;return true;}else return false;}
ll median(ll a,ll b, ll c){return a+b+c-max({a,b,c})-min({a,b,c});}
void ans1(bool x){if(x) cout<<"Yes"<<endl;else cout<<"No"<<endl;}
void ans2(bool x){if(x) cout<<"YES"<<endl;else cout<<"NO"<<endl;}
void ans3(bool x){if(x) cout<<"Yay!"<<endl;else cout<<":("<<endl;}
template<typename T1,typename T2>void ans(bool x,T1 y,T2 z){if(x)cout<<y<<endl;else cout<<z<<endl;}  
template<typename T>void debug(T &v,ll h,ll w,string sv=" "){for(ll i=0;i<h;i++){cout<<v[i][0];for(ll j=1;j<w;j++)cout<<sv<<v[i][j];cout<<endl;}};
template<typename T>void debug(T &v,ll n,string sv=" "){if(n!=0)cout<<v[0];for(ll i=1;i<n;i++)cout<<sv<<v[i];cout<<endl;};
template<typename T>vector<vector<T>>vec(ll x, ll y, T w){vector<vector<T>>v(x,vector<T>(y,w));return v;}
ll gcd(ll x,ll y){ll r;while(y!=0&&(r=x%y)!=0){x=y;y=r;}return y==0?x:y;}
vector<ll>dx={1,-1,0,0,1,1,-1,-1};vector<ll>dy={0,0,1,-1,1,-1,1,-1};
template<typename T>vector<T> make_v(size_t a,T b){return vector<T>(a,b);}
template<typename... Ts>auto make_v(size_t a,Ts... ts){return vector<decltype(make_v(ts...))>(a,make_v(ts...));}
template<typename T1, typename T2>ostream &operator<<(ostream &os, const pair<T1, T2>&p){return os << p.first << " " << p.second;}
template<typename T>ostream &operator<<(ostream &os, const vector<T> &v){for(auto &z:v)os << z << " ";cout<<"|"; return os;}
template<typename T>void rearrange(vector<int>&ord, vector<T>&v){
  auto tmp = v;
  for(int i=0;i<tmp.size();i++)v[i] = tmp[ord[i]];
}
template<typename Head, typename... Tail>void rearrange(vector<ll>&ord,Head&& head, Tail&&... tail){
  rearrange(ord, head);
  rearrange(ord, tail...);
}
template<typename T> vector<int> ascend(const vector<T>&v){
  vector<int>ord(v.size());iota(ord.begin(),ord.end(),0);
  sort(ord.begin(),ord.end(),[&](int i,int j){return v[i]<v[j];});
  return ord;
}
template<typename T> vector<int> descend(const vector<T>&v){
  vector<int>ord(v.size());iota(ord.begin(),ord.end(),0);
  sort(ord.begin(),ord.end(),[&](int i,int j){return v[i]>v[j];});
  return ord;
}
ll FLOOR(ll n,ll div){return n>=0?n/div:(n-div+1)/div;}
ll CEIL(ll n,ll div){return n>=0?(n+div-1)/div:n/div;}
//mt19937 mt(chrono::steady_clock::now().time_since_epoch().count());
int popcount(ll x){return __builtin_popcountll(x);};
int poplow(ll x){return __builtin_ctzll(x);};
int pophigh(ll x){return 63 - __builtin_clzll(x);};

template< typename T = int >
struct edge {
  int to;
  T cost;
  int id;
  edge() = default;
  edge(int to, T cost = 1, int id = -1):to(to), cost(cost), id(id){}
  operator int() const { return to; }
};

template<typename T>
using Graph = vector<vector<edge<T>>>;
template<typename T>
Graph<T> readGraph(int n,int m,int indexed=1,bool directed=false,bool weighted=false){
  Graph<T> ret(n);
  for(int es = 0; es < m; es++){
    int u,v,w=1;
    cin>>u>>v;u-=indexed,v-=indexed;
    if(weighted)cin>>w;
    ret[u].emplace_back(v,w,es);
    if(!directed)ret[v].emplace_back(u,w,es);
  }
  return ret;
}
template<typename T>
Graph<T> readParent(int n,int indexed=1,bool directed=true){
  Graph<T>ret(n);
  for(int i=1;i<n;i++){
    int p;cin>>p;
    p-=indexed;
    ret[p].emplace_back(i);
    if(!directed)ret[i].emplace_back(p);
  }
  return ret;
}

int main(){
  cin.tie(nullptr);
  ios_base::sync_with_stdio(false);
  ll res=0,buf=0;
  bool judge = true;
  ll n,k,m;cin>>n>>k>>m;
  map<ll,ll>mp;
  using I=__int128_t;
  for(ll i=1;i<=n;i++){
    if(I(i)*(i+k)*(i+k+k)>n)break;
    I now=I(i)*(i+k);
    I mul=i+k;
    while(now<=n){
      mul+=k;
      now*=mul;
      if(now<=n)mp[now]++;
    }
  }
  auto f=[&](ll lim){
    ll ok=0,ng=1e12;
    while(ng-ok>=2){
      I mid=(ok+ng)/2;
      if(mid*(mid+k)>lim)ng=mid;
      else ok=mid;
    }
    return ok;
  };
  ll rem=f(n);
  for(auto z:mp){
    ll x=f(z.fi);
    if(I(x)*(x+k)==z.fi){
      z.se++;
      rem--;
    }
    if(z.se==m)res++;
  }
  if(m==1)res+=rem;
  cout<<res<<endl;
  return 0;
}

yukicoder

結果

テストケース

ソースコード