[考试反思]1005csp-s模拟测试61：休止

连续不知道多少场了，都是一场10名以内一场20以外。。。波动极大。。。还极有规律。。。

拿到这套题，看到T1大模拟无话可说。

然后考场上我觉得T2很简单。。。。然后就码了两个半小时。

T3数据水了暴力70。。。

T1：砖块

大模拟。其实也不大。。。

记录上下左右前后边界然后滚就是了

闲得慌打hash_map

 #include<cstdio>
 #include<iostream>
 using namespace std;
 #define P(a,b) make_pair(a,b)
 struct hash_map{
     int fir[],l[],tox[],toy[],v[],cnt;
     #define mod 100003
     int& operator[](pair<int,int>p){
         int x=p.first,y=p.second,hsh=(x*y%mod+mod)%mod;
         for(int i=fir[hsh];i;i=l[i])if(tox[i]==x&&toy[i]==y)return v[i];
         l[++cnt]=fir[hsh];fir[hsh]=cnt;tox[cnt]=x;toy[cnt]=y;v[cnt]=;return v[cnt];
     }
     void clear(){for(int i=;i<mod;++i)fir[i]=;cnt=;}
 }m;
 int mxt,u,l,r,f,b;char s[];
 void n_turn(){
     int U=u,L=l,R=r,F=f,B=b;
     u=B-F;f=B;b=B+U;
     for(int i=l;i<r;++i)for(int j=f;j<b;++j){int &x=m[P(i,j)];x++,mxt=max(mxt,x);}
 }
 void s_turn(){
     int U=u,L=l,R=r,F=f,B=b;
     u=B-F;b=F;f=F-U;
     for(int i=l;i<r;++i)for(int j=f;j<b;++j){int &x=m[P(i,j)];x++,mxt=max(mxt,x);}
 }
 void w_turn(){
     int U=u,L=l,R=r,F=f,B=b;
     u=R-L;r=L;l=L-U;
     for(int i=l;i<r;++i)for(int j=f;j<b;++j){int &x=m[P(i,j)];x++,mxt=max(mxt,x);}
 }
 void e_turn(){
     int U=u,L=l,R=r,F=f,B=b;
     u=R-L;l=R;r=R+U;
     for(int i=l;i<r;++i)for(int j=f;j<b;++j){int &x=m[P(i,j)];x++,mxt=max(mxt,x);}
 }
 int main(){//freopen("ex_block2.in","r",stdin);
     int t;scanf("%d",&t);
     while(t--){
         l=,r=,f=,b=,mxt=;scanf("%d%s",&u,s);m.clear();m[P(,)]++;
         for(int i=;s[i];++i)
             if(s[i]=='N')n_turn();
             else if(s[i]=='E')e_turn();
             else if(s[i]=='W')w_turn();
             else if(s[i]=='S')s_turn();
         for(int i=l;i<r;++i)for(int j=f;j<b;++j)printf("%d ",i);puts("");
         for(int i=l;i<r;++i)for(int j=f;j<b;++j)printf("%d ",j);puts("");
         printf("%d\n",mxt);
     }
 }

T2：数字

比较神仙。

高精还是要打的，但是只需要读入，除2和5的倍数的低精数，模同理，以及大于号。。。

模的话直接取最低位返回就行。。。emm我是万进制的

设Ext(a,b)是b除去所有因子a后剩下的数，那么题目就是求Ext(10,n!)%10^k

直觉，可以CRT。分成2和5处理。

绝大多数情况下因子2都严格比5多，所以Ext(10,n!)%2==0

现在只需要求出Ext(5,n!)%10^k

那么分解1000就是8×125。对于8其因子数也很多所以最后乘逆元。

问题就变成求Ext mod125的答案了。

可以把阶乘里能整除5的和不能的分开考虑，而能被整除的部分就是Ext(5,(n/5)!),递归处理。

可以发现Ext之后阶乘操作其实是循环的，因为mod125所以每125位就循环了。

循环的部分可以快速幂处理。剩余部分可以预处理。

