CQOI2016 day2 模拟赛总结

T1
N=pq
r=(p-1)(q-1)
ed=1(mod r)
c^d=n(mod N)
第一步rho，第二步直接算
第三步exgcd,第四步快速幂
强行算就可以了
exgcd忘开longlong 100->30
T2
蜜汁题意
读懂过后发现建字典树
然后随便维护个单调栈搞搞就OK了
时间nlogn
T3
每次考虑把最大的出堆，把次大的入堆
hash去重
这样是31*k*log的
很慢对不对，我们考虑把31搞到可持久化线段树里面
这样每次修改就只用修改一条链了= =
可是要hash去重。。。开个hashmap就好
然后常数卡不下来了QAQ bzoj只给10s不让人活啦
所以还是要好好学其他人的神奇方法。。。


day2 30+100+60
第一个30纯属脑抽

//Copyright(c)2016 liuchenrui#include<bits/stdc++.h>#define LL long longusing namespace std;LL e,N,c,p,q,n,d,r;LL ksc(LL a,LL b,LL c){LL ans=0;while(b){if(b&1LL){ans=ans+a;if(ans>=c)ans-=c;}b>>=1LL;a=a+a;if(a>=c)a-=c;}return ans;}LL ksm(LL a,LL b,LL c){LL ans=1;while(b){if(b&1LL){ans=ksc(ans,a,c);}b>>=1LL;a=ksc(a,a,c);}return ans;}LL exgcd(LL a,LL b,LL &x,LL &y){if(b==0){x=1;y=0;return a;}LL z=exgcd(b,a%b,x,y);LL t=x;x=y;y=t-a/b*y;return z;}int windows_rand(){return rand()&32767;}LL gcd(LL a,LL b){return b==0?a:gcd(b,a%b);}LL prho(LL c,LL n){LL i,x,y,k,d;i=1;x=y=windows_rand()%n;k=2;do{i++;d=gcd(n+y-x,n);if(d>1&&d<n)return d;if(i==k)y=x,k<<=1LL;x=(ksc(x,x,n)+n-c)%n;}while(y!=x);return n;}void rho(LL n){LL t=n;while(t>=n){t=prho(windows_rand()%(n-1)+1,n);}p=t;q=N/t;return;}int main(){freopen("crack.in", "r",stdin);freopen("crack.out","w",stdout);cin>>e>>N>>c;rho(N);r=(p-1)*(q-1);LL p;exgcd(e,r,d,p);d=((d%r)+r)%r;cout<<d<<" "<<ksm(c,d,N)<<endl;}

#include<bits/stdc++.h>using namespace std;void splay(int &v){char s=getchar();v=0;while(s>'9'||s<'0')s=getchar();while(s>='0'&&s<='9')v=v*10+s-'0',s=getchar();}int son[32000000][2],t[32000000],tot,cnt;int main(){freopen("route.in","r",stdin);freopen("route.out","w",stdout);int q;splay(q);for(;q--;){static char op[4];static int s[35],dl[35];scanf("%s",op+1);if(op[1]=='Q'){int a,b,c,d,l,r;splay(a),splay(b),splay(c),splay(d),splay(l),splay(r);for(int i=8 ;i>=1 ;i--)s[i]=a&1,a>>=1;for(int i=16;i>=9 ;i--)s[i]=b&1,b>>=1;for(int i=24;i>=17;i--)s[i]=c&1,c>>=1;for(int i=32;i>=25;i--)s[i]=d&1,d>>=1;int now=0,ans=0,top=0;for(int i=1;i<=32;i++){now=son[now][s[i]];if(t[now]){while(top&&t[now]<dl[top])top--;dl[++top]=t[now];}if(now==0)break;}for(int i=1;i<=top;i++){if(l<=dl[i] && dl[i]<=r)ans++;}printf("%d\n",ans);}else{cnt++;int a,b,c,d,p;splay(a),splay(b),splay(c),splay(d),splay(p);for(int i=8 ;i>=1 ;i--)s[i]=a&1,a>>=1;for(int i=16;i>=9 ;i--)s[i]=b&1,b>>=1;for(int i=24;i>=17;i--)s[i]=c&1,c>>=1;for(int i=32;i>=25;i--)s[i]=d&1,d>>=1;int now=0;for(int i=1;i<=p;i++){if(!son[now][s[i]])son[now][s[i]]=++tot;now=son[now][s[i]];}t[now]=cnt;}}}

