GYM 102439 Equal Mod Segments

题意是，给你一个序列，问有多少区间 [l,r] 满足 a[l]%a[l+1]%a[l+2].....%a[r]==a[r]%a[r-1]%a[r-2]...%a[l]?

我们考虑一个区间 [l,r] 假设区间内最小值在pos处取得，那么区间的合法性可以表示为

a[l]%a[l+1]%a[l+2].....%a[pos]==a[r]%a[r-1]%a[r-2]...%a[pos].因为在mod a[pos]之后，值是小于apos的，继续取余比较大的数也没有什么意义。

如果pos==l || pos==r 呢？如果区间内没有等于apos的数，那么区间一定是不合法的。只有区间内有多个等于apos的数才有可能合法，这里需要特判一下。

所以只要统计一个每个最小值左右相等的数的乘积，注意不要计算重复。

可以分治解决这个问题。我们每次找到区间内最小值（多个取最左边的），递归左区间，递归右区间。

对于每一个分治的区间L,R，我们用set去维护所有数向左取模和所有数向右取模。

那么计算的时候，先分别把左区间向右的mod apos，右区间向左的mod pos，这里注意，取模的时候先二分找到第一个大于apos的值，再向后依次取模，这样保证每次操作都会使一个数变小。（这里的复杂度是nlognlogn的，因为一个最多有效取模log次，那么n个数最多取模nlog次，再算上set的复杂度）。

然后统计答案，合并左右区间。

为了保证分治的复杂度，就是要启发式分治嘛。

统计答案的时候要遍历小的区间，二分大的区间。合并的时候也是把小区间插入到大的区间里去。

这样复杂度也是nloglog的嘛。

