- 上午
- 模拟考试
- Prob.1（AC）一道简单的博弈题，找到必胜态，反推普遍情况是否可以达到必胜态即可。
- Prob.2（AC）做到原题了呢。入门OJ 2092: [Noip模拟题]舞会
- Prob.3（WA了3个点）一道高精度的dp题，为减少状态，我把数据按某一权值分为了两类，对于一类反着dp，状态很少，只用到高精度加法。对于另一类（30分），要用到高精度乘法，然后进位时少打了一个+号，然后就没有然后了……
- 然后整理了三个模板，贴上！
//maybe have no mistakes,^_^ //1.Big_Int //only support + , * , = , < , <= , > , >= between Big_Int and Big_Int as well as Big_Int and int #define MAXN 1000 struct Big_Int{ int len,a[30]; Big_Int(){len=1; memset(a,0,sizeof(a));} void operator =(int val){ len=0; do{ a[++len]=val%MAXN; val/=MAXN; }while(val); } Big_Int operator +(const Big_Int &rtm) const{ Big_Int now; int l=max(len,rtm.len); for(int i=1;i<=l;i++){ now.a[i]+=a[i]+rtm.a[i]; now.a[i+1]=now.a[i]/MAXN; now.a[i]%=MAXN; } if(now.a[l+1]) l++; now.len=l; return now; } Big_Int operator +(const int &val) const{ Big_Int now,rtm; rtm=val; now=(*this)+rtm; return now; } Big_Int operator *(const Big_Int &rtm) const{ Big_Int now; int l=len+rtm.len; for(int i=1;i<=len;i++) for(int j=1;j<=rtm.len;j++){ now.a[i+j-1]+=a[i]*rtm.a[j]; now.a[i+j]+=now.a[i+j-1]/MAXN; now.a[i+j-1]%=MAXN; } while(!now.a[l]) l--; now.len=l; return now; } Big_Int operator *(const int &val) const { Big_Int now,rtm; rtm=val; now=(*this)*rtm; return now; } bool operator ==(const Big_Int & rtm) const{ if(len!=rtm.len) return 0; for(int i=len;i>=1;i--) if(a[i]!=rtm.a[i]) return 0; return 1; } bool operator <(const Big_Int &rtm) const{ if(len!=rtm.len) return len<rtm.len; for(int i=len;i>=1;i--) if(a[i]!=rtm.a[i]) return a[i]<rtm.a[i]; return 0; } bool operator >(const Big_Int &rtm) const{ return rtm<(*this); } bool operator <=(const Big_Int &rtm) const{ return (*this)<rtm||(*this)==rtm; } bool operator >=(const Big_Int &rtm) const{ return (*this)>rtm||(*this)==rtm; } bool operator ==(const int &val) const{ Big_Int rtm; rtm=val; return (*this)==rtm; } bool operator <(const int &val) const{ Big_Int rtm; rtm=val; return (*this)<rtm; } bool operator >(const int &val) const{ Big_Int rtm; rtm=val; return (*this)<rtm; } bool operator <=(const int &val) const{ Big_Int rtm; rtm=val; return (*this)<=rtm; } bool operator >=(const int &val) const{ Big_Int rtm; rtm=val; return (*this)>=rtm; } void print(){ printf("%d",a[len]); for(int i=len-1;i>=1;i--){ printf("%03d",a[i]); } } }; //2.Fraction //only support + , * and reduction for fraction struct Fraction{ ll Numerator,Denominator; ll Gcd(ll a,ll b){ while(a^=b^=a^=b%=a); return b; } bool Zero(){ return Numerator==0; } void Clear(){ Numerator=Denominator=0; } void Reduction(){ ll g=Gcd(Numerator,Denominator); Numerator/=g; Denominator/=g; } Fraction operator +(Fraction &rtm){ Fraction New; if(!Numerator&&!Denominator) New.Numerator=rtm.Numerator,New.Denominator=rtm.Denominator; else{ ll g=Gcd(Denominator,rtm.Denominator); ll a=Denominator/g,b=rtm.Denominator/g; New.Denominator=a*b*g; New.Numerator=Numerator*b+rtm.Numerator*a; } New.Reduction(); return New; } Fraction operator *(Fraction &rtm){ Fraction New; New.Denominator=Denominator*rtm.Denominator; New.Numerator=Numerator*rtm.Numerator; New.Reduction(); return New; } void Read(){ cin>>Numerator; Denominator=1; } void Write(){ cout<<Numerator; if(Denominator!=1&&Denominator!=0) cout<<"/"<<Denominator<<endl; } } //3.Maxtrix //noly supprt * , ^(quick pow) for Maxtrix struct Maxtrix{ int r,c; int val[75][75]; Maxtrix(){ r=c=0; memset(val,0,sizeof(val)); } void identity(int l){ r=c=l; for(int i=1;i<=l;i++) val[i][i]=1; } Maxtrix operator * (Maxtrix const b){ // c==b.r Maxtrix now; now.r=r; now.c=b.c; for(int i=1;i<=r;i++) for(int j=1;j<=b.c;j++) for(int k=1;k<=c;k++) now.val[i][j]=(now.val[i][j]+val[i][k]*b.val[k][j])%mod; return now; } Maxtrix operator ^ (int b){ Maxtrix base,now; base=*this; now.identity(this->r); while(b){ if(b&1) now=now*base; base=base*base; b>>=1; } return now; } } //____________________________________________________________________by *ZJ
- 模拟考试
- 下午
- 入门OJ 2069: [Noip模拟题]贪吃的CWT
对拍了一下，改出来了呢。
知道了错误在哪儿的我想打人……
这是我原来错误的dis函数：
（真不知道是什么鬼样例数据，dis没开根都过了）
代码：
#include<cmath> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> using namespace std; struct edge{ int u,v; double val; bool operator<(const edge &rtm) const{ return val<rtm.val; } }e[1005*1005]; struct hamburger{ int x,y,val; }p[1005]; int fa[1005],bh[1005]; int n,ent; double tot,ans; double dis(int i,int j){ return sqrt((p[i].x-p[j].x)*(p[i].x-p[j].x)+(p[i].y-p[j].y)*(p[i].y-p[j].y)); } int find(int x){ return fa[x]==x?x:fa[x]=find(fa[x]); } int main(){ scanf("%d",&n); for(int i=1,a,b,c;i<=n;i++){ scanf("%d%d%d",&a,&b,&c); p[i]=(hamburger){a,b,c}; } for(int i=1;i<n;i++) for(int j=i+1;j<=n;j++) e[++ent]=(edge){i,j,dis(i,j)}; sort(e+1,e+ent+1); for(int i=1;i<=n;i++) fa[i]=i; int fu,fv,u,v,need=n-1; for(int i=1;i<=ent&&need;i++){ u=e[i].u; v=e[i].v; fu=find(u); fv=find(v); if(fu==fv) continue; tot+=e[i].val; fa[fv]=fu; need--; } for(int i=1;i<=n;i++) fa[i]=i; for(int i=1;i<=n;i++) bh[i]=p[i].val; need=n-1; for(int i=1;i<=ent&&need;i++){ u=e[i].u; v=e[i].v; fu=find(u); fv=find(v); if(fu==fv) continue; ans=max(ans,1.0*(bh[fu]+bh[fv])/(tot-e[i].val)); bh[fu]=max(bh[fu],bh[fv]); fa[fv]=fu; need--; } printf("%.2lf",ans); return 0; }
- 入门OJ 2069: [Noip模拟题]贪吃的CWT
- 入门OJ 2073: [Noip模拟题]古巴比伦
一个皮皮题。
我直接跑了个n2就过了，应该是数据水了……，
但erge说if语句和++操作常数很小，跑108可以过！
他的做法是：因为数据完全随机，所以当两个串的重叠部分长度小于某个值时，就可以直接跳出，
至于小于哪个值，他说：“就自己估摸吧，我定的不能小于90%的长度。”
……
代码：
#include<cstdio> #include<cstring> #include<iostream> using namespace std; char s1[10005],s2[10005]; int ls1,ls2,cnt,ans; void doit(char *A,int LA,char *B,int LB){ for(int i=0;i<LA;i++){ cnt=0; for(int j=0;j<LB&&i+j<LA;j++) if(A[i+j]==B[j]) cnt++; ans=max(ans,cnt); } } int main(){ scanf("%d",&ls1); scanf(" %s",s1); scanf("%d",&ls2); scanf(" %s",s2); doit(s1,ls1,s2,ls2); doit(s2,ls2,s1,ls1); printf("%d",ans); return 0; }
- 晚上
- 回家了，lalala
- End
-
没有end，只有放假！
-