oj.uz

Submission #798310

# Submission time Handle Problem Language Result Execution time Memory
798310 2023-07-30T15:10:21 Z khshg Parrots (IOI11_parrots) C++14
81 / 100
 2000 ms 17128 KB 
#include "encoder.h"
#include "encoderlib.h"
#include <bits/stdc++.h>
using namespace std;
 
const int BASE = 31622;
int LEN;
 
struct BIGNUM {
	vector<int> d;
}dp[41 * 64 / 10][256];
 
BIGNUM add(const BIGNUM& A, const BIGNUM& B) {
	BIGNUM res;
	int C = 0;
	for(int i = 0; i < max((int)A.d.size(), (int)B.d.size()) || C; ++i) {
		int cur = C;
		if(i < (int)A.d.size()) cur += A.d[i];
		if(i < (int)B.d.size()) cur += B.d[i];
		C = cur >= BASE;
		if(C) cur -= BASE;
		res.d.push_back(cur);
	}
	return res;
}
 
void sub(BIGNUM& A, const BIGNUM& B) {
	int C = 0;
	for(int i = 0; i < (int)B.d.size() || C; ++i) {
		A.d[i] -= C + (i < (int)B.d.size() ? B.d[i] : 0);
		C = A.d[i] < 0;
		if(C) A.d[i] += BASE;
	}
	while((int)A.d.size() && A.d.back() == 0) A.d.pop_back();
}
 
BIGNUM mul(const int& A, const BIGNUM& B) {
	BIGNUM res;
	int C = 0;
	for(int i = 0; i < (int)B.d.size() || C; ++i) {
		int cur = C;
		if(i < (int)B.d.size()) cur += A * B.d[i];
		C = cur / BASE;
		res.d.push_back(cur - C * BASE);
	}
	while((int)res.d.size() && res.d.back() == 0) res.d.pop_back();
	return res;
}
 
bool cmp(const BIGNUM& A, const BIGNUM& B) {
	if((int)A.d.size() != (int)B.d.size()) return (int)A.d.size() < (int)B.d.size();
	for(int i = (int)A.d.size() - 1; i >= 0; --i) {
		if(A.d[i] != B.d[i]) return A.d[i] < B.d[i];
	}
	return 0;
}
 
void encode(int N, int M[]) {
	LEN = 41 * N / 10;
	for(int i = 0; i < 256; ++i) dp[0][i].d = {1};// ONE
	for(int i = 1; i < LEN; ++i) dp[i][0].d = {1};// ONE
	for(int i = 1; i < LEN; ++i) {
		for(int j = 1; j < 256; ++j) {
			dp[i][j] = add(dp[i][j - 1], dp[i - 1][j]);
		}
	}
	BIGNUM val, pwi;
	pwi.d = {1};
	for(int i = N - 1; i >= 0; --i) {
		val = add(val, mul(M[i], pwi));
		if(i > 0) pwi = mul(256, pwi);
	}
	for(int i = LEN - 1; i >= 0; --i) {
		int j;
		for(j = 0; j < 256; ++j) {
			if(cmp(val, dp[i][j])) {
				break;
			}
			sub(val, dp[i][j]);
		}
		send(j);
	}
}

#include "decoder.h"
#include "decoderlib.h"
#include <bits/stdc++.h>
using namespace std;
 
const int BASE = 31622;
int LEN;
 
struct BIGNUM {
	vector<int> d;
}dp[41 * 64 / 10][256];
 
BIGNUM add(const BIGNUM& A, const BIGNUM& B) {
	BIGNUM res;
	int C = 0;
	for(int i = 0; i < max((int)A.d.size(), (int)B.d.size()) || C; ++i) {
		int cur = C;
		if(i < (int)A.d.size()) cur += A.d[i];
		if(i < (int)B.d.size()) cur += B.d[i];
		C = cur >= BASE;
		if(C) cur -= BASE;
		res.d.push_back(cur);
	}
	return res;
}
 
