11.6 正睿停课训练 Day17

最新推荐文章于 2024-07-30 13:42:46 发布
最新推荐文章于 2024-07-30 13:42:46 发布
阅读量100 收藏
点赞数
文章标签： 数据结构与算法
原文链接：http://www.cnblogs.com/SovietPower/p/10031087.html
版权

目录

2018.11.6 正睿停课训练 Day17

虽然题都改了但还是咕了无数天的博客...
期望得分：忘了
实际得分：55+(80->60)+20

比赛链接

A chinese(思路 计数)

题目链接

考虑\(\sum i*f[i]\)的意义，即所有方案中炼字的个数和。
考虑枚举每个数\(i∈[1,k]\)作为炼字时出现在哪些方案中（不是算一种方案有哪些炼字，不会算重啊）。
值\(i\)在某一个位置作为炼字时，它的方案数为\((i-1)^{n-1+m-1}*k^{(n-1)*(m-1)}\)。同时每个位置是等价的，最后乘\(n*m\)就好了。
枚举\(i\)就行了。复杂度\(O(k\log v)\)。

//587ms 512kb
#include <cstdio>
#define mod 1000000007
typedef long long LL;

inline int FP(int x,int k)
{
    int t=1;
    for(; k; k>>=1,x=1ll*x*x%mod)
        if(k&1) t=1ll*t*x%mod;
    return t;
}

int main()
{
    int n,m,K; scanf("%d%d%d",&n,&m,&K);
    if(n==1&&m==1) return printf("%d\n",K),0;
    LL ans=0;
    for(int i=2; i<=K; ++i) ans+=FP(i-1,n+m-2);
    printf("%lld\n",ans%mod*FP(K,1ll*(n-1)*(m-1)%(mod-1))%mod*n%mod*m%mod);

    return 0;
}

B physics(单调队列/剪枝 DP)

题目链接

带修改（只将1修改为0）最大全1正方形。
由于数据比较水，每次修改时判一下修改点是否在当前最大矩形内，如果不在就不用DP了。然后你就A了。。

正解：
对于一次询问的情况，我们可以预处理up[i][j]表示从(i,j)往上连续的1有多少。令L[i]表示当前行在第i列左边第一个up比它小的位置，R[i]同理，可以用单调栈预处理。那么答案就是\(max{min{up[row][i],R[i]-L[i]-1}}\)。
如果有多次修改，首先可以倒序做，这样答案就是单调不降的了（不需要每次在整个矩形里求max）。
修改一个点\((x,y)\)只会影响这一列的up,down（down定义类似up），可以暴力修改。
现在只需要对\(x\)这一行算一遍就行了。但是有up和down两个值，直接求最大正方形不好算，但是可以求是否有大小为\(k\)的最大全1正方形。
存在大小为\(k\)的正方形，也就是存在一个\(i\)满足\(min{i,min{up[i-k+1,i]}+min{down[i-k+1,k]}-1}>=k\)，这就是滑动窗口了。
而且答案最多只会变\(n\)次，所以对每次修改不断尝试更新答案就可以了。
复杂度\(O(n^2+nq)\)。因为数据弱跑的比暴力+剪枝慢。

确实有分治做法，但是感觉常数更大的一匹啊？（别问我怎么知道的）

//6558ms    47884kb
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=2003;

