結果

問題 No.1233 割り切れない気持ち
ユーザー hamamuhamamu
提出日時 2020-09-18 23:29:05
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 5,916 bytes
コンパイル時間 2,679 ms
コンパイル使用メモリ 208,032 KB
実行使用メモリ 12,212 KB
最終ジャッジ日時 2023-09-04 19:44:39
合計ジャッジ時間 11,864 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
8,768 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 3 ms
4,380 KB
testcase_03 AC 14 ms
4,384 KB
testcase_04 AC 14 ms
4,376 KB
testcase_05 AC 3 ms
4,380 KB
testcase_06 AC 3 ms
4,380 KB
testcase_07 AC 1,313 ms
4,548 KB
testcase_08 AC 2,556 ms
5,624 KB
testcase_09 TLE -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
using     ll=long long;
using    vll=vector<   ll>;
using   vvll=vector<  vll>;
using  vvvll=vector< vvll>;
using vvvvll=vector<vvvll>;
constexpr ll INF = 1LL << 60;
struct Fast{ Fast(){ cin.tie(0); ios::sync_with_stdio(false); cout<<fixed<<setprecision(numeric_limits<double>::max_digits10); } } fast;
#define REPS(i, S, E) for (ll i = (S); i <= (E); i++)
#define REP(i, N) REPS(i, 0, (N)-1)
#define DEPS(i, S, E) for (ll i = (E); i >= (S); i--)
#define DEP(i, N) DEPS(i, 0, (N)-1)
#define rep(i, S, E)  for (ll i = (S); i <= (E); i++)
#define dep(i, E, S)  for (ll i = (E); i >= (S); i--)
#define each(e, v) for (auto&& e : v)
#define ALL(v) (v).begin(), (v).end()
#define RALL(v) (v).rbegin(), (v).rend()
template<class T> inline bool chmax(T &a, T b) { if (a < b) { a = b; return true; }return false; }
template<class T> inline bool chmin(T &a, T b) { if (a > b) { a = b; return true; }return false; }
template<class T> inline T MaxE(vector<T>&v,ll S,ll E){ T m=v[S]; rep(i,S,E)chmax(m,v[i]); return m; }
template<class T> inline T MinE(vector<T>&v,ll S,ll E){ T m=v[S]; rep(i,S,E)chmin(m,v[i]); return m; }
template<class T> inline T MaxE(vector<T> &v) { return MaxE(v,0,(ll)v.size()-1); }
template<class T> inline T MinE(vector<T> &v) { return MinE(v,0,(ll)v.size()-1); }
template<class T> inline T Sum(vector<T> &v,ll S,ll E){ T s=T(); rep(i,S,E)s+=v[i]; return s; }
template<class T> inline T Sum(vector<T> &v) { return Sum(v,0,v.size()-1); }
template<class T> inline ll sz(T &v){ return (ll)v.size(); }
inline ll CEIL(ll a,ll b){ return (a<0) ? -(-a/b) : (a+b-1)/b; }
inline ll FLOOR(ll a,ll b){ return -CEIL(-a,b); }

//vector用テンプレート
template<class T> inline vector<T>& operator+=(vector<T> &a,const vector<T> &b){ for (ll i=0; i<(ll)a.size(); i++) a[i]+=b[i];  return a; }
template<class T> inline vector<T>& operator-=(vector<T> &a,const vector<T> &b){ for (ll i=0; i<(ll)a.size(); i++) a[i]-=b[i];  return a; }
template<class T> inline vector<T>& operator*=(vector<T> &a,const vector<T> &b){ for (ll i=0; i<(ll)a.size(); i++) a[i]*=b[i];  return a; }
template<class T> inline vector<T>& operator/=(vector<T> &a,const vector<T> &b){ for (ll i=0; i<(ll)a.size(); i++) a[i]/=b[i];  return a; }
template<class T> inline vector<T>& operator%=(vector<T> &a,const vector<T> &b){ for (ll i=0; i<(ll)a.size(); i++) a[i]%=b[i];  return a; }
template<class T,class S> inline vector<T>& operator+=(vector<T> &a,S b){ for (T &e: a) e+=b;  return a; }
template<class T,class S> inline vector<T>& operator-=(vector<T> &a,S b){ for (T &e: a) e-=b;  return a; }
template<class T,class S> inline vector<T>& operator*=(vector<T> &a,S b){ for (T &e: a) e*=b;  return a; }
template<class T,class S> inline vector<T>& operator/=(vector<T> &a,S b){ for (T &e: a) e/=b;  return a; }
template<class T,class S> inline vector<T>& operator%=(vector<T> &a,S b){ for (T &e: a) e%=b;  return a; }
template<class T,class S> inline vector<T> operator+(const vector<T> &a,S b){ vector<T> c=a; return c+=b; }
template<class T,class S> inline vector<T> operator-(const vector<T> &a,S b){ vector<T> c=a; return c-=b; }
template<class T,class S> inline vector<T> operator*(const vector<T> &a,S b){ vector<T> c=a; return c*=b; }
template<class T,class S> inline vector<T> operator/(const vector<T> &a,S b){ vector<T> c=a; return c/=b; }
template<class T,class S> inline vector<T> operator%(const vector<T> &a,S b){ vector<T> c=a; return c%=b; }
template<class T,class S> inline vector<T> operator-(S b,const vector<T> &a){ vector<T> c=-a; return c+=b; }
template<class T> inline vector<T> operator-(const vector<T> &a){ vector<T> c=a; return c*=(-1); }
template<class T> inline ostream &operator<<(ostream &os,const vector<T> &a){ for (ll i=0; i<(ll)a.size(); i++) os<<(i>0?" ":"")<<a[i];  return os; }

