機房測試10.6
阿新 • • 發佈：2018-10-07
stdout == ash 學習 spa unsigned cst 構造 函數
磚塊
很簡單的水題，寫起來也很簡單。
模擬寫的很短，是否也是碼力的體現？
#include<cstdio> #include<cstring> #include<map> #define FN "block" int mp[1105][1105],dao,ans; char s[105]; int x,y,dx[4]={0,0,-1,1},dy[4]={1,-1,0,0}; std::map<char,int> m; void init() { memset(mp,0,sizeof(mp)); mp[550][550]=1; ans=x=y=0; dao=-1; } int main() { freopen(FN".in","r",stdin); freopen(FN".out","w",stdout); m[‘N‘]=0,m[‘S‘]=1,m[‘W‘]=2,m[‘E‘]=3; int T;scanf("%d",&T); while(T--) { init(); int k;scanf("%d",&k); scanf("%s",s+1); int len=strlen(s+1); for(int i=1;i<=len;i++) { x+=dx[m[s[i]]]; y+=dy[m[s[i]]]; if(!~dao) { for(int j=0;j<k;j++) if(++mp[x+dx[m[s[i]]]*j+550][y+dy[m[s[i]]]*j+550]>ans) ans=mp[x+dx[m[s[i]]]*j+550][y+dy[m[s[i]]]*j+550]; dao=m[s[i]]; } else if((dao^m[s[i]])==1) { if(++mp[x+550][y+550]>ans) ans=mp[x+550][y+550]; dao=-1; } else if(dao==m[s[i]]) { x+=dx[m[s[i]]]*(k-1),y+=dy[m[s[i]]]*(k-1); if(++mp[x+550][y+550]>ans) ans=mp[x+550][y+550]; dao=-1; } else { for(int j=0;j<k;j++) if(++mp[x+dx[dao]*j+550][y+dy[dao]*j+550]>ans) ans=mp[x+dx[dao]*j+550][y+dy[dao]*j+550]; } } if(!~dao) printf("%d\n%d\n",x,y); else if(dao==0 || dao==3) { for(int j=0;j<k;j++) printf("%d%c",x+dx[dao]*j,j==k-1?‘\n‘:‘ ‘); for(int j=0;j<k;j++) printf("%d%c",y+dy[dao]*j,j==k-1?‘\n‘:‘ ‘); } else { for(int j=k-1;j>=0;j--) printf("%d%c",x+dx[dao]*j,j?‘ ‘:‘\n‘); for(int j=k-1;j>=0;j--) printf("%d%c",y+dy[dao]*j,j?‘ ‘:‘\n‘); } printf("%d\n",ans); } return 0; }
可以不用map，但是我懶得改了。
翻轉有好幾種情況，只要去想，還是不會漏的（如果漏了過不了樣例）
數字
真的堪稱神仙題，題解至今沒有看懂。
威爾遜定理加上balabala可以解決\(k=1\)的情況。
預處理1e7階乘%1e7可以50分，%1e5只有20分。
std的代碼（有高精度）：
#include<cstdio> #include<cstdlib> #include<cstring> #include<algorithm> using namespace std; typedef long long ll; const int con=100000000; const int wei=8; class Int{ public:long long a[100]; void read(){ memset(a,0,sizeof(a)); char S[105]={0}; int len=0; scanf("%s",S+1); len=strlen(S+1); for (int i=1;i<=len/2;i++) swap(S[i],S[len-i+1]); for (int i=1;i<=len;i++) S[i]-=‘0‘; a[0]=0; for (int i=1;i<=len;i+=8){ a[0]++; for (int j=i+7;j>=i;j--) a[a[0]]=a[a[0]]*10+S[j]; } } void getdata(int x){memset(a,0,sizeof(a));while (x){a[++a[0]]=x%con;x=x/con;}} Int(int x=0){ getdata(x); } void pri(bool flag){ if (a[0]==0||(a[0]==1&&a[1]==0)){printf("0");if (flag)printf("\n");return;} printf("%lld",a[a[0]]); for (int i=a[0]-1;i;i--) printf("%08lld",a[i]); if (flag)printf("\n"); } bool operator <(const Int &X){ if (a[0]<X.a[0])return true;if (a[0]>X.a[0])return false; for (int i=a[0];i;i--){if (a[i]<X.a[i])return true;if (a[i]>X.a[i])return false;} return false; } bool operator >(const Int &X){ if (a[0]<X.a[0])return false;if (a[0]>X.a[0])return true; for (int i=a[0];i;i--){if (a[i]<X.a[i])return false;if (a[i]>X.a[i])return true;} return false; } bool operator <=(const Int &X){ if (a[0]<X.a[0])return true;if (a[0]>X.a[0])return false; for (int i=a[0];i;i--){if (a[i]<X.a[i])return true;if (a[i]>X.a[i])return false;} return true; } bool operator >=(const Int &X){ if (a[0]<X.a[0])return false;if (a[0]>X.a[0])return true; for (int i=a[0];i;i--){if (a[i]<X.a[i])return false;if (a[i]>X.a[i])return true;} return true; } bool operator ==(const Int &X){ if (a[0]!=X.a[0])return false;for (int i=a[0];i;i--)if (a[i]!=X.a[i])return false; return true; } Int operator +(const Int &X){ Int c;memset(c.a,0,sizeof(c.a)); for (int i=1;i<=a[0]||i<=X.a[0];i++) {c.a[i]=c.a[i]+a[i]+X.a[i];c.a[i+1]+=c.a[i]/con;c.a[i]%=con;} c.a[0]=max(a[0],X.a[0]);if (c.a[c.a[0]+1])c.a[0]++; return c; } Int operator +(int num){ Int c;memcpy(c.a,a,sizeof(c.a));c.a[1]+=num; for (int i=1;i<=c.a[0]&&c.a[i]>=con;i++)c.a[i]-=con,c.a[i+1]++; while (c.a[c.a[0]+1])c.a[0]++; return c; } Int operator -(const Int &X){ Int c;memcpy(c.a,a,sizeof(c.a)); for (int i=1;i<=a[0];i++){c.a[i]=c.a[i]-X.a[i];if (c.a[i]<0){c.a[i+1]--;c.a[i]+=con;}} while (c.a[0]&&!c.a[c.a[0]])c.a[0]--; return c; } Int operator -(int num){ Int c;memcpy(c.a,a,sizeof(c.a));c.a[1]-=num; for (int i=1;i<=c.a[0]&&c.a[i]<0;i++)c.a[i]+=con,c.a[i+1]--; while (c.a[0]&&!c.a[c.a[0]])c.a[0]--; return c; } Int operator *(const Int &X){ Int c;memset(c.a,0,sizeof(c.a)); for (int i=1;i<=a[0];i++)for (int j=1;j<=X.a[0];j++) {c.a[i+j-1]+=a[i]*X.a[j];c.a[i+j]+=c.a[i+j-1]/con;c.a[i+j-1]%=con;} c.a[0]=max(a[0]+X.a[0]-1,0ll);if (c.a[a[0]+X.a[0]]>0)c.a[0]++; return c; } Int operator *(int num){ Int c;memset(c.a,0,sizeof(c.a)); for (int i=1;i<=a[0];i++){c.a[i]+=a[i]*num;if (c.a[i]>=con){c.a[i+1]+=c.a[i]/con;c.a[i]%=con;}} c.a[0]=a[0];if (c.a[c.a[0]+1]>0)c.a[0]++; return c; } Int operator /(int num){ Int c;memset(c.a,0,sizeof(c.a)); long long x=0;for (int i=a[0];i;i--){x=x*con+a[i];c.a[i]=x/num;x=x%num;} c.a[0]=a[0];if (c.a[0]&&!c.a[c.a[0]])c.a[0]--; return c; } int operator %(int num){ int ret=0; for (int i=a[0];i;i--) ret=((ll)ret*con+a[i])%num; return ret; } }; inline int Pow(int a,int b,int p){ int ret=1; for (int i=1;i<=b;i++) ret=ret*a%p; return ret; } Int N,K,K0; int k; int S[1005],m; inline void Pre(){ Int t=N; K0.getdata(0); while (!(t==0)) K0=K0+t/5,t=t/5; } int ans[200][200]; int a[200]; inline int Solve(){ int phi=k==1?4:(k==2?4*5:4*5*5); int mod2=0,mod5,k2=Pow(2,k,2333),k5=Pow(5,k,2333); for (int i=0;i<k2*k5;i++) ans[i%k2][i%k5]=i; int inv=Pow(2,phi-1,k5),Inv=Pow(inv,K0%phi,k5); a[0]=1; for (int i=1;i<k5;i++){ int t=i; if (t%5==0) t=1; a[i]=(a[i-1]*t)%k5; } Int T=N; mod5=1; while (!(T==0)){ mod5=mod5*Pow(a[k5-1],(T/k5)%phi,k5)%k5; mod5=mod5*a[T%k5]%k5; T=T/5; } mod5=mod5*Inv%k5; return ans[mod2][mod5]; } inline void print(long long t,int k){ if (k==1) printf("%lld\n",t%10); else if (k==2) printf("%02lld\n",t%100); else if (k==3) printf("%03lld\n",t%1000); } int main(){ freopen("num.in","r",stdin); freopen("num.out","w",stdout); int T; scanf("%d",&T); while (T--){ N.read(); scanf("%d",&k); if (N<=15){ long long t=1; for (int i=1;Int(i)<=N;i++) t*=i; while (t%10==0) t/=10; print(t,k); }else{ Pre(); print(Solve(),k); } } return 0; }
不寫高精度有80分，但至今不知怎麽做。
暴力分給足了的。
因為合法的順序對應序列是唯一的，於是進行Hash，操作後如果對它所有操作的Hash值與正確的相同，那麽就正確了。
所以用兩個Tag的線段樹維護Hash值的加法和乘法。
std寫的很妙，我學習了一波。
#include<cstdio> #include<cstdlib> #include<algorithm> using namespace std; typedef unsigned int uint; inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++; } inline void read(int &x){ char c=nc(),b=1; for (;!(c>=‘0‘ && c<=‘9‘);c=nc()) if (c==‘-‘) b=-1; for (x=0;c>=‘0‘ && c<=‘9‘;x=x*10+c-‘0‘,c=nc()); x*=b; } const int N=200005; const uint seed=233333; struct abcd{ uint a,b; abcd(uint a=1,uint b=0):a(a),b(b) { } bool one() { return a==1 && b==0; } void add(abcd B){ uint _a=a,_b=b; a=B.a*_a; b=B.a*_b+B.b; } }T[N<<2]; inline void modify(int x,int l,int r,int ql,int qr,abcd t){ if (ql<=l && r<=qr){ T[x].add(t); return; } if (!T[x].one()) T[x<<1].add(T[x]),T[x<<1|1].add(T[x]),T[x]=abcd(); int mid=(l+r)>>1; if (ql<=mid) modify(x<<1,l,mid,ql,qr,t); if (qr>mid) modify(x<<1|1,mid+1,r,ql,qr,t); } uint Hash[N]; inline void query(int x,int l,int r){ if (l==r){ Hash[l]=T[x].b; return; } if (!T[x].one()) T[x<<1].add(T[x]),T[x<<1|1].add(T[x]),T[x]=abcd(); int mid=(l+r)>>1; query(x<<1,l,mid); query(x<<1|1,mid+1,r); } int n,K; int main(){ int T,l,r,x; freopen("deco.in","r",stdin); freopen("deco.out","w",stdout); read(n); read(K); read(T); while (T--){ read(l); read(r); read(x); modify(1,1,n,l,r,abcd(seed,x)); } query(1,1,n); uint h=0; for (int i=1;i<=K;i++) h=h*seed+i; int ans=0; for (int i=1;i<=n;i++) if (Hash[i]==h) ans++; printf("%d\n",ans); return 0; }
用結構體存兩個Tag，用構造函數運算出值。
妙啊~
機房測試10.6