Pagini recente » Cod sursa (job #2632844) | Cod sursa (job #837039) | Cod sursa (job #2236117) | Cod sursa (job #2266543) | Cod sursa (job #1134591)
#include <fstream>
#include <vector>
#include <algorithm>
using namespace std;
const int N = 16001;
struct Punct{
int x, y;
};
class Aib{
vector<int> aib[N];
int size;
inline int step(int x){
return x & -x;
}
int find(vector<int>& v, int x){
if (v.empty() || x < v[0])
return 0;
int ans = 0;
for (int step = 1 << ( 31 - __builtin_clz( v.size() ) ) ; step ; step >>= 1)
if ( (ans ^ step) < v.size() && v[ans ^ step] <= x )
ans ^= step;
return ans + 1;
}
public:
void init(int n){
size = n;
}
void insert(int x, int y){
for (int i = x ; i <= size ; i += step(i))
aib[i].push_back(y);
}
int query(int x, int y){
int ans = 0;
for (int i = x ; i ; i -= step(i))
ans += find(aib[i], y);
return ans;
}
int query(int x, int y, int z, int t){
return query(z, t) - query(z, y) - query(x, t) + query(x, y);
}
void pregen(){
for (int i = 1 ; i <= size ; i++)
sort( aib[i].begin(), aib[i].end() );
}
};
int X[N], n;
Punct P[N];
Aib A;
ifstream in("zoo.in");
ofstream out("zoo.out");
int norm(int x){
int ans = 0;
for (int step = 1 << 13 ; step ; step >>= 1)
if ( (ans ^ step) <= X[0] && X[ans ^ step] <= x )
ans ^= step;
return ans;
}
int main(){
int x, y, z, t, nrQ;
in >> n;
for (int i = 1 ; i <= n ; i++){
in >> P[i].x >> P[i].y;
X[i] = P[i].x;
}
X[0] = 1;
for (int i = 2 ; i <= n ; i++)
if (X[i] != X[i - 1])
X[ ++X[0] ] = X[i];
A.init(X[0]);
for (int i = 1 ; i <= n ; i++)
A.insert( norm(P[i].x), P[i].y );
A.pregen();
in >> nrQ;
while (nrQ--){
in >> x >> y >> z >> t;
out << A.query( norm(x - 1), y - 1, norm(z), t ) << "\n";
}
return 0;
}