#include <bits/stdc++.h>using namespace std;#define N 200025
#define M 610
#define ll long long
#define mod 1000000007
#define go(i,a,b) for(int i=(a);i<=(b);i++)
#define dep(i,a,b) for(int i=(a);i>=(b);i--)
#define pb push_back
#define inf 0x3f3f3f3f
#define ld long double
#define pii pair<int,int>
#define pdd pair<double,double>
#define vi vector<int>
#define madd(a,b) (a+=(b)%mod)%=mod
#define lowb(c,len,x) lower_bound(c+1,c+len+1,x)-c
#define uppb(c,len,x) upper_bound(c+1,c+len+1,x)-c
#define ls i*2+1
#define rs i*2+2
//#define mid (l+r)/2
#define lson l,mid,ls
#define rson mid+1,r,rs
//#define root 1,cnt,0
#define ms(a,b) memset(a,b,sizeof a)
#define muti int T,cas=1;cin>>T;while(T--)
#define lll __int128
#define si short int
#define fi first
#define se second
#define l(x) (x&-x)
#define G(x) for( int ii=h[x];ii;ii=eg[ii].n )int mm[2*N],n,las[N];
int dp[2*N][20],a[N],cnt;
set<pii>st[N];
set<pii>::iterator it,it1,itt;
ll ans=0;
void init()
{mm[0]=-1;go(i,1,n){dp[i][0]=i;mm[i]=((i&(i-1))==0)?mm[i-1]+1:mm[i-1];}go(j,1,mm[n])for(int i=1; i+(1<<j)-1<=n; i++)dp[i][j]=a[dp[i][j-1]]<=a[dp[i+(1<<(j-1))][j-1]]?dp[i][j-1]:dp[i+(1<<(j-1))][j-1];
}
int quy(int x,int y){if(x>y)return n+1;int k=mm[y-x+1];return a[dp[x][k]]<=a[dp[y-(1<<k)+1][k]]?dp[x][k]:dp[y-(1<<k)+1][k];
}
int join(int x,int y,int z){ //合并两个set，插入a[mid]if(st[x].size()>st[y].size())swap(x,y);//if(!st[x].size())return y;for(it=st[x].begin();it!=st[x].end();it++){it1=st[y].lower_bound(pii(it->fi,0));if(it1==st[y].end()||it1->fi!=it->fi)st[y].insert(*it);else {pii tmp=pii(it1->fi,it1->se+it->se);st[y].erase(it1);st[y].insert(tmp);}}st[x].clear();it=st[y].lower_bound(pii(z,0));if(it!=st[y].end()&&it->fi==z){pii tmp=pii(it->fi,it->se+1);st[y].erase(it);st[y].insert(tmp);}else st[y].insert(pii(z,1));return y;
}
void cal(int x,int y){ //统计答案if(st[x].size()>st[y].size())swap(x,y);if(!st[x].size())return ;for(it=st[x].begin();it!=st[x].end();it++){it1=st[y].lower_bound(pii(it->fi,0));if(it1!=st[y].end()&&it1->fi==it->fi)ans+=1ll*(it->se)*(it1->se);}
}
//void pri(int x){
//    cout<<":::::::::::::"<<endl;
//    for(it=st[x].begin();it!=st[x].end();it++){
//        cout<<it->fi<<'-'<<it->se<<' ';
//    }
//    cout<<endl<<endl;
//}
void upd(int x,int y){ //对set所有的值 mod a[mid]for(it=st[x].lower_bound(pii(y,0));it!=st[x].end();){pii tmp=pii(it->fi%y,it->se);it1=it;it++;st[x].erase(it1);itt=st[x].lower_bound(pii(tmp.fi,0));if(itt!=st[x].end()&&itt->fi==tmp.fi){pii tmp1=pii(itt->fi,itt->se+tmp.se);st[x].erase(itt);st[x].insert(tmp1);}else st[x].insert(tmp);}
}pii solve(int l,int r,int fl=0){ if(fl)ans++;if(l>=r){if(l>r)return pii(++cnt,++cnt);//cout<<fl<<endl;st[++cnt].insert(pii(a[l],1));st[++cnt].insert(pii(a[l],1));return pii(cnt-1,cnt);}int mid=quy(l,r);//递归左右区间pii pl=solve(l,mid-1),pr=solve(mid+1,r,a[quy(mid+1,r)]==a[mid]?1:0),res;//分别取模a[mid]upd(pr.fi,a[mid]);upd(pl.se,a[mid]);//特判区间内有多个相同的最小值it=st[pr.fi].lower_bound(pii(0,0));if(it!=st[pr.fi].end()&&it->fi==0)ans+=fl*(it->se);//统计答案cal(pl.se,pr.fi);//合并左右区间，a[mid]res=pii(join(pl.fi,pr.fi,a[mid]),join(pl.se,pr.se,a[quy(mid+1,r)]==a[mid]?0:a[mid]));return res;
}int main()
{cin>>n;go(i,1,n)scanf("%d",&a[i]);a[n+1]=inf;init();solve(1,n);cout<<ans+n<<endl;
}
/*
10
10 9 8 7 1 2 3 4 5 6
*/

GYM 102439 Equal Mod Segments相关推荐

Diffie–Hellman key exchange
原文地址:https://en.wikipedia.org/wiki/Diffie%E2%80%93Hellman_key_exchange Diffie–Hellman key exchange F ...
MOD - Power Modulo Inverted(SPOJ3105) + Clever Y(POJ3243) + Hard Equation (Gym 101853G ) + EXBSGS
思路: 前两题题面相同,代码也相同,就只贴一题的题面了.这三题的意思都是求A^X==B(mod P),P可以不是素数,EXBSGS板子题. SPOJ3105题目链接:https://www.spoj. ...
Codeforces gym 2013-2014 Samara SAU ACM ICPC Quarterfinal Qualification Contest
题目地址:http://codeforces.com/gym/100247 2014.10.27版 26号训练的题,今天发现A题有更高效的解法,用到了前缀最大值和后缀最大值 A. The Power ...
Codeforces gym 100685 C. Cinderella 水题
C. Cinderella Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/gym/100685/problem/ ...
codeforces Gym 100187B B. A Lot of Joy
B. A Lot of Joy Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/gym/100187/proble ...
TVM/Relay 的 PartitionGraph()(mod) 函数讨论整理
TVM/Relay 的 PartitionGraph()(mod) 函数讨论整理 TVM/Relay 的图形分区功能.以下简单示例,错误信息. PartitionGraph() 函数指定图形是用带有 ...
codeforces 610D D. Vika and Segments(离散化+线段树+扫描线算法)
题目链接: D. Vika and Segments time limit per test 2 seconds memory limit per test 256 megabytes input s ...
Codeforces Gym 100513G G. FacePalm Accounting 暴力
G. FacePalm Accounting Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/gym/100513 ...
Codeforces Round #337 (Div. 2) D. Vika and Segments 线段树扫描线
D. Vika and Segments 题目连接: http://www.codeforces.com/contest/610/problem/D Description Vika has an i ...

GYM 102439 Equal Mod Segments

GYM 102439 Equal Mod Segments相关推荐

最新文章

热门文章