# | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
---|---|---|---|---|---|---|---|
919988 | Lobo | 앵무새 (IOI11_parrots) | C++17 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include "decoder.h"
#include "decoderlib.h"
#include <bits/stdc++.h>
using namespace std;
#define mp make_pair
#define pb push_back
#define fr first
#define sc second
#define int long long
// const int maxn = 1e5+10;
// const int inf = 1e18+10;
static int R1;
static int k1;
static int base1 = 1e4;
bool geq1(vector<int> s, vector<int> t) {
if(s.size() > t.size()) return 1;
if(t.size() > s.size()) return 0;
for(int i = (int) s.size() -1; i >= 0; i--) {
if(s[i] > t[i]) return 1;
if(t[i] > s[i]) return 0;
}
return 1;
}
vector<int> mult1(vector<int> s, vector<int> t) {
vector<int> ans((int) s.size()+(int) t.size(),0);
if(s == vector<int>{0} or t == vector<int>{0}) return vector<int>{0};
for(int i = 0; i < s.size(); i++) {
for(int j = 0; j < t.size(); j++) {
ans[i+j]+= s[i]*t[j];
}
}
for(int i = 0; i < ans.size(); i++) {
if(ans[i] >= base1) {
if(i+1 == (int) ans.size()) ans.pb(0);
ans[i+1]+= ans[i]/base1;
ans[i]%= base1;
}
}
while(ans.size() > 1 && ans.back() == 0) ans.pop_back();
return ans;
}
vector<int> neg(vector<int> s) {
for(int i = 0; i < s.size(); i++) {
s[i]*= -1;
}
return s;
}
vector<int> add1(vector<int> s, vector<int> t) {
vector<int> ans = s;
for(int i = 0; i < t.size(); i++) {
if(i == (int) ans.size()) ans.pb(0);
ans[i]+= t[i];
}
for(int i = 0; i < ans.size(); i++) {
if(ans[i] < 0) {
if(i+1 == (int) ans.size()) ans.pb(0);
int qt = (ans[i]-1)/base1 + 1;
ans[i]+= qt*base1;
ans[i+1]-= qt;
}
if(ans[i] >= base1) {
if(i+1 == (int) ans.size()) ans.pb(0);
ans[i+1]+= ans[i]/base1;
ans[i]%= base1;
}
}
while(ans.size() > 1 && ans.back() == 0) ans.pop_back();
return ans;
}
vector<int> trans1(int x) {
if(x == 0) return vector<int>{0};
vector<int> ans;
while(x != 0) {
ans.pb(x%base1);
x/= base1;
}
return ans;
}
int trans1back1(vector<int> s) {
int ans = 0;
int p10 = 1;
for(int i = 0; i < s.size(); i++) {
ans+= p10*s[i];
p10*= base1;
}
return ans;
}
vector<int> dp1[260][10*64+10];
void precomp1() {
// dp(i,sum) = # de (q_0,q_1,...,q_i) tal que a soma é no maximo sum
for(int i = 0; i < R1; i++) {
for(int j = 0; j <= k1; j++) {
dp1[i][j] = vector<int>{0};
}
}
for(int j = 0; j <= k1; j++) {
dp1[0][j] = trans1(j+1);
}
for(int i = 1; i < R1; i++) {
dp1[i][0] = trans1(1);
for(int j = 1; j <= k1; j++) {
for(int x = 0; x <= j; x++) {
dp1[i][j] = add1(dp1[i][j],dp1[i-1][j-x]);
}
}
}
}
vector<int> querysmaller1(vector<int> q) {
vector<int> ans = {0};
int sumq = 0;
for(int i = R1-1; i >= 0; i--) {
for(int j = 0; j <= q[i]-1 && k1-j-sumq >= 0; j++) {
if(i == 0) ans = add1(ans,trans1((int)1));
else ans = add1(ans,dp1[i-1][k1-j-sumq]);
}
sumq+= q[i];
}
if(sumq <= k1) ans = add1(ans,trans1(1));
return ans;
}
vector<int> findg1(vector<int> g) {
vector<int> ans(R1,0);
int sumans = 0;
for(int i = R1-1; i >= 0; i--) {
ans[i] = 0;
for(int j = 1; j+sumans <= k1; j++) {
ans[i] = j;
if(!(geq1(g,querysmaller1(ans)))) {
ans[i] = j-1;
break;
}
}
sumans+= ans[i];
}
return ans;
}
void decode(int32_t n, int32_t L, int32_t X[])
{
k1 = 4*n;
R1 = 200;
precomp1();
vector<int> q(R1,0);
for(int i = 0; i < L; i++) {
q[X[i]]++;
}
vector<int> vm = add1(querysmaller1(q),neg(trans1(1)));
vector<int> m(n);
for(int i = n-1; i >= 0; i--) {
vector<int> p256 = trans1(1);
for(int j = 0; j < i; j++) {
p256 = mult1(p256,trans1(256));
}
m[i] = 0;
while(geq1(vm,p256)) {
m[i]++;
vm = add1(vm,neg(p256));
}
}
for(int i =0; i < n; i++) {
output((int32_t) m[i]);
}
}