遍地都是坑,一个一个填,由于本人懒数据结构弱的一笔所以就慢慢填
由于日常%了Claris,于是得到了Claris的课件,题都是课件里找的_(:зゝ∠)_
你要课件?你tm真的以为我会给你吗?
现在做了多少: 10/10
$3.7$ $qiancl:$喵的填坑填辣么慢
$3.10$ $qiancl:$嘿嘿嘿填完了
遍地都是坑,一个一个填,由于本人懒数据结构弱的一笔所以就慢慢填
由于日常%了Claris,于是得到了Claris的课件,题都是课件里找的_(:зゝ∠)_
你要课件?你tm真的以为我会给你吗?
现在做了多少: 10/10
$3.7$ $qiancl:$喵的填坑填辣么慢
$3.10$ $qiancl:$嘿嘿嘿填完了
之前的坑一直没填,因为懒找不到题
当然坑还是要填的,但是由于本人懒实在不太能找到kd-tree的题了于是把坑从10道减到了5道_(:зゝ∠)_
先放出最后一题的链接吧,慢慢填
http://www.lydsy.com/JudgeOnline/problem.php?id=3815
我有一种预感,这坑要填很久
现在: 4/5
传送门-->http://www.lydsy.com/JudgeOnline/problem.php?id=3944
本来想留着给初三爷出题。。然而都没几个会欧拉函数的,所以就写个题解等我忘了看……
最开始给我灵感的是codevs月赛的一道题,它的一个子问题是求
$\sum_{i=1}^{n}\sigma \left ( i \right )$,其中$\sigma \left ( n \right )= \sum_{d|n}d$
介个东西怎么做呢,当然并不能直接做喽,当时比利直接秒出%%%,he said :"This is a foolish thing." 他是考虑每个数对答案的贡献,我还是暴力推了推公式
$\sum_{i=1}^{n} \sigma \left ( i \right )= \sum_{i=1}^{n}\sum_{j=1}^{n}[j|i]\cdot j= \sum_{i=1}^{n}i\cdot [\frac{n}{i}]$
传送门-->http://www.lydsy.com/JudgeOnline/problem.php?id=3545
sb题
离线做,把边和询问都排序,然后一条条加边,维护连通块,于是这就是个平衡树合并裸题。。然而只有我这种都比被卡了辣么久- =
一发5s,感觉跑好慢……,一看榜,发现rank3而且代码也很短2333
#include<cstring> #include<cstdio> #include<algorithm> using namespace std; #define nc() getchar() inline int read(){ int x=0;char ch=nc();for(;ch<'0'||ch>'9';ch=nc()); for(;ch<='9'&&ch>='0';x=x*10+ch-48,ch=nc());return x; } #define N 100010 int i,j,k,m,n,x,y,Q,h[N],cnt,fa[N],pp[N],ans[N*5],top; struct edge{ int x,y,z; inline void in(){x=read(),y=read(),z=read();} bool operator <(const edge&a)const{return z<a.z;} }e[N*5]; struct ques{ int v,x,k,id; inline void in(){v=read(),x=read(),k=read(),id=i;} bool operator <(const ques&a)const{return x<a.x;} }q[N*5]; int find(int x){return x==fa[x]?x:fa[x]=find(fa[x]);} struct node{ node*ch[2];int r,sz,h,id;inline void up(){sz=ch[0]->sz+ch[1]->sz+1;} inline int cmp(int x)const{return x>h;} }C[N],Null,*null=&Null; inline void rotate(node*&o,int d){node*k=o->ch[d^1];o->ch[d^1]=k->ch[d],k->ch[d]=o,k->sz=o->sz,o->up(),o=k;} int insert(node*&o,int x,int id){ if(o==null)o=&C[id],o->id=id,o->h=x,o->ch[0]=o->ch[1]=null,o->sz=1,o->r=rand(); else{ int d=o->cmp(x);++o->sz; insert(o->ch[d],x,id);if(o->ch[d]->r>o->r)rotate(o,d^1); } return o->id; } void Add(node*x,node*y){if(x==null)return;Add(x->ch[0],y),Add(x->ch[1],y),top=insert(y=&(C[top]),x->h,x->id);} inline void link(int u,int v){ if(u==v)return; if((C+u)->sz<(C+v)->sz)swap(u,v); fa[v]=u,top=pp[u],Add(C+pp[v],C+pp[u]),pp[u]=top; } int kth(node*o,int k){ if(o==null||k<=0||k>o->sz)return 0x7f7f7f7f; if(k<=o->ch[0]->sz)return kth(o->ch[0],k); if(k>o->ch[0]->sz+1)return kth(o->ch[1],k-o->ch[0]->sz-1); return o->h; } int main(){ for(n=read(),m=read(),Q=read(),i=1;i<=n;++i)h[i]=read(),fa[i]=i,pp[i]=i,(C+i)->id=i,(C+i)->h=h[i],(C+i)->sz=1,(C+i)->ch[0]=(C+i)->ch[1]=null,(C+i)->r=rand(); for(i=1;i<=m;++i)e[i].in();for(i=1;i<=Q;++i)q[i].in(); for(cnt=1,sort(e+1,e+m+1),sort(q+1,q+Q+1),i=1;i<=Q;++i){ while(e[cnt].z<=q[i].x&&cnt<=m)link(find(e[cnt].x),find(e[cnt].y)),++cnt; ans[q[i].id]=kth(C+pp[find(q[i].v)],(C+pp[find(q[i].v)])->sz-q[i].k+1); } for(i=1;i<=Q;++i)if(ans[i]==0x7f7f7f7f)puts("-1");else printf("%d\n",ans[i]); }
有人问我为什么有些题那么快_(:зゝ∠)_
所以我来解释一下我加的黑科技,题目编号是乱的= =,有想要我解释的题请评论
首先是fastio,bzoj上有了这个就可以轻易rank1吧。。
#define S 100000 char bf[S],*p1=bf,*p2=bf; #define nc() (p1==p2&&(p2=(p1=bf)+fread(bf,1,S,stdin),p2==p1)?-1:*p1++) inline int read(){ int x=0;char ch=nc();for(;ch<'0'||ch>'9';ch=nc()); for(;ch<='9'&&ch>='0';x=x*10+ch-48,ch=nc());return x; } char pf[S],*o1=pf,*o2=pf+S; #define ot(x) (o1==o2?fwrite(pf,1,S,stdout),*(o1=pf)++=x:*o1++=x) inline void print(int x){static char s[15],*b;b=s;if(!x)*b++=48;for(;x;*b++=x%10+48,x/=10);for(;b--!=s;ot(*b));} fwrite(pf,1,o1-pf,stdout); //for(i=1;i<=m;++i)print(ans[i]),ot('\n');
所以在这里我就不解释用fastio艹上去的题了
这个题解我是在你们做这套题的时候写的♪(^∇^*)
先讲第二题吧-。- 原题是bzoj4080
首先如果有人这题没分。。
。。。。。。。。。。。。。。
好,我们来看前60%的数据,$n<=15$,这个东西只要$2^{15}$枚举就行了我不多说
对于100%的数据,有很多种做法-。-,不知道你们认不认识我给你们的关怀里的那个 “大力”
Host by is-Programmer.com | Power by Chito 1.3.3 beta | Theme: Aeros 2.0 by TheBuckmaker.com