A -?Sqrt Bo

簽到題啦，直接套一個大數模板搞定。

不過我現在使用的大數模板不支持直接定義并賦值，bignum a = 1 這樣是不行的。

　　
　　//查看了下代碼發現，模板中bignum的構造函數，只有這一個，也就是定義一個大數，初始化為0(num[0]非0位數)
　　inline bignum(){memset(num,0,sizeof(bignum_t));num[0]=1 ;sgn=0 ;}

有因為模板中重載了=

inline bignum&operator=(const int a){memset(num,0,sizeof(bignum_t));num[0]=1 ;sgn=SGN (a);add(num,sgn*a);return*this ;}

所以可以采用：

bignum a;

a = 1233;

這種方法來定義一個靜態大數。

// // main.cpp // multi2016.2.firtst // // Created by New_Life on 16/8/1. // Copyright ? 2016年 chenhuan001. All rights reserved. // #include <iostream> #include <math.h> #include <stdio.h> #include <algorithm> #include <string.h> using namespace std;#define DIGIT 4 //四位隔開,即萬進制 #define DEPTH 10000 //萬進制 #define MAX 30 //題目最大位數/4，要不大直接設為最大位數也行 typedef int bignum_t[MAX+1];/************************************************************************/ /* 讀取操作數，對操作數進行處理存儲在數組里 */ /************************************************************************/ int read(bignum_t a,istream&is=cin) {char buf[MAX*DIGIT+1],ch ;int i,j ;memset((void*)a,0,sizeof(bignum_t));if(!(is>>buf))return 0 ;for(a[0]=strlen(buf),i=a[0]/2-1;i>=0;i--)ch=buf[i],buf[i]=buf[a[0]-1-i],buf[a[0]-1-i]=ch ;for(a[0]=(a[0]+DIGIT-1)/DIGIT,j=strlen(buf);j<a[0]*DIGIT;buf[j++]='0');for(i=1;i<=a[0];i++)for(a[i]=0,j=0;j<DIGIT;j++)a[i]=a[i]*10+buf[i*DIGIT-1-j]-'0' ;for(;!a[a[0]]&&a[0]>1;a[0]--);return 1 ; }void write(const bignum_t a,ostream&os=cout) {int i,j ;for(os<<a[i=a[0]],i--;i;i--)for(j=DEPTH/10;j;j/=10)os<<a[i]/j%10 ; }int comp(const bignum_t a,const bignum_t b) {int i ;if(a[0]!=b[0])return a[0]-b[0];for(i=a[0];i;i--)if(a[i]!=b[i])return a[i]-b[i];return 0 ; }int comp(const bignum_t a,const int b) {int c[12]={};for(c[1]=b;c[c[0]]>=DEPTH;c[c[0]+1]=c[c[0]]/DEPTH,c[c[0]]%=DEPTH,c[0]++);return comp(a,c); }int comp(const bignum_t a,const int c,const int d,const bignum_t b) {int i,t=0,O=-DEPTH*2 ;if(b[0]-a[0]<d&&c)return 1 ;for(i=b[0];i>d;i--){t=t*DEPTH+a[i-d]*c-b[i];if(t>0)return 1 ;if(t<O)return 0 ;}for(i=d;i;i--){t=t*DEPTH-b[i];if(t>0)return 1 ;if(t<O)return 0 ;}return t>0 ; } /************************************************************************/ /* 大數與大數相加 */ /************************************************************************/ void add(bignum_t a,const bignum_t b) {int i ;for(i=1;i<=b[0];i++)if((a[i]+=b[i])>=DEPTH)a[i]-=DEPTH,a[i+1]++;if(b[0]>=a[0])a[0]=b[0];elsefor(;a[i]>=DEPTH&&i<a[0];a[i]-=DEPTH,i++,a[i]++);a[0]+=(a[a[0]+1]>0); } /************************************************************************/ /* 大數與小數相加 */ /************************************************************************/ void add(bignum_t a,const int b) {int i=1 ;for(a[1]+=b;a[i]>=DEPTH&&i<a[0];a[i+1]+=a[i]/DEPTH,a[i]%=DEPTH,i++);for(;a[a[0]]>=DEPTH;a[a[0]+1]=a[a[0]]/DEPTH,a[a[0]]%=DEPTH,a[0]++); } /************************************************************************/ /* 大數相減(被減數>=減數) */ /************************************************************************/ void sub(bignum_t a,const bignum_t b) {int i ;for(i=1;i<=b[0];i++)if((a[i]-=b[i])<0)a[i+1]--,a[i]+=DEPTH ;for(;a[i]<0;a[i]+=DEPTH,i++,a[i]--);for(;!a[a[0]]&&a[0]>1;a[0]--); } /************************************************************************/ /* 大數減去小數(被減數>=減數) */ /************************************************************************/ void sub(bignum_t a,const int b) {int i=1 ;for(a[1]-=b;a[i]<0;a[i+1]+=(a[i]-DEPTH+1)/DEPTH,a[i]-=(a[i]-DEPTH+1)/DEPTH*DEPTH,i++);for(;!a[a[0]]&&a[0]>1;a[0]--); }void sub(bignum_t a,const bignum_t b,const int c,const int d) {int i,O=b[0]+d ;for(i=1+d;i<=O;i++)if((a[i]-=b[i-d]*c)<0)a[i+1]+=(a[i]-DEPTH+1)/DEPTH,a[i]-=(a[i]-DEPTH+1)/DEPTH*DEPTH ;for(;a[i]<0;a[i+1]+=(a[i]-DEPTH+1)/DEPTH,a[i]-=(a[i]-DEPTH+1)/DEPTH*DEPTH,i++);for(;!a[a[0]]&&a[0]>1;a[0]--); } /************************************************************************/ /* 大數相乘，讀入被乘數a，乘數b，結果保存在c[] */ /************************************************************************/ void mul(bignum_t c,const bignum_t a,const bignum_t b) {int i,j ;memset((void*)c,0,sizeof(bignum_t));for(c[0]=a[0]+b[0]-1,i=1;i<=a[0];i++)for(j=1;j<=b[0];j++)if((c[i+j-1]+=a[i]*b[j])>=DEPTH)c[i+j]+=c[i+j-1]/DEPTH,c[i+j-1]%=DEPTH ;for(c[0]+=(c[c[0]+1]>0);!c[c[0]]&&c[0]>1;c[0]--); } /************************************************************************/ /* 大數乘以小數，讀入被乘數a，乘數b，結果保存在被乘數 */ /************************************************************************/ void mul(bignum_t a,const int b) {int i ;for(a[1]*=b,i=2;i<=a[0];i++){a[i]*=b ;if(a[i-1]>=DEPTH)a[i]+=a[i-1]/DEPTH,a[i-1]%=DEPTH ;}for(;a[a[0]]>=DEPTH;a[a[0]+1]=a[a[0]]/DEPTH,a[a[0]]%=DEPTH,a[0]++);for(;!a[a[0]]&&a[0]>1;a[0]--); }void mul(bignum_t b,const bignum_t a,const int c,const int d) {int i ;memset((void*)b,0,sizeof(bignum_t));for(b[0]=a[0]+d,i=d+1;i<=b[0];i++)if((b[i]+=a[i-d]*c)>=DEPTH)b[i+1]+=b[i]/DEPTH,b[i]%=DEPTH ;for(;b[b[0]+1];b[0]++,b[b[0]+1]=b[b[0]]/DEPTH,b[b[0]]%=DEPTH);for(;!b[b[0]]&&b[0]>1;b[0]--); } /**************************************************************************/ /* 大數相除,讀入被除數a，除數b，結果保存在c[]數組 */ /* 需要comp()函數 */ /**************************************************************************/ void div(bignum_t c,bignum_t a,const bignum_t b) {int h,l,m,i ;memset((void*)c,0,sizeof(bignum_t));c[0]=(b[0]<a[0]+1)?(a[0]-b[0]+2):1 ;for(i=c[0];i;sub(a,b,c[i]=m,i-1),i--)for(h=DEPTH-1,l=0,m=(h+l+1)>>1;h>l;m=(h+l+1)>>1)if(comp(b,m,i-1,a))h=m-1 ;else l=m ;for(;!c[c[0]]&&c[0]>1;c[0]--);c[0]=c[0]>1?c[0]:1 ; }void div(bignum_t a,const int b,int&c) {int i ;for(c=0,i=a[0];i;c=c*DEPTH+a[i],a[i]=c/b,c%=b,i--);for(;!a[a[0]]&&a[0]>1;a[0]--); } /************************************************************************/ /* 大數平方根，讀入大數a，結果保存在b[]數組里 */ /* 需要comp()函數 */ /************************************************************************/ void sqrt(bignum_t b,bignum_t a) {int h,l,m,i ;memset((void*)b,0,sizeof(bignum_t));for(i=b[0]=(a[0]+1)>>1;i;sub(a,b,m,i-1),b[i]+=m,i--)for(h=DEPTH-1,l=0,b[i]=m=(h+l+1)>>1;h>l;b[i]=m=(h+l+1)>>1)if(comp(b,m,i-1,a))h=m-1 ;else l=m ;for(;!b[b[0]]&&b[0]>1;b[0]--);for(i=1;i<=b[0];b[i++]>>=1); } /************************************************************************/ /* 返回大數的長度 */ /************************************************************************/ int length(const bignum_t a) {int t,ret ;for(ret=(a[0]-1)*DIGIT,t=a[a[0]];t;t/=10,ret++);return ret>0?ret:1 ; } /************************************************************************/ /* 返回指定位置的數字，從低位開始數到第b位，返回b位上的數 */ /************************************************************************/ int digit(const bignum_t a,const int b) {int i,ret ;for(ret=a[(b-1)/DIGIT+1],i=(b-1)%DIGIT;i;ret/=10,i--);return ret%10 ; } /************************************************************************/ /* 返回大數末尾0的個數 */ /************************************************************************/ int zeronum(const bignum_t a) {int ret,t ;for(ret=0;!a[ret+1];ret++);for(t=a[ret+1],ret*=DIGIT;!(t%10);t/=10,ret++);return ret ; }void comp(int*a,const int l,const int h,const int d) {int i,j,t ;for(i=l;i<=h;i++)for(t=i,j=2;t>1;j++)while(!(t%j))a[j]+=d,t/=j ; }void convert(int*a,const int h,bignum_t b) {int i,j,t=1 ;memset(b,0,sizeof(bignum_t));for(b[0]=b[1]=1,i=2;i<=h;i++)if(a[i])for(j=a[i];j;t*=i,j--)if(t*i>DEPTH)mul(b,t),t=1 ;mul(b,t); } /************************************************************************/ /* 組合數 */ /************************************************************************/ void combination(bignum_t a,int m,int n) {int*t=new int[m+1];memset((void*)t,0,sizeof(int)*(m+1));comp(t,n+1,m,1);comp(t,2,m-n,-1);convert(t,m,a);delete[]t ; } /************************************************************************/ /* 排列數 */ /************************************************************************/ void permutation(bignum_t a,int m,int n) {int i,t=1 ;memset(a,0,sizeof(bignum_t));a[0]=a[1]=1 ;for(i=m-n+1;i<=m;t*=i++)if(t*i>DEPTH)mul(a,t),t=1 ;mul(a,t); }#define SGN(x) ((x)>0?1:((x)<0?-1:0)) #define ABS(x) ((x)>0?(x):-(x))int read(bignum_t a,int&sgn,istream&is=cin) {char str[MAX*DIGIT+2],ch,*buf ;int i,j ;memset((void*)a,0,sizeof(bignum_t));if(!(is>>str))return 0 ;buf=str,sgn=1 ;if(*buf=='-')sgn=-1,buf++;for(a[0]=strlen(buf),i=a[0]/2-1;i>=0;i--)ch=buf[i],buf[i]=buf[a[0]-1-i],buf[a[0]-1-i]=ch ;for(a[0]=(a[0]+DIGIT-1)/DIGIT,j=strlen(buf);j<a[0]*DIGIT;buf[j++]='0');for(i=1;i<=a[0];i++)for(a[i]=0,j=0;j<DIGIT;j++)a[i]=a[i]*10+buf[i*DIGIT-1-j]-'0' ;for(;!a[a[0]]&&a[0]>1;a[0]--);if(a[0]==1&&!a[1])sgn=0 ;return 1 ; } struct bignum {bignum_t num ;int sgn ;public :inline bignum(){printf("!!!\n");memset(num,0,sizeof(bignum_t));num[0]=1 ;sgn=0 ;}inline int operator!(){return num[0]==1&&!num[1];}inline bignum&operator=(const bignum&a){memcpy(num,a.num,sizeof(bignum_t));sgn=a.sgn ;return*this ;}inline bignum&operator=(const int a){memset(num,0,sizeof(bignum_t));num[0]=1 ;sgn=SGN (a);add(num,sgn*a);return*this ;};inline bignum&operator+=(const bignum&a){if(sgn==a.sgn)add(num,a.num);else if(sgn&&a.sgn){int ret=comp(num,a.num);if(ret>0)sub(num,a.num);else if(ret<0){bignum_t t ;memcpy(t,num,sizeof(bignum_t));memcpy(num,a.num,sizeof(bignum_t));sub (num,t);sgn=a.sgn ;}else memset(num,0,sizeof(bignum_t)),num[0]=1,sgn=0 ;}else if(!sgn)memcpy(num,a.num,sizeof(bignum_t)),sgn=a.sgn ;return*this ;}inline bignum&operator+=(const int a){if(sgn*a>0)add(num,ABS(a));else if(sgn&&a){int ret=comp(num,ABS(a));if(ret>0)sub(num,ABS(a));else if(ret<0){bignum_t t ;memcpy(t,num,sizeof(bignum_t));memset(num,0,sizeof(bignum_t));num[0]=1 ;add(num,ABS (a));sgn=-sgn ;sub(num,t);}else memset(num,0,sizeof(bignum_t)),num[0]=1,sgn=0 ;}else if(!sgn)sgn=SGN(a),add(num,ABS(a));return*this ;}inline bignum operator+(const bignum&a){bignum ret ;memcpy(ret.num,num,sizeof (bignum_t));ret.sgn=sgn ;ret+=a ;return ret ;}inline bignum operator+(const int a){bignum ret ;memcpy(ret.num,num,sizeof (bignum_t));ret.sgn=sgn ;ret+=a ;return ret ;}inline bignum&operator-=(const bignum&a){if(sgn*a.sgn<0)add(num,a.num);else if(sgn&&a.sgn){int ret=comp(num,a.num);if(ret>0)sub(num,a.num);else if(ret<0){bignum_t t ;memcpy(t,num,sizeof(bignum_t));memcpy(num,a.num,sizeof(bignum_t));sub(num,t);sgn=-sgn ;}else memset(num,0,sizeof(bignum_t)),num[0]=1,sgn=0 ;}else if(!sgn)add (num,a.num),sgn=-a.sgn ;return*this ;}inline bignum&operator-=(const int a){if(sgn*a<0)add(num,ABS(a));else if(sgn&&a){int ret=comp(num,ABS(a));if(ret>0)sub(num,ABS(a));else if(ret<0){bignum_t t ;memcpy(t,num,sizeof(bignum_t));memset(num,0,sizeof(bignum_t));num[0]=1 ;add(num,ABS(a));sub(num,t);sgn=-sgn ;}else memset(num,0,sizeof(bignum_t)),num[0]=1,sgn=0 ;}else if(!sgn)sgn=-SGN(a),add(num,ABS(a));return*this ;}inline bignum operator-(const bignum&a){bignum ret ;memcpy(ret.num,num,sizeof(bignum_t));ret.sgn=sgn ;ret-=a ;return ret ;}inline bignum operator-(const int a){bignum ret ;memcpy(ret.num,num,sizeof(bignum_t));ret.sgn=sgn ;ret-=a ;return ret ;}inline bignum&operator*=(const bignum&a){bignum_t t ;mul(t,num,a.num);memcpy(num,t,sizeof(bignum_t));sgn*=a.sgn ;return*this ;}inline bignum&operator*=(const int a){mul(num,ABS(a));sgn*=SGN(a);return*this ;}inline bignum operator*(const bignum&a){bignum ret ;mul(ret.num,num,a.num);ret.sgn=sgn*a.sgn ;return ret ;}inline bignum operator*(const int a){bignum ret ;memcpy(ret.num,num,sizeof (bignum_t));mul(ret.num,ABS(a));ret.sgn=sgn*SGN(a);return ret ;}inline bignum&operator/=(const bignum&a){bignum_t t ;div(t,num,a.num);memcpy (num,t,sizeof(bignum_t));sgn=(num[0]==1&&!num[1])?0:sgn*a.sgn ;return*this ;}inline bignum&operator/=(const int a){int t ;div(num,ABS(a),t);sgn=(num[0]==1&&!num [1])?0:sgn*SGN(a);return*this ;}inline bignum operator/(const bignum&a){bignum ret ;bignum_t t ;memcpy(t,num,sizeof(bignum_t));div(ret.num,t,a.num);ret.sgn=(ret.num[0]==1&&!ret.num[1])?0:sgn*a.sgn ;return ret ;}inline bignum operator/(const int a){bignum ret ;int t ;memcpy(ret.num,num,sizeof(bignum_t));div(ret.num,ABS(a),t);ret.sgn=(ret.num[0]==1&&!ret.num[1])?0:sgn*SGN(a);return ret ;}inline bignum&operator%=(const bignum&a){bignum_t t ;div(t,num,a.num);if(num[0]==1&&!num[1])sgn=0 ;return*this ;}inline int operator%=(const int a){int t ;div(num,ABS(a),t);memset(num,0,sizeof (bignum_t));num[0]=1 ;add(num,t);return t ;}inline bignum operator%(const bignum&a){bignum ret ;bignum_t t ;memcpy(ret.num,num,sizeof(bignum_t));div(t,ret.num,a.num);ret.sgn=(ret.num[0]==1&&!ret.num [1])?0:sgn ;return ret ;}inline int operator%(const int a){bignum ret ;int t ;memcpy(ret.num,num,sizeof(bignum_t));div(ret.num,ABS(a),t);memset(ret.num,0,sizeof(bignum_t));ret.num[0]=1 ;add(ret.num,t);return t ;}inline bignum&operator++(){*this+=1 ;return*this ;}inline bignum&operator--(){*this-=1 ;return*this ;};inline int operator>(const bignum&a){return sgn>0?(a.sgn>0?comp(num,a.num)>0:1):(sgn<0?(a.sgn<0?comp(num,a.num)<0:0):a.sgn<0);}inline int operator>(const int a){return sgn>0?(a>0?comp(num,a)>0:1):(sgn<0?(a<0?comp(num,-a)<0:0):a<0);}inline int operator>=(const bignum&a){return sgn>0?(a.sgn>0?comp(num,a.num)>=0:1):(sgn<0?(a.sgn<0?comp(num,a.num)<=0:0):a.sgn<=0);}inline int operator>=(const int a){return sgn>0?(a>0?comp(num,a)>=0:1):(sgn<0?(a<0?comp(num,-a)<=0:0):a<=0);}inline int operator<(const bignum&a){return sgn<0?(a.sgn<0?comp(num,a.num)>0:1):(sgn>0?(a.sgn>0?comp(num,a.num)<0:0):a.sgn>0);}inline int operator<(const int a){return sgn<0?(a<0?comp(num,-a)>0:1):(sgn>0?(a>0?comp(num,a)<0:0):a>0);}inline int operator<=(const bignum&a){return sgn<0?(a.sgn<0?comp(num,a.num)>=0:1):(sgn>0?(a.sgn>0?comp(num,a.num)<=0:0):a.sgn>=0);}inline int operator<=(const int a){return sgn<0?(a<0?comp(num,-a)>=0:1):(sgn>0?(a>0?comp(num,a)<=0:0):a>=0);}inline int operator==(const bignum&a){return(sgn==a.sgn)?!comp(num,a.num):0 ;}inline int operator==(const int a){return(sgn*a>=0)?!comp(num,ABS(a)):0 ;}inline int operator!=(const bignum&a){return(sgn==a.sgn)?comp(num,a.num):1 ;}inline int operator!=(const int a){return(sgn*a>=0)?comp(num,ABS(a)):1 ;}inline int operator[](const int a){return digit(num,a);}friend inline istream&operator>>(istream&is,bignum&a){read(a.num,a.sgn,is);return is ;}friend inline ostream&operator<<(ostream&os,const bignum&a){if(a.sgn<0)os<<'-' ;write(a.num,os);return os ;}friend inline bignum sqrt(const bignum&a){bignum ret ;bignum_t t ;memcpy(t,a.num,sizeof(bignum_t));sqrt(ret.num,t);ret.sgn=ret.num[0]!=1||ret.num[1];return ret ;}friend inline bignum sqrt(const bignum&a,bignum&b){bignum ret ;memcpy(b.num,a.num,sizeof(bignum_t));sqrt(ret.num,b.num);ret.sgn=ret.num[0]!=1||ret.num[1];b.sgn=b.num[0]!=1||ret.num[1];return ret ;}inline int length(){return :: length(num);}inline int zeronum(){return :: zeronum(num);}inline bignum C(const int m,const int n){combination(num,m,n);sgn=1 ;return*this ;}inline bignum P(const int m,const int n){permutation(num,m,n);sgn=1 ;return*this ;} };/*int main(){bignum a,b,c;cin>>a>>b;cout<<"加法:"<<a+b<<endl;cout<<"減法:"<<a-b<<endl;cout<<"乘法:"<<a*b<<endl;cout<<"除法:"<<a/b<<endl;c=sqrt(a);cout<<"平方根:"<<c<<endl;cout<<"a的長度:"<<a.length()<<endl;cout<<"a的末尾0個數:"<<a.zeronum()<<endl<<endl;cout<<"組合: 從10個不同元素取3個元素組合的所有可能性為"<<c.C(10,3)<<endl;cout<<"排列: 從10個不同元素取3個元素排列的所有可能性為"<<c.P(10,3)<<endl;return 0 ;}*/int main(int argc, const char * argv[]) {bignum a;while(cin>>a){int cnt=0;int flag=0;bignum b;b = 1;//cout<<b<<endl;//b += 1;while(a != b){a = sqrt(a);//cout<<a<<endl;cnt++;if(cnt>5){flag = 1;break;}}if(flag == 0) printf("%d\n",cnt);else printf("TAT\n");}return 0; } View Code