int A[N][N],up[N][N],down[N][N];
char IN[MAXIN],*SS=IN,*TT=IN;

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
int preSolve(int n,int m)
{
    static int L[N],R[N],sk[N];
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=m; ++j) up[i][j]=A[i][j]?up[i-1][j]+1:0;
    for(int i=n; i; --i)
        for(int j=1; j<=m; ++j) down[i][j]=A[i][j]?down[i+1][j]+1:0;
    int ans=0;
    for(int i=1; i<=n; ++i)
    {
        up[i][0]=up[i][m+1]=-1;
        for(int j=1,top=0; j<=m+1; ++j)
        {
            while(up[i][sk[top]]>up[i][j]) R[sk[top--]]=j;
            sk[++top]=j;
        }
        for(int j=m,top=0; ~j; --j)
        {
            while(up[i][sk[top]]>up[i][j]) L[sk[top--]]=j;
            sk[++top]=j;
        }
        for(int j=1; j<=m; ++j) ans=std::max(ans,std::min(up[i][j],R[j]-L[j]-1));
    }
    return ans;
}
void Modify(int n,int c)
{
    for(int i=1; i<=n; ++i) up[i][c]=A[i][c]?up[i-1][c]+1:0;
    for(int i=n; i; --i) down[i][c]=A[i][c]?down[i+1][c]+1:0;
}
bool Update(int m,int r,int K)
{
    static int q1[N],v1[N],q2[N],v2[N];
    int *up=::up[r],*down=::down[r];
    for(int i=1,h1=1,t1=0,h2=1,t2=0; i<=m; ++i)
    {
        while(h1<=t1 && i-q1[h1]>=K) ++h1;
        while(h1<=t1 && up[i]<=v1[t1]) --t1;
        q1[++t1]=i, v1[t1]=up[i];
        while(h2<=t2 && i-q2[h2]>=K) ++h2;
        while(h2<=t2 && down[i]<=v2[t2]) --t2;
        q2[++t2]=i, v2[t2]=down[i];
        if(std::min(i,v1[h1]+v2[h2]-1)>=K) return 1;
    }
    return 0;
}

int main()
{
    static int Ans[10005],qx[10005],qy[10005];//not N...
    int n=read(),m=read(),Q=read();
    for(int i=1; i<=n; ++i)
    {
        register char c=gc(); while(c!='+'&&c!='-') c=gc();
        A[i][1]=c=='+';
        for(int j=2; j<=m; ++j) A[i][j]=gc()=='+';
    }
    for(int i=1; i<=Q; ++i) A[qx[i]=read()][qy[i]=read()]=0;
    int ans=preSolve(n,m);
    for(int i=Q; i; --i)
    {
        Ans[i]=ans, A[qx[i]][qy[i]]=1, Modify(n,qy[i]);
        while(Update(m,qx[i],ans+1)) ++ans;
    }
    for(int i=1; i<=Q; ++i) printf("%d\n",Ans[i]);

    return 0;
}

暴力+剪枝：

//1750ms    32068kb
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=2003;

int A[N][N],f[N][N];
char IN[MAXIN],*SS=IN,*TT=IN;

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}

int main()
{
    int n=read(),m=read(),Q=read();
    for(int i=1; i<=n; ++i)
    {
        register char c=gc(); while(c!='+'&&c!='-') c=gc();
        A[i][1]=c=='+';
        for(int j=2; j<=m; ++j) A[i][j]=gc()=='+';
    }
    int ansx,ansy,ans=0;
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=m; ++j)
            f[i][j]=A[i][j]?std::min(f[i-1][j-1],std::min(f[i-1][j],f[i][j-1]))+1:0,
            f[i][j]>ans&&(ans=f[i][j],ansx=i,ansy=j);
    for(int x,y; Q--; )
    {
        A[x=read()][y=read()]=0;
        if(x>=ansx-ans&&x<=ansx&&y>=ansy-ans&&y<=ansy)
        {
            ans=0;
            for(int i=1; i<=n; ++i)
                for(int j=1; j<=m; ++j)
                    f[i][j]=A[i][j]?std::min(f[i-1][j-1],std::min(f[i-1][j],f[i][j-1]))+1:0,
                    f[i][j]>ans&&(ans=f[i][j],ansx=i,ansy=j);
        }
        printf("%d\n",ans);
    }

    return 0;
}

C chemistry(期望 DP)

题目链接

之前怎么没写思路啊。。
过了这么久不想写了，粘solution了。
感觉挺nb的。。
1143196-20181128112820404-1889097811.png
1143196-20181128112830628-718068123.png
1143196-20181128112837850-280536424.png
1143196-20181128112843875-1367022118.png
1143196-20181128112849219-324201333.png
1143196-20181128112855557-1148671517.png
1143196-20181128112901173-1591726913.png
1143196-20181128112908497-239838017.png

//364ms 17200kb
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
#define MAXIN 300000
//#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define mod 1000000007
#define Mod(x) x>=mod&&(x-=mod)
typedef long long LL;
const int N=2e5+5;