//liuchenrui#include<bits/stdc++.h>#include<hash_map>#define LL long long#define H 1000001007ULL#define ULL unsigned long long#define ld long doubleusing namespace std;LL n;int k;bool ispri[140];int pri[140],tot;ULL h[35];int root[2000000],ls[20000000],rs[20000000],cnt;LL ji[20000010];ULL hash[20000000];int num[55];struct node{int id;friend bool operator < (const node &a,const node &b){return ji[root[a.id]]<ji[root[b.id]];}};namespace __gnu_cxx{template<> struct hash<ULL>{size_t operator()(ULL x) const{return x;}};}__gnu_cxx::hash_map<ULL,bool>s;void build(int &now,int l,int r,int pos,int v,LL f){now=++cnt;if(l==r){if(l==pos)hash[now]=v,ji[now]=f;else ji[now]=1;return;}int mid=l+r>>1;build(ls[now],l,mid,pos,v,f),build(rs[now],mid+1,r,pos,v,f);hash[now]=hash[ls[now]]*(h[r-mid])+hash[rs[now]];ji[now]=ji[ls[now]]*ji[rs[now]];}void getval(int now,int l,int r){if(l==r){num[l]=hash[now];return;}int mid=l+r>>1;getval(ls[now],l,mid);getval(rs[now],mid+1,r);}void insert(int pre,int &now,int l,int r,int pos,int v){if(!now)now=++cnt;if(l==r){hash[now]=v;LL t=1;for(int i=1;i<=v;i++){t*=(LL)pri[pos];}ji[now]=t;return;}int mid=l+r>>1;if(pos<=mid){if(ls[pre]==ls[now])ls[now]=0;insert(ls[pre],ls[now],l,mid,pos,v);if(!rs[now])rs[now]=rs[pre];}else{if(rs[pre]==rs[now])rs[now]=0;insert(rs[pre],rs[now],mid+1,r,pos,v);if(!ls[now])ls[now]=ls[pre];}hash[now]=hash[ls[now]]*(h[r-mid])+hash[rs[now]];ji[now]=ji[ls[now]]*ji[rs[now]];}ULL gethash(ULL prehash,int pos){prehash-=h[31-pos];if(pos!=1)prehash+=h[32-pos];return prehash;}int main(){freopen("smooth.in","r",stdin);freopen("smooth.out","w",stdout);cin>>n>>k;for(int i=2;i<=128;i++){if(!ispri[i]){for(int j=i+i;j<=128;j+=i){ispri[j]=1;}pri[++tot]=i;}}h[0]=1;for(int i=1;i<=32;i++){h[i]=h[i-1]*H;}priority_queue<node>que;for(int i=1;i<=31;i++){ld x=1;int lim;for(int j=1;;j++){x*=(ld)pri[i];if(x>n){lim=j-1;break;}}LL y=1;for(int j=1;j<=lim;j++){y*=(LL)pri[i];}build(root[i],1,31,i,lim,y);que.push((node){i});s[hash[root[i]]]=1;}int rcnt=31;for(int T=1;T<k;T++){int f=que.top().id;que.pop();ULL nh=hash[root[f]];getval(root[f],1,31);for(int i=1;i<=31;i++){if(num[i]&&!s[gethash(nh,i)]){s[gethash(nh,i)]=1;++rcnt;insert(root[f],root[rcnt],1,31,i,num[i]-1);if(i!=1)insert(root[f],root[rcnt],1,31,i-1,num[i-1]+1);que.push((node){rcnt});}}}printf("%lld\n",ji[root[que.top().id]]);cerr<<clock()<<endl<<rcnt<<endl;}

最后那个会T，所以你懂得