?

B -?Permutation Bo

這題主要根據期望的可加性，將總的期望轉變為求每一個位置的期望和。

然后推下公式就可以發現在最左邊和最右邊概率為(1/2)

在其他位置的為(1/3). 還得注意n＝＝1 的特殊情況。

// // main.cpp // multi2016.3.B // // Created by New_Life on 16/8/1. // Copyright ? 2016年 chenhuan001. All rights reserved. // #include <iostream> #include <stdio.h> #include <string.h> #include <algorithm> using namespace std;int main(int argc, const char * argv[]) {int n;while(cin>>n){double ans = 0;for(int i=1;i<=n;i++){int tmp;scanf("%d",&tmp);if(i==1 || i==n) ans += (0.5)*tmp;else ans += (1.0/3.0)*tmp;}if(n==1){ans *= 2;}printf("%.6lf\n",ans);}return 0; } View Code

?

C -?Life Winner Bo

不得不說這題就是個代碼量稍微大一點的模擬題。

但是看了大家的題解后發現，怪不得大家做出來的這么多，原來直接用幾種經典的博弈可以直接算出來，仔細一想確實把三種經典博弈全部都用上了。

我用的dp的方法，對于king 和 knight 直接dp搞就行了。

對于車和皇后，則則需要用額外的數組紀錄下。

