結果

問題 No.493 とても長い数列と文字列(Long Long Sequence and a String）
ユーザー e869120e869120
提出日時 2017-03-08 20:40:27
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 800 ms
コード長 2,276 bytes
コンパイル時間 842 ms
コンパイル使用メモリ 72,056 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-23 19:32:50
合計ジャッジ時間 3,430 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 115
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<algorithm>
#include<string>
using namespace std;
long long K, L, R, U, V; string Y;
long long sums(long long P, long long Q) {
	if (P >= (1LL << Q))return (1LL << 62);
	if (P < 0)return -(1LL << 62);
	long long ret = 0, U = P;
	for (int i = 1; i <= Q; i++) {
		ret += ((U + 1) / 2) * to_string(i*i).size();
		U -= (U + 1) / 2;
	}
	return ret;
}
long long modpow(long long a, long long b, long long m) {
	long long p = 1, q = a;
	for (int i = 0; i < 63; i++) {
		if ((b / (1LL << i)) % 2 == 1)p *= q;
		q *= q;
		p %= m; q %= m;
	}
	return p;
}
void search() {
	//左の境界を探索
	long long L1 = -1, R1 = (1LL << K) + 1, M1;
	while (true) {
		M1 = (L1 + R1) / 2;
		long long Y1 = sums(M1 - 1, K), Y2 = sums(M1, K);
		if (Y1 < L && Y2 >= L) {
			U = M1; U++;
			if (Y2 - L >= 1) {
				long long I = M1, J = 1;
				while (I % 2 == 0) { I /= 2; J++; }
				Y += to_string(J*J).substr(to_string(J*J).size() - (Y2 - L), (Y2 - L));
			}
			break;
		}
		if (Y2 < L)L1 = M1;
		if (Y1 >= L)R1 = M1;
	}
	//右の境界を探索
	L1 = -1; R1 = (1LL << K) + 1;
	while (true) {
		M1 = (L1 + R1) / 2;
		long long Y1 = sums(M1, K), Y2 = sums(M1 + 1, K);
		if (Y1 <= R && Y2 > R) {
			V = M1;
			if (R - Y1 >= 1) {
				long long I = M1, J = 1; I++;
				while (I % 2 == 0) { I /= 2; J++; }
				Y += to_string(J*J).substr(0, (R - Y1));
			}
			break;
		}
		if (Y2 <= R)L1 = M1;
		if (Y1 > R)R1 = M1;
	}
}
int main() {
	cin >> K >> L >> R; L--; if (K >= 62)K = 62;
	if (R > sums((1LL << K) - 1, K)) { cout << "-1" << endl; return 0; }

	search();

	//端数部分の和と積を求める
	long long wa = 0, seki = 1;
	for (int i = 0; i < Y.size(); i++) {
		long long G = (Y[i] - '0'); if (G == 0)G = 10;
		wa += G; seki *= G;
	}
	//区間の部分の積と和を求める
	U--;
	for (int i = 1; i <= K; i++) {
		long long U1 = (U + 1) / 2;
		long long U2 = (V + 1) / 2;

		string T = to_string(i*i);
		long long w1 = 0, w2 = 1;
		for (int j = 0; j < T.size(); j++) {
			long long G = (T[j] - '0'); if (G == 0)G = 10;
			w1 += G; w2 *= G;
		}
		//cout << U2 - U1 << ' ' << U << ' ' << V << endl;
		wa += w1*(U2 - U1);
		seki *= modpow(w2, U2 - U1, 1000000007);
		seki %= 1000000007;
		U -= U1; V -= U2;
	}
	//出力
	cout << wa << ' ' << seki << endl;
	return 0;
}
0