我知道这作为一个题解过于草率，没时间了溜了。

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 using namespace std;
 char s[];int l,k;const int mod[]={,,,},m5[]={,,,};int fac[],FAC[];
 struct bigint{
     int a[],ws;
     void get(){
         ws=(l-)/+;for(int i=;i<=;++i)a[i]=;
         for(int i=;i<ws;++i)a[i]=s[l-*i]*+s[l-*i+]*+s[l-*i+]*+s[l-*i+];
         for(int i=;i<=(l-)%+;++i)a[ws]=a[ws]*+s[i-];
     }
     void operator/=(int k){
         for(int i=ws;i;--i)a[i-]+=a[i]%k*,a[i]/=k;
         while(!a[ws]&&ws)ws--;
     }
     int operator%(int k){return a[]%k;}
     bool operator>(int k){return ws>||a[]>k;}
 }n,m;
 int pow(int b,int t,int mod,int a=){for(;t;t>>=,b=b*b%mod)if(t&)a=a*b%mod;return a;}
 int Ext(bigint n){
     if(n>);else return fac[n.a[]];
     int x=n%;n/=;int K=Ext(n);n/=;
     return K*pow(FAC[],n%,)%*FAC[x]%;
 }
 void put(int x,int k){for(int i=k;i;--i)if(x%mod[i]/mod[i-]==)putchar();else break;printf("%d\n",x%mod[k]);}
 int main(){//freopen("t2.in","r",stdin);freopen("t2.out","w",stdout);
     fac[]=FAC[]=;for(int i=;i<=;++i){int j=i;while(j%==)j/=;fac[i]=fac[i-]*j%;FAC[i]=FAC[i-]*(i%==?:i)%;}
     int ans,t,tot;scanf("%d",&t);
     while(t--){
         scanf("%s%d",s,&k);l=strlen(s);
         for(int i=;i<l;++i)s[i]-=;
         n.get();ans=Ext(n);
         if(l==&&n.a[]<=){int a=;for(int i=;i<=n.a[];++i)a*=i;if(a%==)a/=;put(a,k);continue;}
         m=n;tot=;
         while(m.ws)m/=,(tot+=m%)%=;
         ans*=pow(,tot,);ans%=;
         for(int i=;i<;i+=)if(i%==ans)put(i,k);
     }
 }

T3：甜圈

考察的思路很不错啊。

因为不能多不能少不能乱序，所以很符合字符串的性质。

因为要区间操作，所以方便快捷的hash。

那么就变成了区间加区间乘的线段树啦！

（居然还有人不会线段树板子%%%）

 #include<cstdio>
 int cl[],cr[],ans;
 unsigned long long hsh[],lzm[],lza[],HSH;
 void build(int p,int l,int r){
     cl[p]=l;cr[p]=r;lzm[p]=;
     if(l==r)return;
     build(p<<,l,l+r>>);build(p<<|,(l+r>>)+,r);
 }
 void down(int p){
     hsh[p<<]*=lzm[p];hsh[p<<|]*=lzm[p];
     lzm[p<<]*=lzm[p];lzm[p<<|]*=lzm[p];
     lza[p<<]*=lzm[p];lza[p<<|]*=lzm[p];
     hsh[p<<]+=lza[p];hsh[p<<|]+=lza[p];
     lza[p<<]+=lza[p];lza[p<<|]+=lza[p];
     lza[p]=;lzm[p]=;
 }
 void chg(int p,int l,int r,int w){
     if(l<=cl[p]&&cr[p]<=r){hsh[p]=hsh[p]*+w;lza[p]*=;lzm[p]*=;lza[p]+=w;return;}
     if(lza[p]||lzm[p]!=)down(p);
     if(l<=cr[p<<])chg(p<<,l,r,w);
     if(r>=cl[p<<|])chg(p<<|,l,r,w);
 }
 void ask(int p){
     if(cl[p]==cr[p]){ans+=HSH==hsh[p];return;}
     if(lza[p]||lzm[p]!=)down(p);
     ask(p<<);ask(p<<|);
 }
 int main(){
     int n,t,l,r,w,k;
     scanf("%d%d%d",&n,&k,&t);
     build(,,n);
     while(t--)scanf("%d%d%d",&l,&r,&w),chg(,l,r,w);
     for(int i=;i<=k;++i)HSH=HSH*+i;
     ask();printf("%d\n",ans);
 }