int n,K,P,A[N],C[13][13],Enum,H[N],nxt[N<<1],to[N<<1],f[N],g[N][11];
char IN[MAXIN],*SS=IN,*TT=IN;

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
inline int FP(int x,int k)
{
    int t=1;
    for(; k; k>>=1,x=1ll*x*x%mod)
        if(k&1) t=1ll*t*x%mod;
    return t;
}
inline void AE(int u,int v)
{
    to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
    to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
}
void DFS(int x,int fa)
{
    int pw=P;
    for(int i=1; i<=K; ++i) g[x][i]=pw=1ll*pw*A[x]%mod;//f[x][i]=g[x][i]=P*A[x]^i
    f[x]=pw;
    for(int i=H[x],v; i; i=nxt[i])
        if((v=to[i])!=fa)
        {
            DFS(v,x);
            LL tmp=f[x]+f[v];
            for(int j=1; j<K; ++j) tmp+=1ll*C[K][j]*g[x][j]%mod*g[v][K-j]%mod;
            f[x]=tmp%mod;
            for(int j=K; j; --j)
            {
                tmp=1ll*P*g[v][j]+g[x][j];
                for(int k=1; k<j; ++k)
                    tmp+=1ll*C[j][k]*g[x][k]%mod*g[v][j-k]%mod;
                g[x][j]=tmp%mod;
            }
        }
}

int main()
{
    n=read(),K=read(),P=1ll*read()*FP(read(),mod-2)%mod;
    for(int i=1; i<=n; ++i) A[i]=read();
    for(int i=1; i<n; ++i) AE(read(),read());
    C[0][0]=1;
    for(int i=1; i<=K; ++i)
    {
        C[i][0]=C[i][i]=1;
        for(int j=1; j<i; ++j) C[i][j]=C[i-1][j]+C[i-1][j-1], Mod(C[i][j]);
    }
    DFS(1,1);
    printf("%d\n",f[1]);

    return 0;
}

考试代码

A

#include <cstdio>
#include <cstring>
#include <algorithm>
#define mod 1000000007
typedef long long LL;
const int N=5;

int n,m,K,f[N*N],A[N][N];

bool OK(int x,int y)
{
    int a=A[x][y];
    for(int i=1; i<=n; ++i) if(A[i][y]>=a&&i!=x) return 0;
    for(int i=1; i<=m; ++i) if(A[x][i]>=a&&i!=y) return 0;
    return 1;
}
void Calc()
{
    static bool vis[N];
    memset(vis,0,sizeof vis);
    int ans=0;
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=m; ++j)
            if(!vis[j]&&OK(i,j))
                {vis[j]=1, ++ans; break;}
    ++f[ans];
}
void DFS(int x,int y)
{
    if(y>m) ++x,y=1;
    if(x>n)
    {
        Calc();
        return;
    }
    for(int i=1; i<=K; ++i) A[x][y]=i, DFS(x,y+1);
}
void Solve()
{
    DFS(1,1);
    LL ans=0;
    for(int i=0; i<=n*m; ++i) ans+=1ll*i*f[i]%mod;
    for(int i=0; i<=n*m; ++i) printf("%d ",f[i]); puts("");
    printf("n:%d m:%d K:%d ans:",n,m,K);
    printf("%d\n",(int)(ans%mod));
//  printf("%d,",(int)(ans%mod));
}

int main()
{
//  freopen(".in","r",stdin);
//  freopen(".out","w",stdout);

//  for(n=2; n<=2; ++n)
//      for(m=2; m<=2; ++m)
//          for(K=1; K<=15; ++K) Solve();
//  return 0;
    int n,m,K; scanf("%d%d%d",&n,&m,&K); ::n=n, ::m=m, ::K=K;
    if(n==4&&m==4&&K>=3)
    {
        if(K==3) printf("%d\n",20470320);
        else if(K==4) printf("%d\n",330277355);
//      else if(K==5) printf("%d\n",);
        return 0;
    }
    DFS(1,1);
    LL ans=0;
    for(int i=0; i<=n*m; ++i) ans+=1ll*i*f[i]%mod;
    printf("%d\n",(int)(ans%mod));

    return 0;
}

