Noip2016

<这篇是以前的，不开新的了，借版面来换了个标题>

高二了

开学一周，每天被文化课作业碾压。。。

但是仍然阻挡不了想刷题的心情。。。

对付noip2016的几块：（有点少，以后补）

高精度（c++模板仍未生成）

[bzoj2656]

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int cas;
struct data{
    ],l;
    data(){
        v[]=;l=;
    }
    void clear(){
        memset(v,,;
    }
}n,p,q;
];
void read(data &x){
    data a;scanf();a.l=strlen(ch+);
    ;i<=a.l;i++)a.v[i]=ch[a.l-i+]-';
    ;i<=a.l+;i++)a.v[i]=;//ºóÃæµÄÊý×é×îºÃ¿ÉÒÔÇå¿ÕÒ»ÏÂ
    x.l=(a.l-)/+;
    ;i<=x.l;i++)
    {
       x.v[i]=;
       ;k<=;k++)
       {
          x.v[i]=x.v[i]*+a.v[(i-)*+(-k)];
       }
    }
}
data operator+(data a,data b){
    if(a.l<b.l)swap(a,b);int t=b.l;
    ;i<=t;i++)a.v[i]+=b.v[i];
    ;i<=t;i++)
      ){
          if(i==a.l){
             a.v[i+]=;a.l++;
          }
          a.v[i+]+=a.v[i]/;t=max(t,i+);
          a.v[i]=a.v[i]%;
      }
    return a;
}
data operator+(data a,int p)
{
    data b;b.v[]=p;b.l=;return a+b;
}
data operator/(data a,int p){
    data b;;
    ;i--)
    {
       all=all*+a.v[i];b.v[i]=all/p;all=all-b.v[i]*p;
    }
    ;i<=a.l;i++)a.v[i]=b.v[i];
    &&a.l>)a.l--;return a;
}
void calc(data x)
{
    &&x.v[]==){
        p=x;q.clear();return;
    }
    calc((x+)/);]&)p=p+q;else q=p+q;
}
void print(data a){
    printf("%d",a.v[a.l]);
    ;i>;i--)printf("%04d",a.v[i]);printf("\n");
}
int main()
{
    scanf("%d",&cas);
    while(cas--){
        read(n);calc(n);print(p);
    }
} 

高精度重载+/

图论相关：欧拉路 双联通分量（以及构造）...点分治？

#include<cstdio>
#include<algorithm>
#include<cstring>
#define N 500000
using namespace std;
int u,v,edgenum,nl,nr,m,ans;
int vet[N],head[N],next[N],flag[N],match[N];
void add(int u,int v)
{
  edgenum++;vet[edgenum]=v;next[edgenum]=head[u];
  head[u]=edgenum;
}
bool dfs(int u)
{
  flag[u]=;
  int e=head[u];
  )
  {
    int v=vet[e];
    ||flag[match[v]]==&&dfs(match[v]))
    {
       match[v]=u;match[u]=v;
       return true;
    }
    e=next[e];
  }
  return false;
}
int main()
{
  scanf("%d%d%d",&nl,&nr,&m);
  ;i<=m;i++)
  {
    scanf("%d%d",&u,&v);
    v+=nl;add(u,v);
  }
  ;i<=nl;i++)

    {
       memset(flag,,sizeof(flag));
       if(dfs(i))ans++;
    }
  printf("%d\n",ans);
  ;i<=nl;i++))match[i]-=nl;
  ;i<=nl;i++)printf("%d ",match[i]);
}

二分图匹配

可行遍性问题：

hamilton回路问题（它是一个np问题，目前没有多项式算法，连判定都需要暴搜）：每个点恰好走一次回到出发点。充分条件： 无向连通图中任意2点度数之和大于等于顶点数，则必定存在哈密顿回路。

无向图的欧拉回路：所有顶点为偶点

有向图的欧拉回路：所有顶点入度等于出度

一个图存在欧拉路当且仅当存在两个奇点（度数为奇数的点）或不存在奇点。

暴搜（去年没好好练吃亏了）

各式各样的dp（也包括找规律）

数论相关：联赛层面应该差不多了，以防万一先留着吧。

高斯消元，只剩下背背背

奇奇怪怪的一些写法：归排 差分 2-sat...想不出来的时候试试二分大法

网络流 可能性不大 准备了总是不会错的

kmp:有个模板是字符串重复那道题，代码在暑期集训day3那里。。。下面贴一段很详细的题解：

假设字符串为S，长度为N，子串T重复K次后得到串S，那么T的长度一定为L = N/K(要整除)，则T = S[1...L]，将S拆分成K份，每份长度为L，则有

S[1...L] = S[L+1...2L] = S[2L+1...3L] = ... = S[(K-1)L+1...KL]

由于要保证K最大，势必L要取最小，所以根据Next函数的定义，有Next[KL] = (K-1)L;

即Next[N] = N - L，所以L = N - Next[N];

但是得出的长度L还要保证能被N整除，所以如果不能整除说明L = N，即K = 1；而如果能整除，那么K = N / (N - Next[N]);

因而我们只要对字符串S做一趟KMP，对其求Next数组，剩下的就是上述结论

时间复杂度O(n）

偶然在一篇blog里看到的，借来共勉：

尽人事，听天命。相信自己可以拿出最好的状态。还是一样，考场思路开阔，胆大心细，勇于乱搞，rp++.

愿这是一场属于所有OIer的盛宴.

day1：先简单记录一下，希望明天加油！

T1：模拟，希望不会写挂，rp++

T2：这样看来好像还真不能说出题人超纲了，好多神犇都用大力数据结构维护了，我稍稍想到了一些意思但是没有什么卵用。希望暴力不要写挂。对拍助我。

T3：确实蛮后悔的，不该跑最小生成树，但是现在说什么都晚了。希望数据水一点。。。能给我一些分数吧。。。另外感觉如果这题在T2应该会有更多人AC吧。复述一下nbc的状态——f[i][j]=min{f[i-1][j]+dis[c[i-1]][c[i]]},f[i][j]=min{g[i-1][j]+dis[d[i-1]][c[i]]*p[i-1]+(1-p[i-1])*dis[c[i-1]][c[i]]}……claris说知道期望可加性就好了。。。唉。。。