12.1 模拟赛

（貌似是远古大坑了

T1 bzoj 4737 组合数问题

题目大意：

给定 $n,m$ 求有多少个$C_n^m$整除$k$ $(n,m le 10^{18})$

思路：

考虑如何计算组合数 使用lucas只有在lucas递归过程中n<m才会产生0

因此我们考虑把$n,m$按照k进制分解 只要中间某一位n<m即满足题意

可以数位dp 由于统计补集比较简单 所以统计补集即中间所有数$m le n$

$dp space i space j space k$表示到了第i位 是否卡在了$n,m$的边界上 转移非常好转移

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
#define ll long long
#define inf 2139062143
#define rep(i,s,t) for(int i=(s);i<=(t);++i)
#define dwn(i,s,t) for(int i=(s);i>=(t);--i)
#define ren for(int i=fst[x];i;i=nxt[i])
#define Fill(x,t) memset(x,t,sizeof(x))
#define MAXN 100
#define MOD 1000000007
using namespace std;
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
ll n,m,T,p,dp[MAXN][2][2],a[MAXN],la,b[MAXN],lb,ans;
ll q_pow(ll bas,ll t)
{
    ll res=1;
    for(;t;t>>=1,(bas*=bas)%=MOD)
        if(t&1) (res*=bas)%=MOD;
    return res;
}
inline ll inv(ll x) {return q_pow(x,MOD-2);}
int main()
{
    T=read(),p=read();ll inv2=inv(2);int lx,ly;
    while(T--)
    {
        n=read(),m=read(),m=min(n,m),la=lb=0LL;Fill(dp,0);Fill(a,0);Fill(b,0);
        ans=(((((m%MOD)*((m+1)%MOD))%MOD)*inv2)%MOD+(((n-m+1)%MOD)*((m+1)%MOD))%MOD)%MOD;
        while(n) a[++la]=n%p,n/=p;while(m) b[++lb]=m%p,m/=p;
        dp[la+1][1][1]=1;
        dwn(i,la,1) rep(j,0,1) rep(k,0,1)
            if(dp[i+1][j][k])
            {
                if(j) lx=a[i];else lx=p-1;
                if(k) ly=b[i];else ly=p-1;
                rep(x,0,lx) rep(y,0,min(ly,x)) (dp[i][j&(x==a[i])][k&(y==b[i])]+=dp[i+1][j][k])%=MOD;
            }
        rep(i,0,1) rep(j,0,1) (ans+=MOD-dp[1][i][j])%=MOD;
        printf("%lld
",ans);
    }
}

T2 bzoj 4738 汽水

题目大意：

在树上求一条路径使得这个路径的边权平均值最接近$K$ 求这个最小差值

思路：

这种东西很容易想到点分治，我们对于所有初始边权都减去$K$，这样就是求最接近$0$的路径

对于分治重心内的所有路径，维护三个值：路径总长，路径上边数，来自哪一个子树

我们对于总长排序后可以二分一个值，判断能否存在两个路径拼起来差值满足mid

把式子转化为$sum-num imes x$其中$x$为二分出的mid 这样只要分别$check (pm mid)$

写两个check即可

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<complex>
#include<iomanip>
#define Fill(a,x) memset(a,x,sizeof(a))
#define rep(i,s,t) for(register int i=(s),i__end=(t);i<=i__end;++i)
#define dwn(i,s,t) for(register int i=(s),i__end=(t);i>=i__end;--i)
#define ren for(int i=fst[x];i;i=nxt[i])
#define inf 1LL<<60
#define ll long long
#define ull unsigned long long
#define MAXN 50100
#define MOD 998244353
using namespace std;
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct data {int num,bl;ll val;}g[MAXN],A,B;
bool operator < (data a,data b) {return a.val<b.val;}
int n,fst[MAXN],nxt[MAXN<<1],to[MAXN<<1],cnt,sz[MAXN],mx[MAXN],Sum,Mx,rt,len,vis[MAXN];
ll val[MAXN<<1],K,ans=inf;
void add(int u,int v,ll w) {nxt[++cnt]=fst[u],fst[u]=cnt,to[cnt]=v,val[cnt]=w;}
void getrt(int x,int pa)
{
    mx[x]=0,sz[x]=1;ren if(to[i]^pa&&!vis[to[i]]) getrt(to[i],x),sz[x]+=sz[to[i]],mx[x]=max(mx[x],sz[to[i]]);
    mx[x]=max(mx[x],Sum-sz[x]);if(mx[x]<Mx) Mx=mx[x],rt=x;
}
void Get(int x,int pa,int dep,int frm,ll dis)
{
    g[++len]=(data){dep,frm,dis};
    ren if(to[i]^pa&&!vis[to[i]]) Get(to[i],x,dep+1,frm,dis+val[i]);
}
int pos;
void upd(data x)
{
    if(B<x) return ;
    if(x<A) {if(x.bl!=A.bl) B=A;A=x;}else if(x.bl!=A.bl) B=x;
}
int cheq(ll x)
{
    A=B=(data){0,0,inf};int j=pos-1;
    rep(i,pos,len)
    {
        for(;j&&g[i].val+g[j].val>=0;j--) upd((data){0,g[j].bl,g[j].val-x*g[j].num});
        if((A.bl==g[i].bl?B.val:A.val)<(ll)x*g[i].num-g[i].val) return 1;
        upd((data){0,g[i].bl,g[i].val-x*g[i].num});
    }
    return 0;
}
int cheQ(ll x)
{
    A=B=(data){0,0,inf};int j=pos;
    dwn(i,pos-1,1)
    {
        for(;j<=len&&g[i].val+g[j].val<0;j++) upd((data){0,g[j].bl,(ll)x*g[j].num-g[j].val});
        if((A.bl==g[i].bl?B.val:A.val)<g[i].val-x*g[i].num) return 1;
        upd((data){0,g[i].bl,(ll)x*g[i].num-g[i].val});
    }
    return 0;
}
void calc(int x)
{
    vis[x]=1,len=0;g[++len]=(data){0,0,0};ren if(!vis[to[i]]) Get(to[i],x,1,to[i],val[i]);
    sort(g+1,g+len+1);
    pos=lower_bound(g+1,g+len+1,(data){0,0,0})-g;ll l=1,r=ans;
    for(ll mid;mid=l+r>>1,l<=r;) if(cheq(mid)||cheQ(-mid)) r=mid-1;else l=mid+1;
    ans=min(ans,l-1);
}
void div(int x)
{
    calc(x);ren if(!vis[to[i]]) Sum=sz[to[i]],Mx=n+1,getrt(to[i],x),div(rt);
}
int main()
{
    n=read(),K=read();int a,b;ll c;
    rep(i,2,n) a=read(),b=read(),c=read(),ans=min(ans,abs(c-K)),add(a,b,c-K),add(b,a,c-K);
    Sum=n,Mx=n+1;getrt(1,0);div(rt);printf("%lld
",ans);
}

T3 定向越野（由于是计算几何所以咕咕咕了