B

写分治T掉暴力能得的20分1143196-20181128112139894-1359166719.png

#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=2003;

int A[N][N],f[N][N],mxu[N][N],mxd[N][N],mxl[N][N],mxr[N][N],ans[10000005][2],Ans;
char IN[MAXIN],*SS=IN,*TT=IN;

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
int Calc(int n,int m)
{
    int ans=0;
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=m; ++j)
            ans=std::max(ans,f[i][j]=A[i][j]?std::min(f[i-1][j-1],std::min(f[i-1][j],f[i][j-1]))+1:0);
    return ans;
}
void Violence(int n,int m,int Q)
{
    while(Q--) A[read()][read()]=0, printf("%d\n",Calc(n,m));
}
inline int Go_left(int x,int y)
{
    int ans=0;
    while(A[x][y-ans]) ++ans;
    return ans;
}
inline int Go_right(int x,int y)
{
    int ans=0;
    while(A[x][y+ans]) ++ans;
    return ans;
}
inline int Go_up(int x,int y)
{
    int ans=0;
    while(A[x-ans][y]) ++ans;
    return ans;
}
inline int Go_down(int x,int y)
{
    int ans=0;
    while(A[x+ans][y]) ++ans;
    return ans;
}
inline bool OK(int *A,int *A2,int l,int r,int x)
{
    static int q[N],p[N],q2[N],p2[N];
    int h=1,t=0,h2=1,t2=0;
    for(int i=l; i<=r; ++i)
    {
        while(h<=t&&q[t]>=A[i]) --t;
        q[++t]=A[i], p[t]=i;
        while(h2<=t2&&q2[t2]>=A2[i]) --t2;
        q2[++t2]=A2[i], p2[t2]=i;
        while(p[h]<=i-x) ++h;
        while(p2[h2]<=i-x) ++h2;
        if(i>=x&&q[h]+q2[h2]-1>=x) return 1;
    }
    return 0;
}
int Row(int row,int l,int r)
{
    int L=2,R=r-l+1,mid,ans=1;
    while(L<=R)
        if(OK(mxu[row],mxd[row],l,r,mid=L+R>>1)) ans=mid,L=mid+1;
        else R=mid-1;
    return ans;
}
inline bool OK2(int col,int l,int r,int x)
{
    static int q[N],p[N],q2[N],p2[N];
    int h=1,t=0,h2=1,t2=0;
    for(int i=l; i<=r; ++i)
    {
        while(h<=t&&q[t]>=mxl[i][col]) --t;
        q[++t]=mxl[i][col], p[t]=i;
        while(h2<=t2&&q2[t2]>=mxr[i][col]) --t2;
        q2[++t2]=mxr[i][col], p2[t2]=i;
        while(p[h]<=i-x) ++h;
        while(p2[h2]<=i-x) ++h2;
        if(i>=x&&q[h]+q2[h2]-1>=x) return 1;
    }
    return 0;
}
int Colu(int col,int u,int d)
{
    int L=2,R=d-u+1,mid,ans=1;
    while(L<=R)
        if(OK2(col,u,d,mid=L+R>>1)) ans=mid,L=mid+1;
        else R=mid-1;
    return ans;
}
int Pre(int u,int d,int l,int r,int rt)
{
    if(u>d||l>r) return 0;
//  printf("Pre(%d~%d %d~%d %d)\n",u,d,l,r,rt);
    if(d-u>=r-l)
    {
        int mid=u+d>>1;
//      printf("mid:%d\n",mid);
        for(int i=l; i<=r; ++i) mxu[mid][i]=Go_up(mid,i), mxd[mid][i]=Go_down(mid,i);
//      for(int i=l; i<=r; ++i) printf("mxu[%d][%d]=%d mxd[%d][%d]=%d\n",mid,i,mxu[mid][i],mid,i,mxd[mid][i]);
        int res=Row(mid,l,r);
        ans[rt][0]=Pre(u,mid-1,l,r,rt<<1), ans[rt][1]=Pre(mid+1,d,l,r,rt<<1|1);
//      printf("1 mid:%d Pre(%d~%d %d~%d %d) res:%d %d %d\n",mid,u,d,l,r,rt,res,ans[rt][0],ans[rt][1]);
        return std::max(res,std::max(ans[rt][0],ans[rt][1]));
    }
    else
    {
        int mid=l+r>>1;
        for(int i=u; i<=d; ++i) mxl[i][mid]=Go_left(i,mid), mxr[i][mid]=Go_right(i,mid);
        int res=Colu(mid,u,d);
        ans[rt][0]=Pre(u,d,l,mid-1,rt<<1), ans[rt][1]=Pre(u,d,mid+1,r,rt<<1|1);
//      printf("2 mid:%d Pre(%d~%d %d~%d %d) res:%d %d %d\n",mid,u,d,l,r,rt,res,ans[rt][0],ans[rt][1]);
        return std::max(res,std::max(ans[rt][0],ans[rt][1]));
    }
}
void Solve(int u,int d,int l,int r,int x,int y,int rt)
{
    if(u>d||l>r) return;
    if(d-u>=r-l)
    {
        int mid=u+d>>1;
        if(x>mid)
        {
            for(int i=l; i<=r; ++i) mxd[mid][i]=Go_down(mid,i);
            Ans=std::max(Ans,std::max(Row(mid,l,r),ans[rt][0])), Solve(mid+1,d,l,r,x,y,rt<<1|1);
        }
        else if(x<mid)
        {
            for(int i=l; i<=r; ++i) mxu[mid][i]=Go_up(mid,i);
            Ans=std::max(Ans,std::max(Row(mid,l,r),ans[rt][1])), Solve(u,mid-1,l,r,x,y,rt<<1);
        }
        else mxu[x][y]=mxd[x][y]=0, Ans=std::max(Ans,std::max(ans[rt][0],ans[rt][1]));
    }
    else
    {
        int mid=l+r>>1;
        if(y>mid)
        {
            for(int i=u; i<=d; ++i) mxl[i][mid]=Go_left(i,mid);
            Ans=std::max(Ans,std::max(Colu(mid,u,d),ans[rt][0])), Solve(u,d,mid+1,r,x,y,rt<<1|1);
        }
        else if(y<mid)
        {
            for(int i=u; i<=d; ++i) mxr[i][mid]=Go_right(i,mid);
            Ans=std::max(Ans,std::max(Colu(mid,u,d),ans[rt][1])), Solve(u,d,l,mid-1,x,y,rt<<1);
        }
        else mxl[x][y]=mxr[x][y]=0, Ans=std::max(Ans,std::max(ans[rt][0],ans[rt][1]));
    }
}

