NOIP2015 Revenge

辣鸡出题人，毁我比赛，颓我精神，耗我钱财，废我青春。

去年成绩惨不忍睹就不说了...好像是100+80+0+100+50+60。

大概列一下算法。

幻方：模拟

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <memory.h>
using namespace std;
int n,mat[50][50],px,py;
#define FO(s) {freopen(#s".in","r",stdin);freopen(#s".out","w",stdout);}
int main()
{
cin>>n;
memset(mat,0,sizeof(mat));
px=1; py=(n+1)/2;
mat[px][py]=1;
//By zzq
for(int i=2;i<=n*n;i++)
{
if(px==1&&py!=n)
{
px=n; ++py; mat[px][py]=i;
}
else if(py==n&&px!=1)
{
--px; py=1; mat[px][py]=i;
}
else if(px==1&&py==n)
{
++px; mat[px][py]=i;
}
else
{
if(!mat[px-1][py+1])
{
--px; ++py; mat[px][py]=i;
}
else
{
++px; mat[px][py]=i;
}
}
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
printf("%d",mat[i][j]);
if(j!=n) putchar(' ');
}
putchar('\n');
}
}

信息传递：tarjan/鬼畜dfs（各种写丑写成平方的都有80）

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <stack>
using namespace std;
//tarjan!
int n,t[233333],low[233333],dfn[233333],cur=0,siz=2000000000;
bool instack[233333],ved[233333];
int ss[233333],sn=0;
void tarjan(int s)
{
dfn[s]=low[s]=++cur; ss[++sn]=s;
instack[s]=1; ved[s]=1;
if(instack[t[s]]) low[s]=min(low[s],dfn[t[s]]);
else if(!dfn[t[s]])
{
tarjan(t[s]);
low[s]=min(low[s],low[t[s]]);
}
if(dfn[s]==low[s])
{
int cs=0;
while(1)
{
++cs;
instack[ss[sn]]=0;
int t=ss[sn--];
if(t==s) break;
}
if(cs!=1) siz=min(siz,cs);
}
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",t+i);
for(int i=1;i<=n;i++) if(!ved[i]) tarjan(i);
printf("%d\n",siz);
}

斗地主：欧洲人

跳石头：二分

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <set>
using namespace std;
int D[233333],L,n,m;
bool ok(int x)
{
int cp=0,del=0;
for(int i=1;i<=n;i++)
{
if(D[i]-cp<x)
{
del++; continue;
}
cp=D[i];
}
if(L-cp<x) del++;
return del<=m;
}
int main()
{
scanf("%d%d%d",&L,&n,&m);
for(int i=1;i<=n;i++) scanf("%d",D+i);
int l=0,r=L;
while(l<r)
{
int mid=(l+r+1)>>1;
if(ok(mid)) l=mid; else r=mid-1;
}
printf("%d\n",l);
}

子串：dp

dp[k][i][j]表示最优情况下使用k个子串，用了A串第i个构成B串前j个的方案数。

dp[k][i][j]=qzh[k-1][i-1][j-1] （注意这个前缀和只有第二维）

如果a[i-1]=b[j-1]，那么dp[k][i][j]+=dp[k][i-1][j-1]

考场上不知道是哪里抽风了顺手多记了一维= =

似乎要随便滚动一下

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
using namespace std;
int MOD=1000000007;
int n,m,w,qzh[2][1005][205],dp[2][1005][205];
char a[233333],b[233333];
int main()
{
scanf("%d%d%d%s%s",&n,&m,&w,a+1,b+1);
dp[0][0][0]=1;
for(int i=0;i<=n;i++) qzh[0][i][0]=1;
for(int k=1;k<=w;k++)
{
for(int i=0;i<=n;i++)
{
for(int j=0;j<=m;j++) qzh[k&1][i][j]=dp[k&1][i][j]=0;
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(a[i]!=b[j]) continue;
int&cp=dp[k&1][i][j];
cp=qzh[(k-1)&1][i-1][j-1];
if(a[i-1]==b[j-1]) cp=(cp+dp[k&1][i-1][j-1])%MOD;
}
for(int j=0;j<=m;j++) qzh[k&1][i][j]=(qzh[k&1][i-1][j]+dp[k&1][i][j])%MOD;
}
}
printf("%d\n",qzh[w&1]
[m]);
}

运输计划：乱搞

考虑二分答案，对于每一个超过限制的，我们必然要在其上选择一条边删去，且该条边权>=实际值-二分值。

那么我们暴力路径+1，再暴力for一下每一条边即可。

