Submission #112596

#	Time	Username	Problem	Language	Result	Execution time	Memory
112596		BorisBarca	Ljetopica (COI19_ljetopica)	C++14	100 / 100	460 ms	32504 KiB

ljetopica

/*
	$$$$$$$\                      $$\           $$$$$$$\
	$$  __$$\                     \__|          $$  __$$\
	$$ |  $$ | $$$$$$\   $$$$$$\  $$\  $$$$$$$\ $$ |  $$ | $$$$$$\   $$$$$$\   $$$$$$$\ $$$$$$\
	$$$$$$$\ |$$  __$$\ $$  __$$\ $$ |$$  _____|$$$$$$$\ | \____$$\ $$  __$$\ $$  _____|\____$$\
	$$  __$$\ $$ /  $$ |$$ |  \__|$$ |\$$$$$$\  $$  __$$\  $$$$$$$ |$$ |  \__|$$ /      $$$$$$$ |
	$$ |  $$ |$$ |  $$ |$$ |      $$ | \____$$\ $$ |  $$ |$$  __$$ |$$ |      $$ |     $$  __$$ |
	$$$$$$$  |\$$$$$$  |$$ |      $$ |$$$$$$$  |$$$$$$$  |\$$$$$$$ |$$ |      \$$$$$$$\\$$$$$$$ |
	\_______/  \______/ \__|      \__|\_______/ \_______/  \_______|\__|       \_______|\_______|
*/
#include <bits/stdc++.h>

using namespace std;

#define PB push_back
#define MP make_pair
#define INS insert
#define LB lower_bound
#define UB upper_bound
#define pii pair <int,int>
#define pll pair <long long, long long>
#define si pair<string, int>
#define is pair<int, string>
#define X first
#define Y second
#define _ << " " <<
#define sz(x) (int)x.size()
#define all(a) (a).begin(),(a).end()
#define FOR(i, a, b) for (int i = (a); i < (b); ++i)
#define FORD(i, a, b) for (int i = (a); i > (b); --i)
#define FORR(i, l, r) for (int i = (l); i <= (r); ++i)
#define FORP(i, a, b) for ((i) = (a); (i) < (b); ++i)
#define FORA(i, x) for (auto &i : x)
#define REP(i, n) FOR(i, 0, n)
#define BITS(x) __builtin_popcount(x)
#define MSET memset
#define MCPY memcpy
#define SQ(a) (a) * (a)

typedef long long ll;
typedef long double ld;
typedef vector <int> vi;
typedef vector <pii> vpi;
typedef vector <ll> vll;
typedef vector <pll> vpl;
typedef vector <double> vd;
typedef vector <ld> vld;
typedef vector<si> vsi;
typedef vector<is> vis;
typedef vector<string> vs;
//((float) t)/CLOCKS_PER_SEC

const int MOD = 1e9 + 7;
const double PI = acos(-1);
const int LOG = 19;
const int INF = 1e9 + 10;
const ll INFL = 1e18 + 10;
const int ABC = 30;
const int dx[] = {-1, 1, 0, 0};
const int dy[] = {0, 0, -1, 1};
const int dox[] = {-1, 1, 0, 0, -1, -1, 1, 1};
const int doy[] = {0, 0, -1, 1, -1, 1, -1, 1};

inline int sum(int a, int b){
  if (a + b >= MOD)
  	return a + b - MOD;
  if (a + b < 0)
  	return a + b + MOD;
  return a + b;
}

inline void add(int &a, int b){
  a = sum(a, b);
}

inline int mul(int a, int b){
  return (ll)a * (ll)b % MOD;
}

inline int sub(int a, int b){
	return (a - b + MOD) % MOD;
}

inline int pot(ll pot, int n){
	ll ret = 1;
	if (n == 0)
		return 1;
	while (n){
		if (n & 1)
			ret = (ret * pot) % MOD;
		pot = (pot * pot) % MOD;
		n >>= 1;
	}
	return ret;
}

inline int divide(int a, int b){
	return mul(a, pot(b, MOD - 2));
}

ll lcm(ll a, ll b){
	return abs(a * b) / __gcd(a, b);
}

inline double ccw(pii A, pii B, pii C){
	return (A.X * B.Y) - (A.Y * B.X) + (B.X * C.Y) - (B.Y * C.X) + (C.X * A.Y) - (C.Y * A.X);
}

inline int CCW(pii A, pii B, pii C){
	double val = ccw(A, B, C);
	double eps = max(max(abs(A.X), abs(A.Y)), max(max(abs(B.X), abs(B.Y)), max(abs(C.X), abs(C.Y)))) / 1e9;
	if (val <= -eps)
		return -1;
	if (val >= eps)
		return 1;
	return 0;
}

void to_upper(string &x){
	REP(i, sz(x))
		x[i] = toupper(x[i]);
}