// // main.cpp // multi2016.3.J // // Created by New_Life on 16/8/1. // Copyright ? 2016年 chenhuan001. All rights reserved. // #include <iostream> #include <string.h> #include <stdio.h> #include <math.h> #include <algorithm> using namespace std;int mark[1002][1002][4]; int save[1002][1002][3];//三個方向//處理king和kinght int dfs(int x,int y,int type) {if(mark[x][y][type] != -1) return mark[x][y][type];int flag = 0;int flag2 = 0;if(type == 0){if(x+1<=1000 && dfs(x+1,y,type)==0) flag=1;if(y+1<=1000 && dfs(x,y+1,type)==0) flag=1;if(x+1<=1000 && y+1<=1000 && dfs(x+1,y+1,type)==0) flag=1;return mark[x][y][type]=flag;}else //只有這里會存在。。。平局 {flag = -1;if(x+2<=1000 && y+1<=1000){if(dfs(x+2,y+1,type) == 0){flag = 1;}else if(dfs(x+2,y+1,type) == 2) flag2= 1;else flag = 0;}if(y+2<=1000 && x+1<=1000){if(dfs(x+1,y+2,type)==0) flag = 1;else if(dfs(x+1,y+2,type)==2) flag2 =1;else flag = 0;}if(flag == 1) return mark[x][y][type] = 1;//能必勝else if(flag2 == 1 || flag == -1) return mark[x][y][type] = 2;else{return mark[x][y][type] = 0;}} }int main() {memset(mark,-1,sizeof(mark));for(int i=0;i<4;i++)mark[1000][1000][i] = 0;for(int i=1;i<=1000;i++)for(int j=1;j<=1000;j++)for(int k=0;k<2;k++){dfs(i,j,k);}//然后來搞車int x[1002],y[1002];memset(x, 0, sizeof(x));memset(y, 0, sizeof(y));for(int i=1000;i>=1;i--)for(int j=1000;j>=1;j--){int flag = 0;if( x[i]==1 || y[j]==1 ) flag = 1;mark[i][j][2] = flag;if(flag == 0){x[i] = 1;y[j] = 1;}}//這里還是錯了for(int i=1000;i>=1;i--)for(int j=1000;j>=1;j--){int flag=0;if(i+1<=1000 && save[i+1][j][0]==1) flag = 1;if(j+1<=1000 && save[i][j+1][1]==1) flag = 1;if(i+1<=1000 && j+1<=1000 && save[i+1][j+1][2]==1) flag = 1;mark[i][j][3] = flag;if(flag == 0){save[i][j][0] = 1;save[i][j][1] = 1;save[i][j][2] = 1;}else{save[i][j][0] = save[i+1][j][0];save[i][j][1] = save[i][j+1][1];save[i][j][2] = save[i+1][j+1][2];}}int T;cin>>T;while(T--){int type,n,m;cin>>type>>n>>m;if(type == 1) type = 0;if(type == 2) type = 2;if(type == 3) type = 1;if(type == 4) type = 3;int flag = mark[1000-n+1][1000-m+1][type];if(flag == 0) printf("G");else if(flag == 1) printf("B");else printf("D");printf("\n");}return 0; }/*101 5 52 5 53 5 54 5 53 2 23 2 33 3 5ans:GGDBDBD */ View Code

