T1
好水啊~~~
#include <cstdio> #include <cstring> #include <cstdlib> #include <iostream> #include <algorithm> #define ll long long #define mem(a,b) memset(a,b,sizeof(a)) #define rint register int using namespace std; inline void read(int &x) { x=0; char q=getchar(); while(q<'0'||q>'9') q=getchar(); while(q>='0'&&q<='9') x=x*10+q-'0',q=getchar(); } const int N=100006; int first[N],nt[N*2],ver[N*2],e; void addbian(int u,int v) { ver[e]=v; nt[e]=first[u]; first[u]=e++; } int n,an; int fa[N],xi[N],id[N]; void dfs(int x) { xi[x]=1; int i; for(i=first[x];i!=-1;i=nt[i]) { if(ver[i]==fa[x]) continue; fa[ver[i]]=x; id[ver[i]]=(i>>1)+1; dfs(ver[i]); xi[x]+=xi[ver[i]]; } } int main(){ //freopen("T1.in","r",stdin); //freopen("div.in","r",stdin); //freopen("div.out","w",stdout); rint i; mem(first,-1); read(n); int tin1,tin2; for(i=1;i<n;++i) { read(tin1); read(tin2); addbian(tin1,tin2); addbian(tin2,tin1); } dfs(1); an=-1; for(i=1;i<=n;++i) if(id[i]&&xi[i]*2==n) an=id[i]; cout<<an; fclose(stdin); fclose(stdout); }
T2
题解释哈夫曼树，然而我并不是这样想的...
可以考虑在01trie上行走
那么它一定是在每一个深度都用完了，在到下一层
可以脑补一下...
#include <cstdio> #include <cstring> #include <cstdlib> #include <iostream> #include <algorithm> #define ll long long #define mem(a,b) memset(a,b,sizeof(a)) #define rint register int using namespace std; inline void read(int &x) { x=0; char q=getchar(); while(q<'0'||q>'9') q=getchar(); while(q>='0'&&q<='9') x=x*10+q-'0',q=getchar(); } const int N=10000006; const int mod=1000000007; ll qpow(ll a,int ci) { ll ans=1; while(ci) { if(ci&1) ans=ans*a%mod; a=a*a%mod; ci>>=1; } return ans; } int n; ll jie[N],jieni[N]; void chu() { int i; jie[0]=1; for(i=1;i<N;++i) jie[i]=jie[i-1]*i%mod; jieni[N-1]=qpow(jie[N-1],mod-2)%mod; for(i=N-2;i>=1;--i) jieni[i]=jieni[i+1]*(i+1)%mod; jieni[0]=1; } ll C(int n,int m) { if(n<0||m<0||n<m) return 0; return jie[n]*jieni[m]%mod*jieni[n-m]%mod; } int main(){ //freopen("str.in","r",stdin); //freopen("str.out","w",stdout); chu(); read(n); int k=0; while((1<<k)<=n) ++k; --k; printf("%lld",C((1<<k),(n-(1<<k)))%mod*jie[n]%mod ); fclose(stdin); fclose(stdout); }
T3
(1)
你可以大暴搜
#include <cstdio> #include <cstring> #include <cstdlib> #include <iostream> #include <algorithm> #define ll long long #define mem(a,b) memset(a,b,sizeof(a)) #define rint register int using namespace std; inline void read(int &x) { x=0; char q=getchar(); while(q<'0'||q>'9') q=getchar(); while(q>='0'&&q<='9') x=x*10+q-'0',q=getchar(); } const int N=5006; const int mod=1000000007; int n,K; int ans; void dfs(int order,int now,int mn,int mx) { if(mx-mn>K) return ; if(order==n+1) { ++ans; return ; } dfs(order+1,now+1,mn,max(mx,now+1)); dfs(order+1,now+2,mn,max(mx,now+2)); dfs(order+1,now-1,min(mn,now-1),mx); } int main(){ //freopen("T3.in","r",stdin); read(n); read(K); dfs(1,0,0,0); cout<<ans; }
(2)
$O(nK^3)$的dp
$f_{i,j,k,v}$ 表示当前走到了第i个位置，最大值是j，最小值是k，当前的权值是v的方案数
#include <cstdio> #include <cstring> #include <cstdlib> #include <iostream> #include <algorithm> #define ll long long #define mem(a,b) memset(a,b,sizeof(a)) #define rint register int using namespace std; inline void read(int &x) { x=0; char q=getchar(); while(q<'0'||q>'9') q=getchar(); while(q>='0'&&q<='9') x=x*10+q-'0',q=getchar(); } const int N=5006; const int mod=1000000007; const int BAO=78; int n,K; ll f[2][BAO*3][BAO][BAO*2]; ll dp() { rint i,j,k,p; int now=1,mn=BAO-n,mx=n+n+BAO,nx,tmx,tt; //printf("mn=%d mx=%d ",mn,mx); mem(f,0); f[0][BAO][BAO][0]=1; for(i=0;i<n;++i) { now^=1; nx=now^1; mem(f[nx],0); for(j=mn;j<=mx;++j) for(k=BAO-i;k<=j;++k) for(p=max(j,BAO);p<=mx;++p) if(f[now][j][k][p-BAO]) { //printf("i=%d ",i); tmx=p-BAO; tt=f[now][j][k][tmx]; f[nx][j+1][k][max(tmx,j+1-BAO)]+=tt; f[nx][j+1][k][max(tmx,j+1-BAO)]%=mod; f[nx][j+2][k][max(tmx,j+2-BAO)]+=tt; f[nx][j+2][k][max(tmx,j+2-BAO)]%=mod; f[nx][j-1][min(j-1,k)][tmx]+=tt; f[nx][j-1][min(j-1,k)][tmx]%=mod; } } ll ans=0; for(i=mn;i<=mx;++i) for(j=mn;j<=BAO;++j) { tt=min(j+K,mx); for(k=j;k<=tt;++k) ans=(ans+f[nx][i][j][k-BAO])%mod; } return ans; } int main(){ freopen("T3.in","r",stdin); read(n); read(K); printf("%lld",dp()); }
(3)
$O(nK^2)$ 的dp
$f_{i,j,k}$ 表示当前走到了第i个位置，最大值和最小值的差值是j，当前值比最小值大k的方案数
特别的，差值j是当前走的路线到过了最大值和最小值才行，所以不会重复
#include <cstdio> #include <cstring> #include <cstdlib> #include <iostream> #include <algorithm> #define ll long long #define mem(a,b) memset(a,b,sizeof(a)) #define rint register int using namespace std; inline void read(int &x) { x=0; char q=getchar(); while(q<'0'||q>'9') q=getchar(); while(q>='0'&&q<='9') x=x*10+q-'0',q=getchar(); } const int mod=1000000007; int n,K; ll f[2][306][306]; ll dp() { int i,j,k,tt,now=1,nx,nj,nk; ll ans=0; f[0][0][0]=1; for(i=0;i<n;++i) { now^=1; nx=now^1; mem(f[nx],0); for(j=0;j<=2*i&&j<=K;++j) for(k=0;k<=j;++k) if(f[now][j][k]) { nk=k+1; nj=j; if(nj<nk) nj=nk; if(nj<=K) f[nx][nj][nk]=(f[nx][nj][nk]+f[now][j][k])%mod; nk=k+2; nj=j; if(nj<nk) nj=nk; if(nj<=K) f[nx][nj][nk]=(f[nx][nj][nk]+f[now][j][k])%mod; nk=k-1; nj=j; if(nk<0) nj=nj-nk,nk=0; if(nj<=K) f[nx][nj][nk]=(f[nx][nj][nk]+f[now][j][k])%mod; } } for(i=0;i<=K;++i) for(j=0;j<=i;++j) ans=(ans+f[nx][i][j])%mod; return ans%mod; } int main(){ read(n); read(K); printf("%lld",dp()); }
(4)
$O(nK)$ 的dp
我们考虑在一个n*K大小的网格里从最左边往右走
枚举起点在最左边的位置，就相当于枚举最小值
$f_{i,j,k}$ 表示当前走到了第i个位置，当前的权值是j，到没到过0(也就是最小值)的方案数
这样开第三维也可以保证不重复
#include <cstdio> #include <cstring> #include <cstdlib> #include <iostream> #include <algorithm> #define ll long long #define mem(a,b) memset(a,b,sizeof(a)) #define rint register int using namespace std; inline void read(int &x) { x=0; char q=getchar(); while(q<'0'||q>'9') q=getchar(); while(q>='0'&&q<='9') x=x*10+q-'0',q=getchar(); } const int mod=1000000007; int n,K; ll f[2][5006][2]; ll dp() { int i,j,k,tt,now=1,nx; ll ans=0; f[0][0][1]=1; for(i=1;i<=K;++i) f[0][i][0]=1; for(i=0;i<n;++i) { now^=1; nx=now^1; mem(f[nx],0); for(j=0;j<=K;++j) { tt=j+1; if(tt<=K) f[nx][tt][0]=(f[nx][tt][0]+f[now][j][0])%mod,f[nx][tt][1]=(f[nx][tt][1]+f[now][j][1])%mod; tt=j+2; if(tt<=K) f[nx][tt][0]=(f[nx][tt][0]+f[now][j][0])%mod,f[nx][tt][1]=(f[nx][tt][1]+f[now][j][1])%mod; tt=j-1; if(tt>0) f[nx][tt][0]=(f[nx][tt][0]+f[now][j][0])%mod,f[nx][tt][1]=(f[nx][tt][1]+f[now][j][1])%mod; else if(tt==0) f[nx][tt][1]=(f[nx][tt][1]+f[now][j][0]+f[now][j][1])%mod; } } for(i=0;i<=K;++i) ans=(ans+f[nx][i][1])%mod; return ans%mod; } int main(){ read(n); read(K); printf("%lld",dp()); }