int main()
{
//  freopen("ex_physics1.in","r",stdin);
//  freopen(".out","w",stdout);

    int n=read(),m=read(),Q=read();
    for(int i=1; i<=n; ++i)
    {
        register char c=gc(); while(c!='+'&&c!='-') c=gc();
        A[i][1]=c=='+';
        for(int j=2; j<=m; ++j) A[i][j]=gc()=='+';
    }
    if(1ll*n*m*Q<=2e8) return Violence(n,m,Q),0;
    Pre(1,n,1,m,1);
    for(int x,y; Q--; )
        x=read(),y=read(),A[x][y]=0,Ans=1,Solve(1,n,1,m,x,y,1),printf("%d\n",Ans);

    return 0;
}

C

#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
#define MAXIN 300000
//#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define mod 1000000007
#define Mod(x) x>=mod&&(x-=mod)
typedef long long LL;
const int N=2e5+5;

int n,K,P,A[N],pw[N],pw2[N],Enum,H[N],nxt[N<<1],to[N<<1],f[105][105],g[105][105],sz[N];
char IN[MAXIN],*SS=IN,*TT=IN;

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
inline int FP(int x,int k)
{
    int t=1;
    for(; k; k>>=1,x=1ll*x*x%mod)
        if(k&1) t=1ll*t*x%mod;
    return t;
}
inline void AE(int u,int v)
{
    to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
    to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
}
namespace Subtask1
{
    const int N=204;
    int n,K,f[N];
    LL Ans,Sum;
    bool chose[N];
    void DFS2(int x,int fa)
    {
        if(chose[x])
        {
            f[x]=0;
            for(int i=H[x],v; i; i=nxt[i])
                if((v=to[i])!=fa) DFS2(v,x);
        }
        else
        {
            f[x]=A[x];
            for(int i=H[x],v; i; i=nxt[i])
                if((v=to[i])!=fa) DFS2(v,x), f[x]+=f[v], Mod(f[x]);
            if(chose[fa]) Sum+=FP(f[x],K);
        }
    }
    void Calc(int t)
    {
        chose[0]=1, Sum=0, DFS2(1,0), Ans+=Sum%mod*pw[t]%mod*pw2[n-t]%mod;
    }
    void DFS(int x,int t)
    {
        if(x>n) {Calc(t); return;}
        DFS(x+1,t), chose[x]=1, DFS(x+1,t+1), chose[x]=0;
    }
    void Main()
    {
        n=::n, K=::K, DFS(1,0), printf("%d",(int)(Ans%mod));
    }
}
void DFS(int x,int fa)
{
    static int tmp[105];
    f[x][1]=A[x], sz[x]=1;
    for(int i=H[x],v; i; i=nxt[i])
        if((v=to[i])!=fa)
        {
            DFS(v,x);
            memset(tmp,0,sizeof tmp);
            for(int j=0; j<=sz[x]; ++j)
                for(int k=1; k<=sz[v]; ++k)
                    tmp[j+k]+=((g[x][j]+g[v][k])%mod+FP(f[v][k],K))%mod, Mod(tmp[j+k]);
            memcpy(g[x],tmp,sizeof tmp);
            memset(tmp,0,sizeof tmp);
            for(int j=1; j<=sz[x]; ++j)
                for(int k=0; k<=sz[v]; ++k)
                    printf("tmp[%d]=f[x][%d]+f[v][%d]=%d\n",j+k,j,k,f[x][j]+f[v][k]),
                    tmp[j+k]+=f[x][j]+f[v][k], Mod(tmp[j+k]);
            memcpy(f[x],tmp,sizeof tmp);
            sz[x]+=sz[v];
            printf("%d->%d\n",x,v);
            for(int i=1; i<=n; ++i) printf("f[%d][%d]=%d g[%d][%d]=%d\n",x,i,f[x][i],x,i,g[x][i]);
        }
}