void sub(BIGNUM& A, const BIGNUM& B) {
	int C = 0;
	for(int i = 0; i < (int)B.d.size() || C; ++i) {
		A.d[i] -= C + (i < (int)B.d.size() ? B.d[i] : 0);
		C = A.d[i] < 0;
		if(C) A.d[i] += BASE;
	}
	while((int)A.d.size() && A.d.back() == 0) A.d.pop_back();
}
 
BIGNUM mul(const int& A, const BIGNUM& B) {
	BIGNUM res;
	int C = 0;
	for(int i = 0; i < (int)B.d.size() || C; ++i) {
		int cur = C;
		if(i < (int)B.d.size()) cur += A * B.d[i];
		C = cur / BASE;
		res.d.push_back(cur - C * BASE);
	}
	while((int)res.d.size() && res.d.back() == 0) res.d.pop_back();
	return res;
}
 
bool cmp(const BIGNUM& A, const BIGNUM& B) {
	if((int)A.d.size() != (int)B.d.size()) return (int)A.d.size() < (int)B.d.size();
	for(int i = (int)A.d.size() - 1; i >= 0; --i) {
		if(A.d[i] != B.d[i]) return A.d[i] < B.d[i];
	}
	return 0;
}
 
void decode(int N, int L, int X[]) {
	LEN = 41 * N / 10;
	sort(X, X + L);
	for(int i = 0; i < 256; ++i) dp[0][i].d = {1};// ONE
	for(int i = 1; i < LEN; ++i) dp[i][0].d = {1};// ONE
	for(int i = 1; i < LEN; ++i) {
		for(int j = 1; j < 256; ++j) {
			dp[i][j] = add(dp[i][j - 1], dp[i - 1][j]);
		}
	}
	BIGNUM val;
	for(int i = LEN - 1; i >= 0; --i) {
		for(int j = 0; j < X[i]; ++j) val = add(val, dp[i][j]);
	}
	vector<BIGNUM> p(N);
	p[0].d = {1};
	for(int i = 1; i < N; ++i) {
		p[i] = mul(256, p[i - 1]);
	}
	for(int i = N - 1; i >= 0; --i) {
		int L = 0, R = 255;
		while(L < R) {
			int M = (L + R + 1) / 2;
			if(cmp(val, mul(M, p[i]))) {
				R = M - 1;
			} else {
				L = M;
			}
		}
		output(L);
		sub(val, mul(L, p[i]));
	}
}

Subtask #1
17.0 / 17.0

# Verdict Execution time Memory Grader output
1 Correct 20 ms 4536 KB Output is correct

Subtask #2
17.0 / 17.0

# Verdict Execution time Memory Grader output
1 Correct 86 ms 4752 KB Output is correct
2 Correct 180 ms 5228 KB Output is correct
3 Correct 222 ms 5816 KB Output is correct
4 Correct 247 ms 5864 KB Output is correct

Subtask #3
18.0 / 18.0

# Verdict Execution time Memory Grader output
1 Correct 83 ms 4800 KB Output is correct
2 Correct 162 ms 5252 KB Output is correct
3 Correct 217 ms 5784 KB Output is correct
4 Correct 254 ms 5932 KB Output is correct

Subtask #4
29.0 / 29.0

# Verdict Execution time Memory Grader output
1 Correct 86 ms 4788 KB Output is correct
2 Correct 246 ms 5868 KB Output is correct
3 Correct 336 ms 6468 KB Output is correct
4 Correct 610 ms 8668 KB Output is correct
5 Correct 661 ms 8920 KB Output is correct
6 Correct 679 ms 9140 KB Output is correct

Subtask #5
0.0 / 19.0

# Verdict Execution time Memory Grader output
1 Correct 244 ms 5808 KB Output is correct - P = 4.062500
2 Correct 788 ms 9128 KB Output is correct - P = 4.093750
3 Correct 741 ms 9468 KB Output is correct - P = 4.090909
4 Correct 1531 ms 13564 KB Output is correct - P = 4.100000
5 Correct 1930 ms 16028 KB Output is correct - P = 4.100000
6 Execution timed out 2001 ms 17128 KB Time limit exceeded
7 Execution timed out 2024 ms 17092 KB Time limit exceeded