标签：Sort index 10 int 最小值 lld% rk Strings Revision
-
题意：
一个s串满足s[i] = i % 10,给出p,d数组的构造方法，每次将s串p[i]位换成d[i], 总共得到n + 1个串，对它们进行排名。第i次操作后的排名为r[i];
输出 \((\sum_{i=0}^{n}(r_{i}*10000019^{i}))\%(1e9+7)\) -
题解：
-
很关键的一点是p[]是一个排列，所以对s串每一位只修改一次。我们从低位到高位去修改，如果修改p[i]位要变成的数小于该位原来的数(p[i]%10>d[i])，那么修改p[i]之后位得到的的串字典序一定小于修改p[i]前得到的字符串。所以可以用这一次修改将字符串排名分成两类，一类占字典序大的部分，一类占字典序小的部分。再继续分治对每个部分用其P[i]最小的位进行划分。一直划分下去就能得到所有的排名。
-
关键是当划分成两部分后怎么快速确定这个处理区间的最小值位置，RMQ问题。这里可以用O(n)的笛卡尔树来处理。笛卡尔树根节点是整个区间最小值的位置。左子节点是最小位置左边区间的最小值的位置，右子节点是最小值右边区间最小值所在位置。
-
分治时，如果 \(p[i]\%10>d[i]\) , 其中i是当前区间最小值所在位置，那么其左边的区间操作形成的字符串rk+右边个数(占rk高位)。右边不影响(占rk低位)
dfs(ls[index], l, index, rank + (p[index] % 10 > d[index]) * (r - index)); -
如果 \(p[i]\%10<d[i]\) 那i右边区间操作形成的字符串rk+左边个数(占rk高位)，左边不受影响(占rk低位)
dfs(rs[index], index + 1, r, rank + (p[index] % 10 < d[index]) * (index - l + 1)); -
如果 \(p[i]\%10=d[i]\) 时，这一位操作不会影响任何其他位的rk, 一个技巧，让p[i]=INF, 这样这一位就在迪卡尔树的叶子节点上，最后处理时直接确定rk就行。同时如果处理到边界，及叶子节点，说明这个操作处理这个区间的所能处理的最高位，能影响它的低位都处理了，根据前面的rk累计可以直接确定其rk。
if(p[index] == INF || l >= r) { for(int i = l; i <= r; ++ i) rk[i] = rank + (i - l); return; }
-
-
代码：
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N = 2e6+10; const ll mod = 1e9 + 7; const ll INF = 1 << 30; // int T, n; ll p[N], d[N]; int rk[N]; int ls[N], rs[N]; int st[N], top; void Build() { st[0]=0; for(int i=0;i<n;i++){ int k=st[0]; while(k>0&&p[st[k]]>p[i]) k--; if(k) rs[st[k]]=i; if(k<st[0]) ls[i]=st[k+1]; st[++k]=i; st[0]=k; } } void dfs(int index, int l, int r, int rank){ if(p[index] == INF || l >= r) { for(int i = l; i <= r; ++ i) rk[i] = rank + (i - l); return; } dfs(ls[index], l, index, rank + (p[index] % 10 > d[index]) * (r - index)); dfs(rs[index], index + 1, r, rank + (p[index] % 10 < d[index]) * (index - l + 1)); } int main() { scanf("%d",&T); while(T --){ scanf("%d",&n); ll pseed, pa, pb, pmod; ll dseed, da, db, dmod; scanf("%lld%lld%lld%lld",&pseed, &pa, &pb, &pmod); scanf("%lld%lld%lld%lld",&dseed, &da, &db, &dmod); for(int i = 0; i < n; ++ i) p[i] = i; for(int i = 1; i < n; ++ i){ swap(p[pseed % (i + 1)], p[i]); pseed = (1ll * pseed * pa % pmod + 1ll * pb) % pmod; } for(int i = 0; i < n; ++ i){ d[i] = dseed % 10; dseed = (1ll * dseed * da % dmod + 1ll * db) % dmod; if(p[i] % 10 == d[i]) p[i] = INF; } Build(); for(int i = 0; i <= n; ++ i) rk[i] = 0; dfs(st[1], 0, n, 0); ll res = 0, temp = 1; for(int i = 0; i <= n; ++ i){ res = (res + 1ll * rk[i] * temp % mod) % mod; temp = temp * 10000019ll % mod; } printf("%lld\n",res); } return 0; }
标签：Sort,index,10,int,最小值,lld%,rk,Strings,Revision 来源： https://www.cnblogs.com/A-sc/p/13371368.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。