2018.10.29 正睿停课训练 Day11
一场rating排名从11掉到40+ ==。掉就掉吧
状态很迷 全程写T1的。。随机算法。。（一开始就想错了，误以为它正确性很高）
T2有想法但没调出来
T3有50暴力但是直接没看==
A 线段树什么的最讨厌了(思路 DFS)
容易发现线段树上区间([l,r])的父节点只有(4)种情况：([l,r+len-1],[l,r+len],[l-len-1,r],[l-len,r])。
其实也比较好想，因为只需对父节点区间长度的奇偶性分类讨论一下。
那么我们可以直接从给定区间([l,r])不断枚举父节点，直到找到一个根节点([0,n])？
显然只有(log)层，那这样是(4^{log n}=n^2)的？注意到有限制为((frac{l}{l-r+1})^2leq2000)，所以(n^2)是可过的。（好吧我也不太明白为啥是(frac{l}{l-r+1})？）
注意要剪枝，不去搜明显不属于线段树上的区间（不然还是暴力分）。
//371ms 504kb
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
int Ans;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
void DFS(int l,int r)
{
if(l<0||r>=Ans) return;
if(!l) {Ans=r; return;}
int len=r-l+1;
if(l-len>=2*len||!(l-len))DFS(l-len,r);
if(l-len>=2*(len+1)||!(l-len-1)) DFS(l-len-1,r);
if(l>2*len) DFS(l,r+len);
if(l>2*(len-1)) DFS(l,r+len-1);
}
int main()
{
for(int T=read(); T--; )
{
int l=read(),r=read(),lim=read();
if(l==r) {printf("%d
",r); continue;}
Ans=lim+1, DFS(l,r), printf("%d
",Ans>lim?-1:Ans);
}
return 0;
}
B 已经没有什么好害怕的了(差分 前缀和)
将每个配对的括号表示成([l,r))的形式。
首先可以将区间([l,r))加一。
然后对于左端点(l')在(r)位置的匹配括号对，显然可以将([l,r))的贡献都加到([l',r'))上；
对于右端点(r')在(l)位置的匹配括号对，显然也可以将([l,r))的贡献都加到([l',r'))上。
然后就可以神奇差分了。。然后没看懂。。学套路吧。。
//899ms 21408kb
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
#define mod 1000000007
#define Mod(x) x>=mod&&(x-=mod)
#define Add(x,v) (x+=v)>=mod&&(x-=mod)
typedef long long LL;
const int N=1e6+5;
int sk[N],L[N],R[N],sum[N],delta[N];
LL sum2[N];
char s[N];
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
int main()
{
for(int T=read(),n; T--; )
{
memset(L,0,sizeof L);
memset(R,0,sizeof R);
memset(sum,0,sizeof sum);
memset(delta,0,sizeof delta);
scanf("%s",s+1), n=strlen(s+1);
int top=0;
for(int i=1; i<=n; ++i)
if(s[i]=='(') sk[++top]=i;
else if(top) L[i+1]=sk[top], R[sk[top--]]=i+1;
for(int i=n+1; i; --i) sum[L[i]]+=sum[i]+1;
for(int i=1; i<=n; ++i) delta[R[i]]+=delta[i]+1;
for(int i=1; i<=n; ++i) sum2[i]=sum2[i-1]+sum[i]-delta[i];
LL ans=0;
for(int i=1; i<=n; ++i) ans+=sum2[i]*i%mod;
printf("%lld
",ans);
}
return 0;
}
C 我才不是萝莉控呢(DP 贪心 哈夫曼树)
显然可以(n^2)DP。我们可以从((1,1))DP到((n,1))这样DP就只有两种转移了（不需要考虑下取整），即(f[i][j]=min{f[i-1][j+1], f[i][(j+1)/2]})。
如果你熟悉哈夫曼树，会发现这就类似哈夫曼树的DP转移。（哈夫曼树的DP方法见这儿，也是(n^2)的）
不是类似，把哈夫曼树的DP过程倒过来就是这个题的走路方式了。
所以本题等价于求哈夫曼树。用堆就可以(O(nlog n))解决了。
因为每次合并出的东西单调的，所以二叉哈夫曼树也可以(O(n))解决（见代码）。
//48ms 1552kb
#include <queue>
#include <cstdio>
#include <cctype>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=1e5+5;
int n,tot,pa,pb,A[N],B[N<<1];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline int Get()
{
return (pb>tot||A[pa]<=B[pb])?A[pa++]:B[pb++];
}
int main()
{
for(int T=read(); T--; )
{
n=read(),tot=0,pa=pb=1;
for(int i=n; i; --i) A[i]=read();
LL ans=0; A[n+1]=0x7fffffff;
for(int i=1,x,y; i<n; ++i)
x=Get(),y=Get(),B[++tot]=x+y,ans+=x+y;
printf("%lld
",ans);
}
return 0;
}
考试代码
A
#include <map>
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
#define LIM 5000000
typedef long long LL;
const int N=5e6+5;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
bool Query(int l,int r,int L,int R)
{
if(l==L&&r==R) return 1;
if((L<l&&R>=l)||(L<=r&&R>r)||R<l||L>r||l==r) return 0;
// if(L<=l&&r<=R) return 0;
int m=(LL)l+r>>1;
return Query(l,m,L,R)||Query(m+1,r,L,R);
}
void Violence(int l,int r,int lim)
{
for(int i=r; i<=lim; ++i)
if(Query(0,i,l,r)) {printf("%d
",i); return;}
puts("-1");
}
//inline int Getpos(int x)
//{
// return x>65536?bit[x>>16]:bit[x];
//}
inline int Rand()
{
return (rand()<<16)|rand();
}
bool Solve(int L,int R,int lim)
{
// std::map<int,bool> vis;
int ans=lim+1;
for(int suc=0,mod=lim+1-R,fail=0; suc<=50000&&ans>R; ++suc)
{
int p=rand()%(ans-R)+R;
// if(vis[p])
// {
// if(++fail>=150000) break;
// continue;
// }
// vis[p]=1;
if(Query(0,p,L,R))
{
ans=std::min(ans,p);
for(int i=1; p>>i>=R; ++i)
if(Query(0,p>>i,L,R)) ans=std::min(ans,p>>i), ++suc;
else break;
}
}
if(ans<=lim)
{
printf("%d
",ans);
return 1;
}
return 0;
}/*
3
521 535 1196
608 624 3828
304 314 3650
1
990 1755 1080952540
540 555 2308
464 475 2833
*/
int main()
{
freopen("tree.in","r",stdin);
freopen("my.out","w",stdout);
// for(int i=1; i<=65536; ++i)
// for(int j=16; ~j; --j)
// if(i>>j&1) {bit[i]=j; break;}
int t=rand(); srand(t);
for(int T=read(),t=1; t<=T; ++t)
{
int L=read(),R=read(),lim=read();
if(R>lim) {puts("-1"); continue;}
if(lim<=1000) {Violence(L,R,lim); continue;}
if(!Solve(L,R,lim)) puts("-1");
}
return 0;
}