USTC April AB题

 

USTC April 解题报告

【解题过程】

第一题看了A，直接套了整数求根的模板……其实不是很了解基本的原理，周末看一看。如果自己写的话应该是二分查找和高精度乘法。

A过了以后看B通过得人数很多，应该是水题，看了一下没什么问题就直接写了。系统返回Restricted Function若干次，（去年有次比赛也是在科大的OJ上碰到这个结果）于是包含各种库，再把函数名都加了”usr_”还是RF。最后看了faq才明白RE也可能判成Restricted Function，再读一下数据范围，少了一个0，加一下就过了。本来是可以1Y的代码折腾了这么多次，导致罚时直接罚到同题数倒数。

看了G，方向想偏了，没有做。

写了比较长时间的D题，最后还是没有过。

 

【思路&代码】

A

源代码：

 #include <cstdio>#include <cstdlib>#include <iostream>#include <cstring>using namespace std;#define iszero(t) (t.len==1&&t.s[0]==0)#define setlen(l,t) t.len=l;while(t.len>1&&t.s[t.len-1]==0) t.len--const int maxlen=500;struct bigint{ int len,s[maxlen]; bigint() {*this=0;} bigint(int a) {*this=a;} bigint(const char *a) {*this=a;} bigint operator=(int); bigint operator=(const char*); bigint operator-(const bigint&); bigint operator*(const bigint&); static int cmp(const bigint&,const bigint&); static bigint sqrt(const bigint&); //Require - * cmp};bigint bigint::operator=(int a){ len=0; do{s[len++]=a%10;a/=10;} while(a>0); return *this;}bigint bigint::operator=(const char *a){ len=strlen(a); for(int i=0;i<len;i++) s[i]=a[len-i-1]-'0'; return *this;}//Make sure *this>=abigint bigint::operator-(const bigint &a){ bigint b; for(int i=0;i<len;i++) b.s[i]=s[i]-(i<a.len?a.s[i]:0); for(int i=0;i<len;i++) if(b.s[i]<0) {b.s[i]+=10;b.s[i+1]--;} setlen(len,b); return b;}bigint bigint::operator*(const bigint &a){ bigint b; memset(b.s,0,sizeof(*s)*(len+a.len+1)); for(int i=0;i<len;i++) for(int j=0;j<a.len;j++) b.s[i+j]+=s[i]*a.s[j]; for(int i=0;i<len+a.len;i++) {b.s[i+1]+=b.s[i]/10;b.s[i]%=10;} setlen(len+a.len+1,b); return b;}int bigint::cmp(const bigint &a,const bigint &b){ if(a.len<b.len) return -1; else if(a.len>b.len) return 1; for(int i=a.len-1;i>=0;i--) if(a.s[i]!=b.s[i]) return a.s[i]-b.s[i]; return 0;}bigint bigint::sqrt(const bigint &a){ int n=(a.len-1)/2,p; bigint b,d; b.len=n+1; for(int i=n;i>=0;i--) { if(!iszero(d)) { for(int j=d.len+1;j>1;j--) d.s[j]=d.s[j-2]; d.s[0]=a.s[i*2];d.s[1]=a.s[i*2+1]; d.len+=2; } else d=a.s[i*2]+(i*2+1<a.len?a.s[i*2+1]*10:0); bigint c; c.s[1]=0; for(int j=1;j<=n-i;j++) { c.s[j]+=b.s[i+j]<<1; if(c.s[j]>=10) {c.s[j+1]=1;c.s[j]-=10;} else c.s[j+1]=0; } c.len=n-i+1+c.s[n-i+1]; for(p=1;;p++) { c.s[0]=p; if(cmp(d,c*p)<0) break; } b.s[i]=c.s[0]=p-1; d=d-c*(p-1); } return b;}int main(){ bigint n, ans1, ansori; int t, i, j; char temp[maxlen]; scanf("%d", &t); while (t--){ scanf("%s", temp); //整根，判断小数点用 ansori=temp; ansori=ansori.sqrt(ansori); //末尾加0 for (i=0; temp[i]!='/0'; i++) ; for (j=200; j>=1; j--) temp[i++]='0'; temp[i]='/0'; n=temp; //求根 ans1=n.sqrt(n); //输出 for (i=ans1.len-1, j=0; i>=2; i--, ansori.len--){ printf("%d", ans1.s[i]); if (++j==99) break; if (ansori.len==1) printf("."); } printf("/n"); } return 0;}  

B

根据数据范围其实可以很容易判断出靠谱算法大概的时间和空间复杂度。从”0<=n<=10^6”可以推知大概有O(nlogn)的解法，于是就可以想到二分了。

二分一下结果的串长，然后用一次扫描1~n所有数字。扫描过程中用flag[i]记录一下数字i最后出现的位置。（读入时进行预处理，把不属于1..m范围的数字都改成0）

源代码：

#include <cstdio>#include <cstdlib>#define maxlen 1000010int n,m;int flag[5010];int a[maxlen];bool usrinit(){ scanf("%d%d", &m, &n); if (m==0 && n==0) return false; for (int i=0; i<n; i++){ scanf("%d", &a[i]); if (a[i]<0 || a[i]>m) a[i]=0; } return true;}int usrcheck(int k){ int num; //return rightmost index (if not found return 0) for (int i=0; i<=m; i++) flag[i]=-1; num=0; //init for (int i=0; i<k; i++){ if (a[i]!=0 && flag[a[i]]==-1) num++; flag[a[i]]=i; } if (num==m) return k; for (int i=k; i<n; i++){ int j=i-k; if (a[j]!=0 && flag[a[j]]<=j){num--; flag[a[j]]=-1;} if (a[i]!=0 && flag[a[i]]==-1){num++; flag[a[i]]=i;} else flag[a[i]]=i; if (num==m) return i+1; } return 0;}void usrfind(){ if (m==0 || n<m){ printf("-1/n"); return; } int l=m, r=n, mid, idx; int lastidx, lastok=-1; while (l<=r){ mid=(l+r)>>1; if ((idx=usrcheck(mid))!=0){ lastok=mid; lastidx=idx; r=mid-1; } else l=mid+1; } if (lastok==-1) printf("-1/n"); else printf("%d %d/n", lastok, lastidx);}int main(){ while (usrinit()){ usrfind(); } return 0;}