結果

問題 No.1233 割り切れない気持ち
ユーザー ゆにぽけ
提出日時 2024-01-29 00:28:40
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 414 ms / 3,153 ms
コード長 1,247 bytes
コンパイル時間 1,219 ms
コンパイル使用メモリ 125,152 KB
最終ジャッジ日時 2025-02-19 00:37:36
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
#include <array>
#include <iterator>
#include <string>
#include <cctype>
#include <cstring>
#include <cstdlib>
#include <cassert>
#include <cmath>
#include <ctime>
#include <iomanip>
#include <numeric>
#include <stack>
#include <queue>
#include <map>
#include <unordered_map>
#include <set>
#include <unordered_set>
#include <bitset>
#include <random>
#include <utility>
#include <functional>
using namespace std;
void Main()
{
	int N;
	cin >> N;
	const int M = (int)2e5;
	vector<int> cnt(M + 1);
	for(int i = 0;i < N;i++)
	{
		int a;
		cin >> a;
		cnt[a]++;
	}
	vector<long long> S(M + 1);
	for(int i = 1;i <= M;i++)
	{
		S[i] = S[i - 1];
		S[i] += (long long) cnt[i] * i;
	}
	long long ans = 0;
	for(int i = 1;i <= M;i++) if(cnt[i])
	{
		//a = b * (a / b) + a % b
		//a % b = a - b * (a / b)
		long long cur = 0;
		cur += (long long) N * i;
		for(int j = 1;j * j <= i;j++)
		{
			cur -= (long long) j * (i / j) * cnt[j];
			if(j != i / j)
			{
				long long sum = S[i / j] - S[i / (j + 1)];
				cur -= sum * j;
			}
		}
		ans += cur * cnt[i];
	}
	cout << ans << endl;
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int tt = 1;
	/* cin >> tt; */
	while(tt--) Main();
}
0