# 【Nowcoder】暑期多校day5 VCD （思维 树状数组）

## 题目大意

Kanade has an infinity set H:{{(a,b) | ax,b[y1,y2]}  |  y1R,y2R,xR,y1y2} H : { { ( a , b )   |   a ≥ x , b ∈ [ y 1 , y 2 ] }     |     y 1 ∈ R , y 2 ∈ R , x ∈ R , y 1 ≤ y 2 } $H:\{ \{(a,b)\ |\ a≥ x,b∈ [y1,y2]\} \ \ |\ \ y1∈ R,y2 ∈ R,x∈ R,y1≤ y2\}$
A point set S S $S$ is good if and only if for each subset $T$$T$ of S S $S$ there exist $h\in H$$h ∈ H$ satisfy hS=T h ∩ S = T $h \cap S = T$
Now kanade has n distinct points and she want to know how many non-empty subset of these points is
good.

（真的不是我懒啊，我也转述不明白这个定义orz）

## 代码

#include <bits/stdc++.h>
using namespace std;

const int moder=998244353;
inline int add(const int &a,const int &b) {return (a+b<moder)?a+b:a+b-moder;}
inline int mul(const int &a,const int &b) {return 1ll*a*b%moder;}
inline int les(const int &a,const int &b) {return (a>=b)?a-b:a-b+moder;}

register int val=0, sign=1; char ch;
while(~(ch=getchar()) && (ch<'0' || ch>'9') && ch!='-'); ch=='-'?sign=-1:val=ch-'0';
while(~(ch=getchar()) && (ch>='0' && ch<='9')) val=(val<<1)+(val<<3)+ch-'0';
return val*sign;
}

#define x first
#define y second
const int maxn=int(1e5)+111;
typedef pair<int,int> pii;

int n;
int sta[maxn], top;
pii p[maxn];

bool cmp(const pii &a,const pii &b) {
if(a.x==b.x) return a.y<b.y;
return a.x>b.x;
}

#define lowbit(k) ((k)&(-(k)))
int sum[maxn];

void modify(int pos) {
for(;pos<=n;sum[pos]++,pos+=lowbit(pos));
return;
}
int query(int pos) {
int res=0;
for(;pos>0;res+=sum[pos],pos-=lowbit(pos));
return res;
}

void work() {
register int i;
for(i=1;i<=n;++i) {
sta[i]=p[i].y;
}

sort(p+1,p+1+n,cmp);
sort(sta+1,sta+1+n);
top=unique(sta+1,sta+1+n)-(sta+1);

for(i=1;i<=n;++i)
p[i].y=lower_bound(sta+1,sta+1+top,p[i].y)-sta;

int ans=0, pos=0;
for(i=1;i<=n+1;++i) {
if(i>1 && p[i].x!=p[pos=i-1].x) {
while(pos>=1 && p[pos].x==p[i-1].x)
modify(p[pos--].y);
}
if(i>n) break;
int cnt1=query(p[i].y-1);
int cnt2=query(top)-query(p[i].y);
}

for(i=1;i<=top;++i) {
int cnt=query(i)-query(i-1);
if(cnt==1) continue;
if(cnt&1) ans=les(ans,mul(cnt,(cnt-1)/2));
else ans=les(ans,mul(cnt/2,cnt-1));
}
printf("%d\n",ans);
return;
}

int main() {
#ifndef ONLINE_JUDGE
freopen("input.txt","r",stdin);
#endif
work();

return 0;
}
• 点赞 1
• 评论
• 分享
x

海报分享

扫一扫，分享海报

• 收藏
• 手机看

分享到微信朋友圈

x

扫一扫，手机阅读

• 打赏

打赏

yashem66

你的鼓励将是我创作的最大动力

C币 余额
2C币 4C币 6C币 10C币 20C币 50C币
• 一键三连

点赞Mark关注该博主, 随时了解TA的最新博文
10-24
08-03 180

08-04 104
08-02 150
01-22 1441
02-06 99
07-19 6520
08-01