?

D -?Gambler Bo

用構造，二分圖，dp什么的鬼想法想了一下午。 后面看題解發現是高斯消元裸題，簡直要撞墻。。。 md taicaile

如果用高斯消元的方法就很好想了。

不過這題還是個線性同余方程組。

具體參見。

?

?G -?Explorer Bo

首先這題確定的最小的移動次數。

稍作分析可以發現，當葉節點個數為偶數時，移動次數必為 (葉節點總數/2)

當葉節點個數為奇數時，移動次數比為（葉節點總數+1）/2

對于偶數的情況，很好分析。 對于每一條樹邊，如果兩側的葉節點個數都為偶數，那么這條樹邊必會走兩次，否則只需要走一次。這樣一趟 dfs 就可以搞定

對于奇數的情況，則需要唯一確定一條不用兩端都為葉節點的路徑。用樹形dp和線段樹都可以解決。 我用的線段樹.

// // main.cpp // hdu5758 // // Created by New_Life on 16/8/3. // Copyright &#169; 2016年 chenhuan001. All rights reserved. // #include <iostream> #include <string.h> #include <stdio.h> #include <algorithm> using namespace std; #define N 100100int n; struct node{int to,next; }edge[2*N];int cnt,pre[N]; int leafid; int savel[N],saver[N]; int ans;void add_edge(int u,int v) {edge[cnt].to = v;edge[cnt].next = pre[u];pre[u] = cnt++; }void dfs(int s,int fa,int first) {int flag = 0;//是否是葉子結點int saveid = leafid;for(int p = pre[s];p!=-1;p=edge[p].next){int v = edge[p].to;if(v == fa) continue;flag ++;dfs(v,s,0);}if(first == 1){if(flag == 1)//is leaf {leafid++;}}elseif(flag == 0){leafid++;}savel[s] = saveid+1;saver[s] = leafid; }void dfs1(int s,int fa) {for(int p=pre[s];p!=-1;p=edge[p].next){int v = edge[p].to;if(v == fa) continue;int x1 = saver[v]-savel[v]+1;//int x2 = leafid - x1;ans += (2-x1%2);dfs1(v,s);} }#define max(a,b) (a>b)?a:b #define min(a,b) (a>b)?b:a #define lson l , m , rt << 1 #define rson m + 1 , r , rt << 1 | 1#define LL long longconst int maxn = N; LL lazy[maxn<<2]; LL sum[maxn<<2];void putup(int rt) {sum[rt] = sum[rt<<1] + sum[rt<<1|1]; } void putdown(int rt,int m) {if (lazy[rt]){lazy[rt<<1] += lazy[rt];lazy[rt<<1|1] += lazy[rt];sum[rt<<1] += lazy[rt] * (m - (m >> 1));sum[rt<<1|1] += lazy[rt] * (m >> 1);lazy[rt] = 0;} } void build(int l,int r,int rt) {lazy[rt] = 0;if (l == r){sum[rt] = 0;//scanf("%I64d",&sum[rt]);return ;}int m = (l + r) >> 1;build(lson);build(rson);putup(rt); } void update(int L,int R,int c,int l,int r,int rt) {if (L <= l && r <= R){lazy[rt] += c;sum[rt] += (LL)c * (r - l + 1);return ;}putdown(rt , r - l + 1);int m = (l + r) >> 1;if (L <= m) update(L , R , c , lson);if (m < R) update(L , R , c , rson);putup(rt); } LL query(int L,int R,int l,int r,int rt) {if (L <= l && r <= R){return sum[rt];}putdown(rt , r - l + 1);int m = (l + r) >> 1;LL ret = 0;if (L <= m) ret += query(L , R , lson);if (m < R) ret += query(L , R , rson);return ret; } /*int main(){int n , m;int a , b , c;char str[5];scanf("%d%d",&n,&m);build(1 , n , 1);while (m--){scanf("%s",str);if (str[0] == 'Q'){scanf("%d%d",&a,&b);printf("%I64d\n",query(a , b , 1 , n , 1));}else if(str[0]=='C'){scanf("%d%d%d",&a,&b,&c);update(a , b , c , 1 , n , 1);}}return 0;}*/void dfs2(int s,int fa) {for(int p=pre[s];p!=-1;p=edge[p].next){int v = edge[p].to;if(v == fa) continue;int x1 = saver[v]-savel[v]+1;int x2 = leafid - x1;int tmp1=0,tmp2=0;x1--;if(x1!=0) tmp1 = (2-x1%2);else tmp1 = 1;x1 ++;x2--;if(x2!=0) tmp2 = (2-x2%2);else tmp2 = 1;x2++;update(1,leafid,tmp2,1,leafid,1);update(savel[v],saver[v],tmp1-tmp2,1,leafid,1);dfs2(v,s);} }int main(int argc, const char * argv[]) {int T;cin>>T;while(T--){cin>>n;cnt = 0;memset(pre,-1,sizeof(pre));for(int i=1;i<n;i++){int x,y;scanf("%d%d",&x,&y);add_edge(x,y);add_edge(y,x);}if(n == 1){printf("0\n");continue;}leafid = 0;dfs(1, -1, 1);//第一遍建立線段區間，且得出有多少葉節點. ans = 0;if(leafid%2 == 0)//葉節點為偶數 {//一趟dfs搞定dfs1(1, -1);}else{build(1, leafid, 1);dfs2(1,-1);ans = 1000000000;for(int i=1;i<=leafid;i++){ans = min(ans,(int)query(i, i, 1, leafid, 1));}}printf("%d\n",ans);}return 0; } /*1041 22 33 451 55 22 32 491 22 32 43 53 64 77 87 9*/ View Code

