#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <set>
#include <map>
#include <utility>
#include <queue>
#include <vector>
#include <string>
#include <cstring>
#define REP(a,n) for (int a = 0; a<(n); ++a)
#define FOR(a,b,c) for (int a = (b); a<=(c); ++a)
#define FORD(a,b,c) for (int a = (b); a>=(c); --a)
#define FOREACH(a,v) for (auto a : v)
#define MP make_pair
#define PB push_back
template<class T> inline int size(const T&t) { return t.size(); }
using namespace std;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<string> vs;
typedef long long LL;
///////////////////////////////
#define DEBUG 0
int XS, YS, Q;
string s1, s2;
int podnies[12][3001][3001];
int ans[100000];
int main() {
cin >> XS >> YS >> Q >> s1 >> s2;
FOR(i, 0, 11) {
for (int sr = 1 << i; sr <= XS; sr += (1 << (i + 1))) {
FORD(x, sr - 1, sr - (1 << i)) {
int podn_dod_nad = YS;
//podnies[i][x][YS] = YS;
FORD(y, YS - 1, 0) {
if (s1[x] == s2[y]) {
podnies[i][x][y] = podn_dod_nad;
podn_dod_nad = y;
}
else
podnies[i][x][y] = max(y, min(podnies[i][x + 1][y], podn_dod_nad));
podn_dod_nad = max(podn_dod_nad, podnies[i][x + 1][y]);
}
}
REP(y, YS + 1)
podnies[i][sr][y] = y;
FOR(x, sr + 1, min(sr + (1 << i) - 1, XS)) {
int podn_dod_pod = 0;
//podnies[i][x][0] = 0;
FOR(y, 1, YS) {
if (s1[x - 1] == s2[y - 1]) {
podnies[i][x][y] = podn_dod_pod;
podn_dod_pod = y;
}
else
podnies[i][x][y] = min(y, max(podnies[i][x - 1][y], podn_dod_pod));
podn_dod_pod = min(podn_dod_pod, podnies[i][x - 1][y]);
}
}
}
#if DEBUG
if ((1 << i) <= XS) {
cerr << "i = " << i << endl;
FORD(y, YS, 0) {
REP(x, XS + 1)
cerr << podnies[i][x][y] << " ";
cerr << endl;
}
}
#endif
}
REP(q, Q) {
int x0, x1, y0, y1;
cin >> x0 >> x1 >> y0 >> y1;
--x0;
--y0;
int i = 11;
while ((x0 & (1 << i)) == (x1 & (1 << i))) --i;
#if DEBUG
cerr << "zap: " << x0 << " " << y0 << " " << x1 << " " << y1 << "; i=" << i << endl << "gorne:" << endl;
#endif
//int sr = x1 & ~((1 << i) - 1);
int gorne[YS + 1], dolne[YS + 1], wyn[YS + 1];
memset(gorne, 0, (YS + 1) * sizeof(int));
memset(dolne, 0, (YS + 1) * sizeof(int));
FORD(y, YS - 1, y0)
++gorne[podnies[i][x0][y]];
int ll = 0, best = 0;
FORD(y, YS, y0) {
ll += gorne[y];
wyn[y] = ll;
#if DEBUG
cerr << ll << endl;
#endif
--ll;
}
FOR(y, 1, y1)
++dolne[podnies[i][x1][y]];
ll = 0;
FOR(y, 0, y1) {
ll += dolne[y];
if (y >= y0)
best = max(best, ll + wyn[y]);
--ll;
}
cout << best << endl;
}
}
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 | #include <iostream> #include <cstdio> #include <cstdlib> #include <algorithm> #include <set> #include <map> #include <utility> #include <queue> #include <vector> #include <string> #include <cstring> #define REP(a,n) for (int a = 0; a<(n); ++a) #define FOR(a,b,c) for (int a = (b); a<=(c); ++a) #define FORD(a,b,c) for (int a = (b); a>=(c); --a) #define FOREACH(a,v) for (auto a : v) #define MP make_pair #define PB push_back template<class T> inline int size(const T&t) { return t.size(); } using namespace std; typedef pair<int, int> pii; typedef vector<int> vi; typedef vector<string> vs; typedef long long LL; /////////////////////////////// #define DEBUG 0 int XS, YS, Q; string s1, s2; int podnies[12][3001][3001]; int ans[100000]; int main() { cin >> XS >> YS >> Q >> s1 >> s2; FOR(i, 0, 11) { for (int sr = 1 << i; sr <= XS; sr += (1 << (i + 1))) { FORD(x, sr - 1, sr - (1 << i)) { int podn_dod_nad = YS; //podnies[i][x][YS] = YS; FORD(y, YS - 1, 0) { if (s1[x] == s2[y]) { podnies[i][x][y] = podn_dod_nad; podn_dod_nad = y; } else podnies[i][x][y] = max(y, min(podnies[i][x + 1][y], podn_dod_nad)); podn_dod_nad = max(podn_dod_nad, podnies[i][x + 1][y]); } } REP(y, YS + 1) podnies[i][sr][y] = y; FOR(x, sr + 1, min(sr + (1 << i) - 1, XS)) { int podn_dod_pod = 0; //podnies[i][x][0] = 0; FOR(y, 1, YS) { if (s1[x - 1] == s2[y - 1]) { podnies[i][x][y] = podn_dod_pod; podn_dod_pod = y; } else podnies[i][x][y] = min(y, max(podnies[i][x - 1][y], podn_dod_pod)); podn_dod_pod = min(podn_dod_pod, podnies[i][x - 1][y]); } } } #if DEBUG if ((1 << i) <= XS) { cerr << "i = " << i << endl; FORD(y, YS, 0) { REP(x, XS + 1) cerr << podnies[i][x][y] << " "; cerr << endl; } } #endif } REP(q, Q) { int x0, x1, y0, y1; cin >> x0 >> x1 >> y0 >> y1; --x0; --y0; int i = 11; while ((x0 & (1 << i)) == (x1 & (1 << i))) --i; #if DEBUG cerr << "zap: " << x0 << " " << y0 << " " << x1 << " " << y1 << "; i=" << i << endl << "gorne:" << endl; #endif //int sr = x1 & ~((1 << i) - 1); int gorne[YS + 1], dolne[YS + 1], wyn[YS + 1]; memset(gorne, 0, (YS + 1) * sizeof(int)); memset(dolne, 0, (YS + 1) * sizeof(int)); FORD(y, YS - 1, y0) ++gorne[podnies[i][x0][y]]; int ll = 0, best = 0; FORD(y, YS, y0) { ll += gorne[y]; wyn[y] = ll; #if DEBUG cerr << ll << endl; #endif --ll; } FOR(y, 1, y1) ++dolne[podnies[i][x1][y]]; ll = 0; FOR(y, 0, y1) { ll += dolne[y]; if (y >= y0) best = max(best, ll + wyn[y]); --ll; } cout << best << endl; } } |
English