void to_lower(string &x){
	REP(i, sz(x))
		x[i] = tolower(x[i]);
}

string its(ll x){
	if (x == 0)
		return "0";
	string ret = "";
	while (x > 0){
		ret += (x % 10) + '0';
		x /= 10;
	}
	reverse(all(ret));
	return ret;
}

ll sti(string s){
	ll ret = 0;
	REP(i, sz(s)){
		ret *= 10;
		ret += (s[i] - '0');
	}
	return ret;
}

const int N = 1e3 + 10;

int n, k, memo[2][2][N][N], pom[2][2][N][N];
string a, b, m, l;
map <char, char> nxt = {{'L', 'R'}, {'R', 'L'}};
map <char, int> p = {{'L', 0}, {'R', 1}};

int cnt(int manji, int stav, int x, int y){
	if (x == n)
		return (y == 0);
	int &ret = pom[manji][stav][x][y];
	if (~ret)
		return ret;
	ret = 0;

	int move = p[m[x]];
	if (stav)
		move = !move;
	int bit = l[x] - '0';
	if (manji){
		add(ret, cnt(manji, stav, x + 1, y));
		if (y){
			stav++;
			stav %= 2;
			move = !move;
			add(ret, cnt(manji, stav, x + 1, y - 1));

			stav++;
			stav %= 2;
		}

	} else {
		if (move < bit)
			add(ret, cnt(1, stav, x + 1, y));
		else if (move == bit)
			add(ret, cnt(0, stav, x + 1, y));

		if (y){
			stav++;
			stav %= 2;

			move = !move;
			if (move < bit)
				add(ret, cnt(1, stav, x + 1, y - 1));
			else if (move == bit)
				add(ret, cnt(0, stav, x + 1, y - 1));

			stav++;
			stav %= 2;
		}

	}

	return ret;
}

int dp(int manji, int stav, int x, int y){
	if (x == n)
		return 0;
	int &ret = memo[manji][stav][x][y];
	if (~ret)
		return ret;
	ret = 0;

	int move = p[m[x]];
	if (stav)
		move = !move;
	int bit = l[x] - '0';
	if (manji){
		add(ret, mul(mul(pot(2, n - x - 1), move), cnt(manji, stav, x + 1, y))), add(ret, dp(manji, stav, x + 1, y));
		
		if (y){
			stav++;
			stav %= 2;
			move = !move;
			add(ret, mul(mul(pot(2, n - x - 1), move), cnt(manji, stav, x + 1, y - 1))), add(ret, dp(manji, stav, x + 1, y - 1));

			stav++;
			stav %= 2;
		}

	} else {
		if (move < bit)
			add(ret, mul(mul(pot(2, n - x - 1), move), cnt(1, stav, x + 1, y))), add(ret, dp(1, stav, x + 1, y));
		else if (move == bit)
			add(ret, mul(mul(pot(2, n - x - 1), move), cnt(0, stav, x + 1, y))), add(ret, dp(0, stav, x + 1, y));

		if (y){
			stav++;
			stav %= 2;

			move = !move;
			if (move < bit)
				add(ret, mul(mul(pot(2, n - x - 1), move), cnt(1, stav, x + 1, y - 1))), add(ret, dp(1, stav, x + 1, y - 1));
			else if (move == bit)
				add(ret, mul(mul(pot(2, n - x - 1), move), cnt(0, stav, x + 1, y - 1))), add(ret, dp(0, stav, x + 1, y - 1));

			stav++;
			stav %= 2;
		}

	}
	//cout << manji _ stav _ x _ y _ ret << '\n';

	return ret;
}

int main () {
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);

	cin >> n >> k;
	cin >> m;
	cin >> a >> b;
	a = a.substr(1, sz(a) - 1);
	b = b.substr(1, sz(b) - 1);
	n--;

	l = b;
	memset(pom, -1, sizeof pom);
	memset(memo, -1, sizeof memo);
	int _b = 0;
	add(_b, dp(0, 0, 0, k));
	add(_b, dp(0, 1, 0, k));
	add(_b, mul(sum(cnt(0, 0, 0, k), cnt(0, 1, 0, k)), pot(2, n)));

	l = a;
	memset(pom, -1, sizeof pom);
	memset(memo, -1, sizeof memo);
	int _a = sum(dp(0, 0, 0, k), dp(0, 1, 0, k));
	add(_a, mul(sum(cnt(0, 0, 0, k), cnt(0, 1, 0, k)), pot(2, n)));

	int sol = sub(_b, _a);
	REP(i, n)
		add(sol, mul(pot(2, n - i - 1), a[i] - '0'));
	add(sol, pot(2, n));

	cout << sol << '\n';

	return 0;
}

• Subtask 1: 8 / 8
• Subtask 2: 14 / 14
• Subtask 3: 17 / 17
• Subtask 4: 61 / 61