loj3056/洛谷P5288 [HNOI2019]多边形

最新推荐文章于 2020-07-02 11:30:57 发布

litble

最新推荐文章于 2020-07-02 11:30:57 发布

阅读量285

点赞数 1

分类专栏：瞎搞数学文章标签：性质分析

本文链接：https://blog.csdn.net/litble/article/details/89604491

版权

数学同时被 2 个专栏收录

74 篇文章 4 订阅

订阅专栏

瞎搞

7 篇文章 0 订阅

订阅专栏

题目分析

一条边将多边形分成了“面对 $n$ 点”的一侧和“背对 $n$ 点”的一侧。“背对 $n$ 点的一侧”构成了一个点区间，这就是这条边管理的区间。因为不存在相交的边，所以每条边管理的区间之间存在嵌套关系，可以写成树结构。又发现因为整个多边形被严格分成了 $n - 3$ 个三角形，所以形成的东西是二叉树森林，其中每棵树的根代表转一下就会和 $n$ 相连的边。

最终状态一定是所有划分边（即不是构成多边形的边）都和 $n$ 相连。

最优方案一定是依次将所有划分边都转成和n相连的，并且若还存在不与 $n$ 相连的划分边，一定可以完成一次将一条边转成和n相连的操作，所以最少操作次数是 $n - 3 -$ 与 $n$ 相连的划分边数。

将二叉树森林画出来后，发现一定要转完一个点的所有祖先，转它的时候才能正好把它转成和 $n$ 相连的。

于是计算方案数也很好办了，若当前节点左子树要转 $a$ 次，右子树要转 $b$ 次，则方案数要乘以 $C _ {a+b} ^ a$ 。还要用组合数算一下不同的二叉树之间旋转方式排列好的方案数。

对于多组询问，发现若你要转一条边，它代表的点一定是其父亲的左儿子。若父亲管理的区间是
$(l, r)$ ，它管理的区间是 $(l, k)$ ，它的左儿子管理的区间是 $(l, k^{'})$ （若无左儿子， $k^{'} = k - 1$ ）。那么转一下它后，它管理的区间将变成 $(k^{'}, r)$ 。于是需要重排一下几个点的位置，得到新的树（如图）：

灵魂画手litble要退役啦

特别像一个splay的zig对吧？

然后还有若要转的边是某棵二叉树的根节点，贡献的计算方式要特别处理一下（详见代码）。

代码

#include<bits/stdc++.h>
using namespace std;
#define RI register int
int read() {
    int q=0;char ch=' ';
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9') q=q*10+ch-'0',ch=getchar();
    return q;
}
typedef pair<int,int> PR;
const int N=200005,mod=1e9+7;
int typ,n,Q,ans1,ans2,SZ,sum;
int inv[N],fac[N],ifac[N],s[N][2],sz[N],fa[N],rt[N];
vector<int> e[N];
map<PR,int> mp;

int C(int d,int u) {return 1LL*fac[d]*ifac[u]%mod*ifac[d-u]%mod;}
int ksm(int x,int y) {
    int re=1;
    for(;y;y>>=1,x=1LL*x*x%mod) if(y&1) re=1LL*re*x%mod;
    return re;
}
void init() {
    inv[0]=inv[1]=ifac[0]=fac[0]=1;
    for(RI i=1;i<=n+n;++i) fac[i]=1LL*fac[i-1]*i%mod;
    for(RI i=2;i<=n+n;++i) inv[i]=1LL*(mod-mod/i)*inv[mod%i]%mod;
    for(RI i=1;i<=n+n;++i) ifac[i]=1LL*ifac[i-1]*inv[i]%mod;
}
void prework() {
    int x,y;
    for(RI i=1;i<=n-3;++i)
        x=read(),y=read(),e[x].push_back(y),e[y].push_back(x);
    for(RI i=1;i<n;++i) e[i].push_back(i+1),e[i+1].push_back(i);
    e[1].push_back(n),e[n].push_back(1);
    for(RI i=1;i<=n;++i) sort(e[i].begin(),e[i].end());
}
void orzyyb(int &x,int las,int l,int r) {
    if(r-l<=1) return;
    int k=upper_bound(e[r].begin(),e[r].end(),l)-e[r].begin();
    k=e[r][k],x=++SZ,fa[x]=las,mp[(PR){l,r}]=x;
    orzyyb(s[x][0],x,l,k),orzyyb(s[x][1],x,k,r);
    sz[x]=sz[s[x][0]]+sz[s[x][1]]+1;
    ans2=1LL*ans2*C(sz[x]-1,sz[s[x][0]])%mod;
}

int main()
{
    typ=read(),n=read();
    init(),prework();
    ans1=n-1-(int)e[n].size(),ans2=1;
    for(RI i=1;i<(int)e[n].size();++i) {
        orzyyb(rt[i],0,e[n][i-1],e[n][i]);
        sum+=sz[rt[i]],ans2=1LL*ans2*C(sum,sz[rt[i]])%mod;
    }
    if(!typ) printf("%d\n",ans1);
    else printf("%d %d\n",ans1,ans2);
    Q=read();
    while(Q--) {
        int x=read(),y=read(); if(x>y) swap(x,y);
        int k=mp[(PR){x,y}];
        int kans1=ans1-(fa[k]?0:1)+(y==1),kans2=ans2;
        if(fa[k]) {
            int f=fa[k];
            kans2=1LL*kans2*ksm(C(sz[k]-1,sz[s[k][0]]),mod-2)%mod;
            kans2=1LL*kans2*ksm(C(sz[f]-1,sz[s[f][0]]),mod-2)%mod;
            kans2=1LL*kans2*C(sz[s[k][1]]+sz[s[f][1]],sz[s[k][1]])%mod;
            kans2=1LL*kans2*C(sz[f]-1,sz[s[k][0]])%mod;
        }
        else {
            kans2=1LL*kans2*ksm(C(sz[k]-1,sz[s[k][0]]),mod-2)%mod;
            kans2=1LL*kans2*ksm(C(sum,sz[k]),mod-2)%mod;
            kans2=1LL*kans2*C(sum-sz[k]+sz[s[k][0]],sz[s[k][0]])%mod;
            kans2=1LL*kans2*C(sum-sz[k]+sz[s[k][0]]+sz[s[k][1]],sz[s[k][1]])%mod;
        }
        if(!typ) printf("%d\n",kans1);
        else printf("%d %d\n",kans1,kans2);
    }
    return 0;
}