template<class T,class S> inline ll CeilIdx(vector<T> &v,S x){
	return lower_bound(ALL(v),x)-v.begin();
}
/*
ll idx=CeilIdx(v,x);
xがv[N-1]より大きい時はNが返る。
*/

template<class T,class S> inline ll FloorIdx(vector<T> &v,S x){
	return upper_bound(ALL(v),x)-v.begin()-1;
}
/*
ll idx=FloorIdx(v,x);
xがv[0]未満の時は-1が返る。
*/

template <class T> struct CumulativeSumG{
	ll s=0,e=-1;  vector<T> cum;
	CumulativeSumG(){};
	CumulativeSumG(vector<T> &v): e(v.size()-1),cum(e+2) { Ini(v); }
	CumulativeSumG(vector<T> &v,ll S,ll E): s(S),e(E),cum(e-s+2){ Ini(v); }
	void init(vector<T> &v){ e=v.size()-1; cum.resize(e+2); Ini(v); }
	void init(vector<T> &v,ll S,ll E){ s=S; e=E; cum.resize(e-s+2); Ini(v); }
	void add(T x) { e++; cum.push_back(cum.back() + x); }
	T operator () (ll begin,ll end){ return cum[end+1-s] - cum[begin-s]; }
	void Ini(vector<T> &v) { rep(i,s,e) cum[i+1-s] = cum[i-s] + v[i]; }
};
using CumulativeSum = CumulativeSumG<ll>;


#if 0
#include <atcoder/all>
using namespace atcoder;
#endif


void solve()
{
	ll M=500;
	ll n;  cin >> n;
	vector<ll> a(n);  rep(i,0,n-1){ ll a_; cin>>a_; a[i]=a_; }
	sort(ALL(a));
	vll cna;
	cna.push_back(1);
	rep(i,1,n-1){
		if (a[i-1]==a[i]){
			cna.back()++;
		}
		else{
			a[sz(cna)]=a[i];
			cna.push_back(1);
		}
	}
	a.resize(cna.size());

	vll b=a;
	b*=cna;
	CumulativeSum cm(b);
	CumulativeSum cntcm(cna);
	n=sz(a);

	ll ans=0;
	rep(i,0,n-1){
		ll aa=a[i];
		ll cnaa=cna[i];
		rep(p,0,M){
			ll lw,up;
			lw=(aa/(p+1))+1;
			if (p==0)up=INF;
			else up=aa/p;

			ll lwi=CeilIdx(a,lw);
			ll upi=FloorIdx(a,up);

			//ans+=(upi-lwi+1)*aa-p*cm(lwi,upi);
			ans+= (cntcm(lwi,upi)*aa-p*cm(lwi,upi)) * cnaa;
		}

		ll up=aa/(M+1);
		ll upi=FloorIdx(a,up);
		rep(k,0,upi){
			ll r=aa%a[k];
			ans+=r*cnaa*cna[k];
		}
	}
	cout << ans << '\n';
}


int main(){
#if 1
	solve();
#else
	ll t;  cin >> t;
	rep(i, 0, t-1){
		solve();
	}
#endif
	return 0;
}
0