Submission #718551

#TimeUsernameProblemLanguageResultExecution timeMemory
718551vjudge1Permutation (APIO22_perm)C++17
10 / 100
1076 ms456 KiB
#include "perm.h"
#include <bits/stdc++.h>

using namespace std;

const long long MX = 1e18;

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

int randint(int l, int r){
	return uniform_int_distribution <int> (l, r) (rng);
}

double randdouble(double l, double r){
	return uniform_real_distribution <double> (l, r) (rng);
}

long long count_increasing(const vector<int>& v) {
  vector<long long> dp(v.size() + 1, 0);
  dp[0] = 1;
  for (int x : v)
  {
  	for (int i = 0; i <= x; i++)
  	{
  		dp[x+1] += dp[i];
  		dp[x+1] = min(dp[x+1],MX+1);
  	}
  }
  long long result = 0;
  for (int i = 0; i <= (int)v.size(); i++){
  	result += dp[i];
  	result = min(result,MX+1);
  }
  return result;
}
std::vector<int> construct_permutation(long long k)
{
	vector <int> ans;
	//return ans;
	for(int i = 0;; ++i){
		vector <int> tmp;
		int sz = i + 1;
		if (ans.empty()) ans.push_back(i);
		else {
			for(int p = 0; p < sz; ++p){
				vector<int> cur = ans;
				cur.insert(cur.begin() + p, i);
				if (count_increasing(cur) <= k) tmp = cur;
				else break;
			}
			ans = tmp;
		}

		if (count_increasing(ans) == k){
			break;
		}
	}

	while (1.0 * clock() / CLOCKS_PER_SEC <= 0.6){
		vector <int> ans1;
		for(int i = 0; i < (int)ans.size() - 1; ++i){
			vector <int> tmp;
			int sz = i + 1;
			if (ans1.empty()) ans1.push_back(i);
			else {
				int bound = 0;
				for(int p = 0; p < sz; ++p){
					vector<int> cur = ans1;
					cur.insert(cur.begin() + p, i);
					if (count_increasing(cur) <= k) tmp = cur, bound = p;
					else break;
				}

				int r = randint(0, 2);
				bound = max(0, bound - r);

				ans1.insert(ans1.begin() + bound, i);
				
				ans1 = tmp;
			}

			if (count_increasing(ans1) == k){
				ans = ans1;
				break;
			}
		}
	}

	return ans;
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...