luogu P5288 [HNOI2019]多边形

传送门

这是什么神仙操作...

首先要注意一些性质.首先每一个((x,n))的边可以把当前多边形分成两半,这两半的操作是独立的.然后对于某一个没有((x,n))的边的多边形,最优操作是唯一的.拿样例举例,必须先选((1,5)),然后多边形被分成两半,这两半分别只能选((1,3)),((3,5)).

可以发现,这里每次操作的多边形上的点分别是(l,l+1...,r(lge 1,r<n,r-l>1))和(n),然后最优方案要选((l,r))这条边,把多边形分成({[l,mid],n})和({[mid,r],n})((mid)是当前多边形中的(l)和(r)除了(n)以外的那个公共点)两部分.根据这个我们可以建出若干二叉树,每个点的左右两个儿子是这个多边形操作之后的左右两半多边形(不能是三角形,因为三角形不用操作)

考虑计算答案,第一问答案为这些二叉树节点个数(也就是(n-1-)((x,n))边条数),因为每个节点都是一次操作;第二问答案,因为这个二叉树每个节点要在他的父亲之后才能操作,所以如果把每个点看成一个排列中的数,就是要求满足(x)在(fa_x)后面的排列个数,这个直接把SAO蒯过来就好了对于每个节点的子树,首先这个点要放在第一位,两个儿子子树的序列之间是没有先后关系的,所以一个节点对答案的贡献是(inom{sz_{ch_{x,0}}+sz_{ch_{x,1}}}{sz_{ch_{x,0}}}),把这些贡献乘起来就好了.然后对于所有二叉树的根,我们把他们向一个超级根连边,因为相互没有影响,所以答案还要乘上把所有树合并的方案,类似于树型背包,每次把一个并到总集合,然后类似于上面的组合数即可

现在可能会有一些翻转边的操作.手玩可以发现这相当于对于那条边对应的点进行一次(rotate)操作.然后如果这个点不是二叉树的根,那么只要先把这个点和父亲的贡献扣掉,(rotate)后在加上新的贡献即可;否则这个点是代表着初始情况最优操作中的一个,第一问答案要减去1,然后第二问答案先要扣掉这个点贡献和这个点子树对超级根的贡献,然后因为这个点被操作了,就相当于删掉这个点,两个儿子分别成为一个新二叉树,那么贡献加上这两个儿子子树并到超级根的贡献就好了

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<ctime>
#include<queue>
#include<map>
#include<set>
#define LL long long
#define uLL unsigned long long 

using namespace std;
const int N=1e5+10,mod=1e9+7;
int rd()
{
    int x=0,w=1;char ch=0;
    while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x*w;
}
map<int,int> id[N];
vector<int> e[N];
int fac[N],iac[N];
int fpow(int a,int b){int an=1;while(b){if(b&1) an=1ll*an*a%mod;a=1ll*a*a%mod,b>>=1;} return an;}
int inv(int a){return fpow(a,mod-2);}
int C(int n,int m){return m<0||n<m?0:1ll*fac[n]*iac[m]%mod*iac[n-m]%mod;}
int w,n,fa[N],ch[N][2],sz[N],size,rt[N],tp,ans=1;
void bui(int &x,int l,int r)
{
    if(l+1==r) return;
    x=id[l][r];
    sz[x]=1;
    int mid=*(--lower_bound(e[l].begin(),e[l].end(),r));
    bui(ch[x][0],l,mid),bui(ch[x][1],mid,r);
    if(ch[x][0]) fa[ch[x][0]]=x;
    if(ch[x][1]) fa[ch[x][1]]=x;
    sz[x]+=sz[ch[x][0]]+sz[ch[x][1]];
    ans=1ll*ans*C(sz[ch[x][0]]+sz[ch[x][1]],sz[ch[x][0]])%mod;
}

int main()
{
    fac[0]=1;
    for(int i=1;i<=N-10;++i) fac[i]=1ll*fac[i-1]*i%mod;
    iac[N-10]=fpow(fac[N-10],mod-2);
    for(int i=N-10;i;--i) iac[i-1]=1ll*iac[i]*i%mod;
    w=rd(),n=rd();
    for(int i=1;i<n;++i) e[i].push_back(i+1),e[i+1].push_back(i);
    e[n].push_back(1),e[1].push_back(n);
    for(int i=1;i<=n-3;++i)
    {
        int x=rd(),y=rd();
        e[x].push_back(y),e[y].push_back(x);
        id[x][y]=id[y][x]=i;
    }
    for(int i=1;i<=n;++i) sort(e[i].begin(),e[i].end());
    int nn=e[n].size(),ff=n-1-nn;
    for(int i=0;i<nn-1;++i)
    {
        int x=e[n][i],y=e[n][i+1];
        if(x+1==y) continue;
        ++tp,bui(rt[tp],x,y);
    }
    for(int i=1;i<=tp;++i)
        size+=sz[rt[i]],ans=1ll*ans*C(size,sz[rt[i]])%mod;
    printf("%d ",ff);
    if(w) printf("%d",ans);
    puts("");
    int q=rd();
    while(q--)
    {
        int x=id[rd()][rd()],na=ans;
        if(fa[x])
        {
            int yy=ch[fa[x]][1]==x;
            na=1ll*na*inv(C(sz[ch[x][0]]+sz[ch[x][1]],sz[ch[x][0]]))%mod*inv(C(sz[ch[fa[x]][0]]+sz[ch[fa[x]][1]],sz[ch[fa[x]][0]]))%mod;
            na=1ll*na*C(sz[ch[x][yy^1]]+sz[ch[fa[x]][yy^1]],sz[ch[x][yy^1]])%mod*C(sz[ch[x][yy]]+sz[ch[x][yy^1]]+sz[ch[fa[x]][yy^1]]+1,sz[ch[x][yy]])%mod;
        }
        else
        {
            na=1ll*na*inv(C(sz[ch[x][0]]+sz[ch[x][1]],sz[ch[x][0]]))%mod*inv(C(size,sz[x]))%mod;
            na=1ll*na*C(size-sz[x]+sz[ch[x][0]],sz[ch[x][0]])%mod*C(size-1,sz[ch[x][1]])%mod;
        }
        printf("%d ",ff-(!fa[x]));
        if(w) printf("%d",na);
        puts("");
    }
    return 0;
}