【做题记录】 [JLOI2011]不等式组

P5482 [JLOI2011]不等式组

超烦人的细节题！(本人调了两天 QAQ )

这里介绍一种只用到一只树状数组的写法(离线)。

树状数组的下标是：所有可能出现的数据进行离散化之后的值。

其含义为：当 \(x\) 离散化后值为 \(i\) 时能满足的不等式个数为 \(query(i)\) 个。

处理数据

首先我们先读入所有数据，并对数据处理：

\(\text{Add} ~a_i~b_i~c_i\) ：

若 \(a_i>0\) 将 \(a_ix+b_i>c_i\) 转化成 \(x\ge t_i\) 的形式 。

若 \(a_i<0\) 将 \(a_ix+b_i>c_i\) 转化成 \(x\le t_i\) 的形式 。

并将 \(t_i\) 丢进离散化的序列中。

注意：所有的除法运算都是向 \(0\) 取整，还要注意除法变号问题等等。

\(\text{Del}\) ：

在处理 \(\text{Add}\) 时提前记录第 \(x\) 个 \(\text{Add}\) 操作所对应的输入操作编号。

\(\text{Query}\) ：

将 \(k_i\) 丢进离散化序列中。

之后将序列中的数离散化，给 \(\text{Add}\) 中的 \(t_i\) 和 \(\text{Query}\) 中的\(k_i\) 都附上一个离散化后的值( \(Instead_i\) ) 。

计算答案

\(\text{Add}\) ：

若 \(a_i>0\) 则在 \([t_i,+\infty)\) 区间内的 \(Instead_x\) 都可以使不等式成立。

同理，若 \(a_i<0\) 则在 \((-\infty,t_i]\) 区间内的 \(Instead_x\) 都可以使不等式成立。

在区间内加 \(1\) 即可 。

\(\text{Del}\) 和 \(\text{Add}\) 几乎一致，变为区间减 \(1\) 。

\(\text{Query}~k_i\) 即可直接查询并输出 \(Query(Instead_i)\) 。

最后附上 100pts 代码：

#include<bits/stdc++.h>
using namespace std;
#define Maxn 100005
#define inf 0x7f7f7f7f
typedef long long ll;
inline int rd()
{
	 int x=0;
     char ch,t=0;
     while(!isdigit(ch = getchar())) t|=ch=='-';
     while(isdigit(ch)) x=x*10+(ch^48),ch=getchar();
     return x=t?-x:x;
}
int n,tmp,tot,cnt;
map<int,int> mp;
int Ins_val[Maxn],hist[Maxn];
struct Data
{
	 int opt,t,Ins;
	 int pre,used; 
}a[Maxn];
int tree[Maxn];
inline int lowbit(int x){ return x&(-x); }
void add(int x,int k)
{
	 while(x<=tot+1) tree[x]+=k,x+=lowbit(x);
}
int query(int x)
{
	 int ret=0;
	 while(x) ret+=tree[x],x-=lowbit(x);
	 return ret;
}
int main()
{
     //freopen(".in","r",stdin);
     //freopen(".out","w",stdout);
	 n=rd();
	 string opt;
	 for(int i=1,x,y,z,A,B,C;i<=n;i++)
	 {
	 	 cin>>opt;
	 	 if(opt=="Add")
		 {
		 	 A=rd(),B=rd(),C=rd(),hist[++cnt]=i;
			 a[i].opt=2-(A>=0); // 当 a>=0 时 opt=1 ，否则 opt=2 
			 if(A==0) a[i].t=(B>C)?(-inf+1):inf;
			 else if(A>0) a[i].t=(C-B)/A+(((C-B)>=0)?1:(((C-B)/A*A==(C-B))?1:0));
			 else a[i].t=(C-B)/A-(((C-B)>=0)?1:(((C-B)/A*A==(C-B))?1:0));
			 Ins_val[++tmp]=a[i].t;
		 }
	 	 if(opt=="Del")
		 {
		 	 a[i].pre=hist[rd()];
			 a[i].opt=a[a[i].pre].opt+2; // 当 a>=0 时 opt=3 ，否则 opt=4 
			 a[i].t=a[a[i].pre].t;
		 }
	 	 if(opt=="Query") a[i].opt=5,a[i].t=rd(),Ins_val[++tmp]=a[i].t;
	 }
	 sort(Ins_val+1,Ins_val+1+tmp);
	 Ins_val[0]=-inf;
	 for(int i=1;i<=tmp;i++) if(Ins_val[i]!=Ins_val[i-1]) mp[Ins_val[i]]=++tot;
	 for(int i=1;i<=n;i++) a[i].Ins=mp[a[i].t];
	 for(int i=1;i<=n;i++)
	 {
	 	 if(a[i].t==inf) continue; // a==0 && b<c
	 	 if(a[i].opt==1) add(tot+1,-1),add(a[i].Ins,1); //  a>0 || (a==0 && b>c)
	 	 if(a[i].opt==2) add(a[i].Ins+1,-1),add(1,1); //  a<0 
	 	 if(a[i].opt==3 && !a[a[i].pre].used) add(tot+1,1),add(a[i].Ins,-1),a[a[i].pre].used=1; // a>0 || (a==0 && b>c)
	 	 if(a[i].opt==4 && !a[a[i].pre].used) add(a[i].Ins+1,1),add(1,-1),a[a[i].pre].used=1; // a<0
	 	 if(a[i].opt==5) printf("%d\n",query(a[i].Ins));
	 }
     //fclose(stdin);
     //fclose(stdout);
     return 0;
}