QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#647041 | #7434. 冷たい部屋、一人 | piggy123 | WA | 7ms | 75400kb | C++14 | 6.7kb | 2024-10-17 11:13:29 | 2024-10-17 11:13:29 |
Judging History
answer
#include <bits/stdc++.h>
#define ll long long
using namespace std;
ll a[500005],b[21][500005],c[21][500005],lg[500005],ans[500005],n,m;
struct qry {
ll l,r,id;
} p[500005],p2[500005];
vector<ll> col[500005],CC[500005];
const ll S=0;
ll S2=0;
struct block {
const ll B=600;
ll tag[805],z[500005];
void add(ll pos,ll v) {
z[pos]+=v;
tag[(pos-1)/B+1]+=v;
}
ll query(ll pos) {
ll v=(pos-1)/B+1,r=(n-1)/B+1,ans=0;
for (ll i=v+1; i<=r; i++)ans+=tag[i];
for (ll i=pos; i<=v*B; i++)ans+=z[i];
return ans;
}
} B;
ll getmx(ll l,ll r) {
ll k=lg[r-l+1];
return max(c[k][l],c[k][r-(1ll<<k)+1]);
}
ll getmi(ll l,ll r) {
ll k=lg[r-l+1];
return min(b[k][l],b[k][r-(1ll<<k)+1]);
}
ll v[1000005],cnt[1000005],d[1000005],nxt[1000005],to[1000005],T;
vector<qry> Q[1000005];
struct CHAIN{
ll vis[1000005],sz[1000005],ano[1000005];
array<ll,5> stk[1000005];
ll top,gx;
void init(ll len){
for (ll i=0;i<=len;i++)vis[i]=0,sz[i]=1,ano[i]=i;
top=0,gx=0;
}
void link(ll a){
ll l=ano[a],r=ano[a+1];
ano[l]=r,ano[r]=l;
stk[++top]={l,r,a,sz[a],sz[a+1]};
gx+=sz[a]*sz[a+1];
sz[l]=sz[r]=sz[a]+sz[a+1];
// cout<<"?" <<a<<" "<< sz[a]<<" "<<sz[a+1]<<"\n";
}
void add(ll pos){
vis[d[pos]]++;
if (vis[d[pos]]==2){
link(d[pos]);
}
}
void del(ll pos){
vis[d[pos]]--;
if (vis[d[pos]]==1){
redo();
}
}
void redo(){
auto tp=stk[top--];
gx-=tp[3]*tp[4];
ano[tp[0]]=tp[2];
ano[tp[1]]=tp[2]+1;
sz[tp[0]]=tp[3];
sz[tp[1]]=tp[4];
}
}C;
int main() {
cin.tie(0);
cout.tie(0);
ios::sync_with_stdio(false);
cin >> n >> m;
for (ll i=1; i<=n; i++) {
cin >> a[i];
col[a[i]].push_back(i);
}
for (ll i=1; i<=n; i++) {
if (i>1)lg[i]=lg[i>>1]+1;
cin >> b[0][i];
c[0][i]=b[0][i];
}
for (ll i=1; i<=20; i++) {
for (ll j=1; j+(1ll<<i)-1<=n; j++) {
b[i][j]=min(b[i-1][j],b[i-1][j+(1ll<<(i-1))]);
c[i][j]=max(c[i-1][j],c[i-1][j+(1ll<<(i-1))]);
}
}
for (ll i=1; i<=m; i++) {
cin >> p[i].l >> p[i].r;
p[i].id=i;
}
for (ll i=1; i<=n; i++) {
if (col[i].size()<=S) {
for (ll j=0; j<(ll)col[i].size(); j++) {
for (ll k=j+1; k<(ll)col[i].size(); k++) {
ll mi=getmi(col[i][j],col[i][k]),mx=getmx(col[i][j],col[i][k]);
CC[mx].push_back(mi);
}
}
}
}
sort(p+1,p+1+m,[](qry a,qry b) {
return a.r<b.r;
});
ll ps=0;
for (ll i=1; i<=m; i++) {
while (ps+1<=n&&ps+1<=p[i].r) {
for (ll j:CC[ps+1]) {
B.add(j,1);
}
ps++;
}
ans[p[i].id]+=B.query(p[i].l);
}
for (ll i=1; i<=n; i++) {
if (col[i].size()>S) {
T=i;
ll N=col[i].size();
vector<ll> book;
for (ll j=0; j+1<N; j++) {
ll mi=getmi(col[i][j],col[i][j+1]),mx=getmx(col[i][j],col[i][j+1]);
book.push_back(mi);
book.push_back(mx);
}
sort(book.begin(),book.end());
book.erase(unique(book.begin(),book.end()),book.end());
for (ll j=0; j+1<N; j++) {
ll mi=getmi(col[i][j],col[i][j+1]),mx=getmx(col[i][j],col[i][j+1]);
++cnt[lower_bound(book.begin(),book.end(),mi)-book.begin()+1];
++cnt[lower_bound(book.begin(),book.end(),mx)-book.begin()+1];
}
for (ll j=1; j<=(ll)book.size(); j++)cnt[j]+=cnt[j-1];
for (ll j=0; j+1<N; j++) {
ll mi=getmi(col[i][j],col[i][j+1]),mx=getmx(col[i][j],col[i][j+1]);
ll mip=lower_bound(book.begin(),book.end(),mi)-book.begin()+1;
ll mxp=lower_bound(book.begin(),book.end(),mx)-book.begin()+1;
d[cnt[mip]--]=j;
to[cnt[mip]+1]=mi;
d[cnt[mxp]--]=j;
to[cnt[mxp]+1]=mx;
}
ll ps=0;
for (ll j=1; j<=n; j++) {
while (ps+1<=2*(N-1)&&to[ps+1]<=j)ps++;
nxt[j]=ps;
}
S2=2*(N-1)/sqrt(m)+1;
// cout<<"! "<< i<<"\n";
for (ll j=1; j<=m; j++) {
p2[j].l=nxt[p[j].l-1]+1,p2[j].r=nxt[p[j].r],p2[j].id=p[j].id;
if ((p2[j].l-1)/S2==(p2[j].r-1)/S2) {
C.init(2*N);
for (ll k=p2[j].l; k<=p2[j].r; k++) {
C.add(k);
}
ans[p2[j].id]+=C.gx;
// cout<< p2[j].id<<":"<<C.gx<<"\n";
} else {
Q[(p2[j].l-1)/S2+1].push_back(p2[j]);
}
}
for (ll j=1; j<=(2*(N-1)-1)/S2+1; j++) {
ll lb=j*S2;
C.init(2*N);
ll L=lb+1,R=lb;
for (qry k:Q[j]) {
while (R<k.r)C.add(++R);
while (L>k.l)C.add(--L);
while (L<lb+1)C.del(L++);
ans[k.id]+=C.gx;
}
}
for (ll j=1; j<=2*(N-1); j++)cnt[j]=to[j]=nxt[j]=d[j]=0,Q[j].clear();
}
}
for (ll i=1; i<=m; i++) {
cout<< ans[i]<<"\n";
}
return 0;
}
/*
■■■■■ ■■ ■■■ ■■■ ■ ■ ■ ■■■■ ■■■■
■ ■■ ■■ ■ ■■ ■ ■■ ■ ■ ■■ ■ ■■ ■■ ■
■ ■ ■■ ■ ■ ■ ■ ■ ■ ■■■ ■■ ■■ ■ ■■
■ ■ ■■ ■ ■ ■ ■ ■ ■ ■■ ■ ■■ ■■
■ ■ ■■ ■ ■ ■■ ■■ ■■ ■
■ ■■ ■■ ■ ■■■ ■ ■■■ ■■ ■■ ■■ ■■■
■■■■■ ■■ ■ ■ ■ ■ ■■ ■■ ■■ ■■
■ ■■ ■ ■ ■ ■ ■■ ■■ ■■ ■
■ ■■ ■ ■ ■ ■ ■■ ■■ ■ ■ ■■
■ ■■ ■■ ■■ ■■ ■■ ■■ ■■ ■ ■■ ■■
■ ■■ ■■■■ ■■■■ ■■ ■■ ■■■■■■ ■■■■
*/
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 0
Wrong Answer
time: 7ms
memory: 75400kb
input:
100 100 4 1 5 1 2 1 7 5 3 7 2 3 6 6 5 3 2 2 4 1 6 5 6 2 2 2 7 6 1 3 6 3 5 6 7 6 1 2 3 3 4 2 1 1 5 4 4 3 6 7 1 1 6 1 5 6 2 3 7 4 2 4 6 7 7 3 5 3 7 2 3 3 5 1 4 7 1 2 2 5 2 2 4 3 4 7 2 7 7 3 7 3 6 6 5 4 5 4 7 6 93 52 12 70 25 36 18 37 27 99 68 40 84 3 76 57 60 19 33 41 92 87 58 13 15 43 28 63 64 59 31 ...
output:
0 0 6 35 1 0 3 1 0 0 1 1 1 16 5 6 32 11 1 0 0 1 3 12 0 0 2 0 0 0 2 0 75 0 19 1 1 20 9 15 1 0 0 7 4 0 7 2 0 1 0 1 0 3 0 11 6 12 1 1 0 0 6 1 0 0 1 3 5 0 1 2 0 1 6 1 0 0 0 1 0 0 1 0 29 1 8 2 0 12 12 0 1 0 0 8 0 1 0 14
result:
wrong answer 1st numbers differ - expected: '1', found: '0'