結果

問題 No.1260 たくさんの多項式
ユーザー SSRS
提出日時 2020-10-16 22:03:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 187 ms / 2,000 ms
コード長 1,125 bytes
コンパイル時間 1,379 ms
コンパイル使用メモリ 166,920 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-20 21:47:33
合計ジャッジ時間 7,561 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 61
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
const long long MOD = 1000000007;
long long modpow(long long a, long long b){
	long long ans = 1;
	while (b > 0){
		if (b % 2 == 1){
			ans *= a;
			ans %= MOD;
		}
		a *= a;
		a %= MOD;
		b /= 2;
	}
	return ans;
}
long long modinv(long long a){
	return modpow(a, MOD - 2);
}
long long ap(long long a, long long b){
  return (a + b) % MOD * ((b - a + 1) % MOD) % MOD * modinv(2) % MOD;
}
int main(){
  long long N;
  cin >> N;
  long long ans = 0;
  for (long long i = 2; i <= N; i++){
    long long N2 = N;
    int cnt = 0;
    long long sum = 0;
    while (N2 > 0){
      cnt++;
      sum += N2 % i;
      N2 /= i;
    }
    if (cnt >= 3){
      ans += sum;
      ans %= MOD;
    } else {
      break;
    }
  }
  for (long long i = 1; i <= N; i++){
    long long mn = N / (i + 1) + 1;
    long long mx = N / i;
    if (i >= mn){
      break;
    }
    long long sum = N % MOD * ((mx - mn + 1) % MOD) % MOD - i * ap(mn, mx) % MOD;
    if (sum < 0){
      sum += MOD;
    }
    ans += sum;
    ans += i * ((mx - mn + 1) % MOD) % MOD;
    ans %= MOD;
  }
  cout << ans << endl;
}
0