?

J -?Rower Bo

應該是高中的物理競賽會考的東西。

這種巧妙的公式，沒做過的人好難。 聽說很多人都是直接猜的公式，這場比賽3題都可以找規律。 orz

?

K -?Teacher Bo

因為n*m的圖的曼哈頓距離最多有2*(n+m)種，所以當n>500 時必有相同點，n<500直接暴力。

// // main.cpp // multi2016.3.B // // Created by New_Life on 16/8/1. // Copyright ? 2016年 chenhuan001. All rights reserved. // #include <iostream> #include <stdio.h> #include <string.h> #include <math.h> #include <algorithm> using namespace std;struct point {int x,y; }g[100100];int mark[200200];//僅適合純數字輸入 int Scan() //輸入外掛 {int res=0,ch,flag=0;if((ch=getchar())=='-')flag=1;else if(ch>='0'&&ch<='9')res=ch-'0';while((ch=getchar())>='0'&&ch<='9')res=res*10+ch-'0';return flag?-res:res; } void Out(int a) //輸出外掛 {if(a>9)Out(a/10);putchar(a%10+'0'); }int main(int argc, const char * argv[]) {int T;cin>>T;while(T--){int n,m;scanf("%d%d",&n,&m);for(int i=0;i<n;i++){g[i].x = Scan();g[i].y = Scan();}if(n>500) printf("YES\n");else{memset(mark,0,sizeof(mark));//for(int i=0;i<2*m;i++)int flag=0;for(int i=0;i<n;i++){if(flag) break;for(int j=i+1;j<n;j++){int tmp = abs(g[i].x-g[j].x)+abs(g[i].y-g[j].y);mark[ tmp ]++;if(mark[tmp] >= 2){flag = 1;break;}}}if(flag) printf("YES\n");else printf("NO\n");}}return 0; } View Code

?

總結

以上是生活随笔為你收集整理的2016 Multi-University Training Contest 3的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。