結果

問題 No.757 チャンパーノウン定数 (2)
ユーザー FF256grhyFF256grhy
提出日時 2018-12-05 07:15:02
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 686 ms / 2,000 ms
コード長 3,096 bytes
コンパイル時間 1,824 ms
コンパイル使用メモリ 172,280 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-18 05:51:21
合計ジャッジ時間 12,267 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 51
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

typedef long long   signed int LL;
typedef long long unsigned int LU;

#define incID(i, l, r) for(int i = (l)    ; i <  (r); i++)
#define incII(i, l, r) for(int i = (l)    ; i <= (r); i++)
#define decID(i, l, r) for(int i = (r) - 1; i >= (l); i--)
#define decII(i, l, r) for(int i = (r)    ; i >= (l); i--)
#define  inc(i, n) incID(i, 0, n)
#define inc1(i, n) incII(i, 1, n)
#define  dec(i, n) decID(i, 0, n)
#define dec1(i, n) decII(i, 1, n)

#define inII(v, l, r) ((l) <= (v) && (v) <= (r))
#define inID(v, l, r) ((l) <= (v) && (v) <  (r))

#define PB push_back
#define EB emplace_back
#define MP make_pair
#define FI first
#define SE second
#define PQ priority_queue

#define  ALL(v)  v.begin(),  v.end()
#define RALL(v) v.rbegin(), v.rend()
#define  FOR(it, v) for(auto it =  v.begin(); it !=  v.end(); ++it)
#define RFOR(it, v) for(auto it = v.rbegin(); it != v.rend(); ++it)

template<typename T> bool   setmin(T & a, T b) { if(b <  a) { a = b; return true; } else { return false; } }
template<typename T> bool   setmax(T & a, T b) { if(b >  a) { a = b; return true; } else { return false; } }
template<typename T> bool setmineq(T & a, T b) { if(b <= a) { a = b; return true; } else { return false; } }
template<typename T> bool setmaxeq(T & a, T b) { if(b >= a) { a = b; return true; } else { return false; } }
template<typename T> T gcd(T a, T b) { return (b == 0 ? a : gcd(b, a % b)); }
template<typename T> T lcm(T a, T b) { return a / gcd(a, b) * b; }

// ---- ----

int L, b;

void f(string & s, int v, int p) {
	dec1(i, p) {
		s[i] += v % b;
		v /= b;
		if(s[i] >= '0' + b) { s[i - 1]++; s[i] -= b; }
		else if(v == 0) { return; }
	}
	assert(false);
}

string diff(const string & s, const string & t) { // s - t
	assert(s.size() == t.size());
	
	string ans = string(s.size(), '0');
	dec(i, s.size()) {
		ans[i] += s[i] - t[i];
		if(ans[i] < '0') { ans[i - 1]--; ans[i] += b; }
	}
	return ans;
}

int strmod(const string & s, int x) {
	int v = 0;
	inc(i, s.size()) { v = (b * v + (s[i] - '0')) % x; }
	return v;
}

int si(const string & s) {
	int v = 0;
	inc(i, s.size()) { v = b * v + (s[i] - '0'); }
	return v;
}

string is(int v) {
	string s;
	while(v) { s += '0' + v % b; v /= b; }
	reverse(ALL(s));
	return s;
}

string strdiv(const string & s, int x) {
	string t = s;
	
	string ans = string(t.size(), '0');
	int p = 0;
	inc(i, t.size()) {
		int y = si(t.substr(p, i - p + 1));
		ans[i] += y / x;
		y %= x;
		string u = is(y);
		u = string(i - p + 1 - u.size(), '0') + u;
		incII(j, p, i) { t[j] = u[j - p]; }
		while(p < i && t[p] == '0') { p++; }
	}
	return ans;
}

char g(string & n) {
	L = n.size() + 20;
	
	string m = "1";
	n = string(L - n.size(), '0') + n;
	m = string(L - m.size(), '0') + m;
	
	int l = -1;
	inc1(i, 100001) {
		string m2 = m;
		f(m2, i * (b - 1), L - i);
		if(n < m2) { l = i; break; }
		m = m2;
	}
	
	string d = diff(n, m);
	
	int    mo = strmod(d, l);
	string di = strdiv(d, l);
	di[L - l] += 1;
	
	return di[L - l + mo];
}

int main() {
	string n;
	cin >> b >> n;
	
	cout << g(n) << endl;
	
	return 0;
}
0