什么？怎么+1？前缀和！

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define SZ 666666
#define P 20
typedef pair<int,int> pii;
int n,m,fa[SZ],fv[SZ],vc[SZ],dep[SZ],fst[SZ],vb[SZ],nxt[SZ],M=0,qa[SZ],qb[SZ],ql[SZ],ea[SZ],eb[SZ],ec[SZ],reald[SZ];
void ad_de(int a,int b,int c) {++M; nxt[M]=fst[a]; fst[a]=M; vb[M]=b; vc[M]=c;}
int cc=0,app[SZ],bs[SZ],up[SZ][P],cx[SZ],c2=0;
pii pp[SZ],minn[SZ][P];
void dfs(int x)
{
++cc; app[x]=cc; pp[cc]=pii(dep[x],x); cx[++c2]=x;
for(int e=fst[x];e;e=nxt[e])
{
int b=vb[e]; if(b==fa[x]) continue;
fa[b]=x; fv[b]=fv[x]+vc[e]; dep[b]=dep[x]+1;
dfs(b); pp[++cc]=pii(dep[x],x);
}
}
void build()
{
for(int i=1;i<=n;i++) up[i][0]=fa[i];
for(int p=1;p<P;p++)
{
for(int i=1;i<=n;i++)
{
if(up[i][p-1]) up[i][p]=up[up[i][p-1]][p-1];
}
}
for(int i=1;i<=cc;i++) minn[i][0]=pp[i];
for(int i=1;i<=cc;i++)
{
int g=0;
while((1<<g)<=i) ++g;
bs[i]=g-1;
}
for(int p=1;p<P;p++)
{
for(int i=1;i<=cc;i++)
{
if(i+(1<<p)-1>cc) break;
minn[i][p]=min(minn[i][p-1],minn[i+(1<<(p-1))][p-1]);
}
}
}
int jump(int x,int d)
{
for(int i=P-1;i>=0;i--)
{
if(up[x][i]&&dep[up[x][i]]>=d) x=up[x][i];
}
return x;
}
int lca(int a,int b)
{
a=app[a]; b=app[b];
if(a>b) swap(a,b);
int l2=bs[b-a+1];
return min(minn[a][l2],minn[b-(1<<l2)+1][l2]).second;
}
int dis(int a,int b,int l)
{
return fv[a]+fv[b]-fv[l]*2;
}
int pv[SZ],cv[SZ];
bool ok(int x)
{
int cnt=0,need=0;
for(int i=1;i<=n;i++) pv[i]=0;
for(int i=1;i<=m;i++)
{
if(reald[i]<=x) continue;
++cnt; pv[ql[i]]-=2; pv[qa[i]]++; pv[qb[i]]++;
need=max(need,reald[i]-x);
}
if(!cnt) return 1;
for(int i=n;i>=1;i--) pv[fa[cx[i]]]+=pv[cx[i]];
for(int i=1;i<n;i++)
{
if(ec[i]<need) continue;
if(pv[eb[i]]!=cnt) continue;
return 1;
}
return 0;
}
#define BUFSIZE 300000
namespace fib {char b[BUFSIZE]={},*f=b;}
#define gc ((*fib::f)?(*(fib::f++)):(fgets(fib::b,sizeof(fib::b),stdin)?(fib::f=fib::b,*(fib::f++)):-1))
int g_i()
{
int tmp=0; bool fu=0; char s;
while(s=gc,s!='-'&&(s<'0'||s>'9')) ;
if(s=='-') fu=1; else tmp=s-'0';
while(s=gc,s>='0'&&s<='9') tmp=tmp*10+s-'0';
if(fu) return -tmp; else return tmp;
}
#define gi g_i()
namespace fob {char b[BUFSIZE]={},*f=b,*g=b+BUFSIZE-2;}
#define pob (fwrite(fob::b,sizeof(char),fob::f-fob::b,stdout),fob::f=fob::b,0)
#define pc(x) (*(fob::f++)=(x),(fob::f==fob::g)?pob:0)
struct foce {~foce() {pob; fflush(stdout);}} _foce;
namespace ib {char b[100];}
inline void pint(int x)
{
if(x==0) {pc(48); return;}
if(x<0) {pc('-'); x=-x;}
char *s=ib::b;
while(x) *(++s)=x%10, x/=10;
while(s!=ib::b) pc((*(s--))+48);
}
int main()
{
n=gi,m=gi;
for(int i=1;i<n;i++)
{
int a=gi, b=gi, c=gi;
ea[i]=a; eb[i]=b; ec[i]=c;
ad_de(a,b,c); ad_de(b,a,c);
}
for(int i=1;i<=m;i++) qa[i]=gi, qb[i]=gi;
dfs(1); build();
int l=0,r=0;
for(int i=1;i<=m;i++) ql[i]=lca(qa[i],qb[i]), reald[i]=dis(qa[i],qb[i],ql[i]), r=max(r,reald[i]);
for(int i=1;i<n;i++) if(dep[ea[i]]>dep[eb[i]]) swap(ea[i],eb[i]);
while(l<r)
{
int mid=l+r>>1;
if(ok(mid)) r=mid; else l=mid+1;
}
printf("%d\n",l);
}

然后我们只要兹磁O(1) lca就行了。可能需要卡常~（我卡了三发才过）