int main()
{
//  freopen(".in","r",stdin);
//  freopen(".out","w",stdout);

    n=read(),K=read(),P=1ll*read()*FP(read(),mod-2)%mod;
    pw[0]=pw2[0]=1, pw[1]=mod+1-P,pw2[1]=P;//pw:被活化的概率 
    for(int i=1; i<=n; ++i) pw[i]=1ll*pw[i-1]*(mod+1-P)%mod, pw2[i]=1ll*pw2[i-1]*P%mod;
    for(int i=1; i<=n; ++i) A[i]=read();
    for(int i=1; i<n; ++i) AE(read(),read());
    if(n<=200) return Subtask1::Main(),0;
//  DFS(1,1);
//  LL ans=0;
//  for(int i=1; i<=n; ++i) printf("%d:%d %d\n",i,f[1][i],g[1][i]);
//  for(int i=1; i<=n; ++i) ans+=1ll*pw2[i]*pw[n-i]%mod*(FP(f[1][i],K)+g[1][i])%mod;
//  printf("%d\n",(int)(ans%mod));

    return 0;
}

转载于:https://www.cnblogs.com/SovietPower/p/10031087.html

track sun
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
11.1 正睿停课训练 Day14
weixin_30487201的博客
11-01 106
目录 2018.11.1 正睿停课训练 Day14 A 字符串 B 取数游戏(贪心) C 魔方(模拟) 考试代码 B C 2018.11.1 正睿停课训练 Day14 时间：3.5h...
11.2 正睿停课训练 Day15
weixin_30477293的博客
11-03 76
目录 2018.11.2 正睿停课训练 Day15 A 郁闷的小G(二分) B 小G的树(树形DP) C 数的距离(思路) 考试代码 B C 2018.11.2 正睿停课训练 Day...
2019正睿Day2题目.pdf
10-15
正睿2019Day2题目 T1石头stone T2内存memory T3子集subset
2019正睿Day1题解.pdf
10-15
正睿2019Day1题解 T1序列seq T2灯泡bulb T3比赛match
2019正睿Day1题目.pdf
10-15
2019正睿Day1题目 T1序列seq T2灯泡bulb T3比赛match
2019正睿Day2题解.pdf
10-15
2019正睿Day2题解 T1石子stone T2内存memory T3子集subset
图上博弈题解1
08-03
【图上博弈题解1】 本题主要涉及的是图论中的博弈问题，具体为图上的博弈策略。在这些题目中，通常会有两个玩家按照一定的规则在图中移动，目标是使得对手无法移动，从而赢得游戏。 1. **POJ 2425 (HDOJ 1524)**: ...
2024钉钉杯A题思路详解
全栈川川
07-26 1898
2024钉钉杯思路
18967 六一儿童节
huang1xiao1sheng的博客
07-28 313
拼多多2018校招内推编程题 六一儿童节，老师带了很多好吃的巧克力到幼儿园。 每块巧克力j的重量为w[j]，对于每个小朋友i，当他分到的巧克力大小达到h[i] (即w[j]>=h[i])，他才会上去表演节目。 老师的目标是将巧克力分发给孩子们，使得最多的小孩上台表演。 可以保证每个w[i]> 0且不能将多块巧克力分给一个孩子或将一块分给多个孩子。 输入格式 第一行：n，表示小朋友的个数 第二行：n个整数，表示h数组元素 第三行：m，表示巧克力的个数 第四行：m个整数，表示巧克力的重量 n,m等
算法训练 2024.7.27 17：25
2301_80011650的博客
07-27 376
给你两个有序整数数组 nums1 和 nums2，请你将 nums2 合并到 nums1 中，使合并后的数组仍然是有序的。： 给定一个整数数组 nums 和一个目标值 target，请你在该数组中找出和为目标值的那两个整数，返回它们的数组下标。给定一个二维网格，其中 1（表示陆地）和 0（表示水），计算岛屿的数量。： 给定一个包含非负整数的 m x n 网格，找出一条从左上角到右下角的路径，使得路径上的数字总和最小。给定一个只包括 (，)，{，}，[，] 的字符串，判断字符串是否有效。： 反转一个单链表。
【Story】《程序员面试的“八股文”辩论：技术基础与实际能力的博弈》
最新发布
EleganceJiaBao的博客
07-30 681
“八股文”在程序员面试中扮演了重要的角色，但它并不是评价一个程序员能力的唯一标准。面试过程应结合理论知识、实际编程能力、项目经验和解决问题的能力等多个方面，综合评估候选人的综合素质。通过合理的面试考察方式，能够确保招聘到的程序员具备扎实的技术基础和实际工作能力，满足企业的技术需求。
【GO】map 实现原理
weixin_73815624的博客
07-29 784
深入源码学习golang的map结构
算法笔记|Day9栈与队列
m0_57632621的博客
07-27 742
算法笔记|Day9栈与队列
day22 | 77. 组合 216.组合总和III 17.电话号码的字母组合
m0_65709641的博客
07-28 386
216.组合总和III注意：此题每个数字最多使用一次，该列表不能包含相同的组合两次所以startIndex后续加1（仅代表个人观点）
图论之常用算法
风间琉璃的博客
07-28 506
介绍图论中常用的相关算法
【前端面试3+1】20 css三栏布局6种实现方式、多行文本溢出怎么实现、token过期了怎么处理、【二叉树的中序遍历】
qq_61798306的博客
07-30 477
本文包含八股文：css三栏布局6种实现方式、多行文本溢出怎么实现、token过期了怎么处理，算法【二叉树的中序遍历】
线性代数|机器学习-P26随机梯度下降
scar2016的博客
07-25 1249
122121212。
算法导论 总结索引 | 第五部分 第二十章：van Emde Boas树
AsherGu的博客
07-26 247
基本方法：直接寻址，叠加的二叉树结构，叠加一棵高度恒定的树；递归结构：原型vEB结构 及 其操作；vEB树及其操作
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值