【BZOJ3456】轩辕朗的城市规划 无向连通图计数 CDQ分治 FFT 多项式求逆 多项式ln...

最新推荐文章于 2022-06-01 19:14:24 发布
最新推荐文章于 2022-06-01 19:14:24 发布
阅读量150 收藏
点赞数
文章标签： 数据结构与算法
原文链接：http://www.cnblogs.com/ywwyww/p/8511071.html
版权

题解

分治FFT

　　设\(f_i\)为\(i\)个点组成的无向图个数，\(g_i\)为\(i\)个点组成的无向连通图个数

　　经过简单的推导（枚举\(1\)所在的连通块大小），有：
\[ f_i=2^{\frac{i(i-1)}{2}} \]

\[ \begin{align} g_i&=f_i-\sum_{j=1}^{i-1}\binom{n-1}{j-1}g_jf_{i-j}\\ &=f_i-(i-1)!\sum_{j=1}^{i-1}\frac{g_j}{(j-1)!}\frac{f_{i-j}}{(i-j)!} \end{align} \]

　　用CDQ分治+FFT优化。

　　就是每次先求出\(g_l\cdots g_{mid+1}\)，然后卷上\(f_1\cdots f_{len}\)，加到\(g_{mid+1}\cdots g_r\)上面去。

　　时间复杂度：\(O(n\log^2 n)\)

多项式求逆

\[ \begin{align} f_i&=\sum_{j=1}^n\binom{i-1}{j-1}g_jf_{i-j}\\ \frac{f_i}{(i-1)!}&=\sum_{j=1}^n\frac{g_j}{(j-1)!}\frac{f_{i-j}}{(i-j)!} \end{align} \]

　　设
\[ \begin{align} A&=\sum_{i\geq 1}\frac{f_i}{(i-1)!}x^i\\ B&=\sum_{i\geq 1}\frac{g_i}{(i-1)!}x^i\\ C&=\sum_{i\geq 0}\frac{f_i}{i!}x^i \end{align} \]
　　所以
\[ \begin{align} A&=B\times C~~~~~~(mod~x^{n+1})\\ B&=A\times C^{-1}~~(mod~x^{n+1}) \end{align} \]
　　时间复杂度：\(O(n\log n)\)

多项式求ln

　　设
\[ \begin{align} G(x)&=\sum_{i\geq 0}\frac{2^{\binom{i}{2}}}{i!}x^i\\ F(x)&=\sum_{i\geq 0}\frac{f_i}{i!}x^i \end{align} \]

　　根据指数生成函数和有标号计数的那套理论，
　　由\(0\)个连通块组成的连通图的个数为\(\frac{{F(x)}^0}{0!}\)；
　　由\(1\)个连通块组成的连通图的个数为\(\frac{{F(x)}^1}{1!}\)；
　　由\(2\)个连通块组成的连通图的个数为\(\frac{{F(x)}^2}{2!}\)；
　　\(\vdots\)

　　这些加起来就是无向图的个数\(G(x)\)

　　所以
\[ \begin{align} G(x)&=\frac{{F(x)}^0}{0!}+\frac{{F(x)}^1}{1!}+\frac{{F(x)}^2}{2!}+\frac{{F(x)}^3}{3!}+\cdots\\ &=\sum_{i\geq 0}\frac{{F(x)}^i}{i!}\\ &=e^{F(x)}\\ F(x)&=\ln(G(x)) \end{align} \]

　　直接求ln即可。

　　时间复杂度：\(O(n\log n)\)

代码

分治FFT

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
ll p=1004535809;
ll fp(ll a,ll b)
{
    ll s=1;
    while(b)
    {
        if(b&1)
            s=s*a%p;
        a=a*a%p;
        b>>=1;
    }
    return s;
}
ll inv(ll x)
{
    return fp(x,p-2);
}
namespace ntt
{
    ll w1[300010];
    ll w2[300010];
    int rev[300010];
    int n;
    void init(int m)
    {
        n=1;
        while(n<=m)
            n<<=1;
        int i;
        for(i=2;i<=n;i<<=1)
        {
            w1[i]=fp(3,(p-1)/i);
            w2[i]=inv(w1[i]);
        }
        rev[0]=0;
        for(i=1;i<n;i++)
            rev[i]=(rev[i>>1]>>1)|(i&1?n>>1:0);
    }
    void ntt(ll *a,int t)
    {
        int i,j,k;
        ll u,v,w,wn;
        for(i=0;i<n;i++)
            if(rev[i]<i)
                swap(a[i],a[rev[i]]);
        for(i=2;i<=n;i<<=1)
        {
            wn=(t==1?w1[i]:w2[i]);
            for(j=0;j<n;j+=i)
            {
                w=1;
                for(k=j;k<j+i/2;k++)
                {
                    u=a[k];
                    v=a[k+i/2]*w%p;
                    a[k]=(u+v)%p;
                    a[k+i/2]=(u-v+p)%p;
                    w=w*wn%p;
                }
            }
        }
        if(t==-1)
        {
            u=inv(n);
            for(i=0;i<n;i++)
                a[i]=a[i]*u%p;
        }
    }
};
ll a[300010];
ll b[300010];
ll fac[300010];
ll invfac[300010];
ll in[300010];
ll f[300010];
ll g[300010];
void solve(int l,int r)
{
    if(l==r)
    {
        g[l]=(f[l]-fac[l-1]*g[l]%p+p)%p;
        return;
    }
    int mid=(l+r)>>1;
    solve(l,mid);
    int len=r-l+1;
    ntt::init(len);
    int i;
    a[0]=b[0]=0;
    if(l==3&&r==4)
        int x=1;
    if(l==1&&r==4)
        int x=1;
    for(i=1;i<=mid-l+1;i++)
        a[i]=g[i+l-1]*invfac[i+l-1-1]%p;
    for(i=mid-l+1+1;i<ntt::n;i++)
        a[i]=0;
    for(i=1;i<=len;i++)
        b[i]=f[i]*invfac[i]%p;
    for(i=len+1;i<ntt::n;i++)
        b[i]=0;
    ntt::ntt(a,1);
    ntt::ntt(b,1);
    for(i=0;i<ntt::n;i++)
        a[i]=a[i]*b[i]%p;
    ntt::ntt(a,-1);
    for(i=mid+1;i<=r;i++)
        g[i]=(g[i]+a[i-l+1])%p;
    solve(mid+1,r);
//  fprintf(stderr,"%d %d\n",l,r);
}
int main()
{
//  freopen("bzoj3456.in","r",stdin);
//  freopen("bzoj3456.out","w",stdout);
    int n;
    scanf("%d",&n);
    invfac[0]=fac[0]=1;
    int i;
    in[0]=in[1]=1;
    for(i=2;i<=n;i++)
        in[i]=(-(p/i)*in[p%i]%p+p)%p;
    for(i=1;i<=n;i++)
    {
        fac[i]=fac[i-1]*i%p;
        invfac[i]=invfac[i-1]*in[i]%p;
    }
    for(i=1;i<=n;i++)
        f[i]=fp(2,(ll(i)*(i-1)/2)%(p-1));
    memset(g,0,sizeof g);
    solve(1,n);
    printf("%lld\n",g[n]);
    return 0;
}

多项式求逆

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
ll p=1004535809;
ll g=3;
ll fp(ll a,ll b)
{
    ll s=1;
    while(b)
    {
        if(b&1)
            s=s*a%p;
        a=a*a%p;
        b>>=1;
    }
    return s;
}
namespace ntt
{
    ll w1[300010];
    ll w2[300010];
    int rev[300010];
    int n;
    void init(int m)
    {
        n=1;
        while(n<m)
            n<<=1;
        int i;
        for(i=2;i<=n;i++)
        {
            w1[i]=fp(g,(p-1)/i);
            w2[i]=fp(w1[i],p-2);
        }
        rev[0]=0;
        for(i=1;i<n;i++)
            rev[i]=(rev[i>>1]>>1)|(i&1?n>>1:0);
    }
    void ntt(ll *a,int t)
    {
        int i,j,k;
        ll u,v,w,wn;
        for(i=0;i<n;i++)
            if(rev[i]<i)
                swap(a[i],a[rev[i]]);
        for(i=2;i<=n;i<<=1)
        {
            wn=(t==1?w1[i]:w2[i]);
            for(j=0;j<n;j+=i)
            {
                w=1;
                for(k=j;k<j+i/2;k++)
                {
                    u=a[k];
                    v=a[k+i/2]*w%p;
                    a[k]=(u+v)%p;
                    a[k+i/2]=(u-v+p)%p;
                    w=w*wn%p;
                }
            }
        }
        if(t==-1)
        {
            u=fp(n,p-2);    
            for(i=0;i<n;i++)
                a[i]=a[i]*u%p;
        }
    }
    ll x[300010];
    ll y[300010];
    void copy_clear(ll *a,ll *b,int m)
    {
        int i;
        for(i=0;i<m;i++)
            a[i]=b[i];
        for(i=m;i<n;i++)
            a[i]=0;
    }
    void copy(ll *a,ll *b,int m)
    {
        int i;
        for(i=0;i<m;i++)
            a[i]=b[i];
    }
    void inverse(ll *a,ll *b,int m)
    {
        if(m==1)
        {
            b[0]=fp(a[0],p-2);
            return;
        }
        inverse(a,b,m>>1);
        init(2*m);
        copy_clear(x,a,m);
        copy_clear(y,b,m>>1);
        ntt(x,1);
        ntt(y,1);
        int i;
        for(i=0;i<n;i++)
            x[i]=(2*y[i]%p-x[i]*y[i]%p*y[i]%p+p)%p;
        ntt(x,-1);
        copy(b,x,m);
    }
};
ll a[300010];
ll b[300010];
ll c[300010];
ll fac[300010];
int main()
{
//  freopen("bzoj3456.in","r",stdin);
    int n;
    scanf("%d",&n);
    int i;
    fac[0]=1;
    int m=1;
    while(m<=2*n)
        m<<=1;
    for(i=1;i<=n;i++)
        fac[i]=fac[i-1]*i%p;
    a[0]=0;
    for(i=1;i<=n;i++)
        a[i]=fp(2,(ll(i-1)*i/2)%(p-1))*fp(fac[i-1],p-2)%p;
    b[0]=1;
    for(i=1;i<=n;i++)
        b[i]=fp(2,(ll(i-1)*i/2)%(p-1))*fp(fac[i],p-2)%p;
    ntt::inverse(b,c,m>>1);
    for(i=n+1;i<m;i++)
        c[i]=0;
    ntt::init(m);
    ntt::ntt(a,1);
    ntt::ntt(c,1);
    for(i=0;i<m;i++)
        a[i]=a[i]*c[i]%p;
    ntt::ntt(a,-1);
    printf("%d\n",a[n]*fac[n-1]%p);
    return 0;
}

多项式ln

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
#include<cmath>
#include<functional>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
void sort(int &a,int &b)
{
    if(a>b)
        swap(a,b);
}
void open(const char *s)
{
#ifndef ONLINE_JUDGE
    char str[100];
    sprintf(str,"%s.in",s);
    freopen(str,"r",stdin);
    sprintf(str,"%s.out",s);
    freopen(str,"w",stdout);
#endif
}
int rd()
{
    int s=0,c;
    while((c=getchar())<'0'||c>'9');
    do
    {
        s=s*10+c-'0';
    }
    while((c=getchar())>='0'&&c<='9');
    return s;
}
int upmin(int &a,int b)
{
    if(b<a)
    {
        a=b;
        return 1;
    }
    return 0;
}
int upmax(int &a,int b)
{
    if(b>a)
    {
        a=b;
        return 1;
    }
    return 0;
}
const ll p=1004535809;
const ll g=3;
ll fp(ll a,ll b)
{
    ll s=1;
    while(b)
    {
        if(b&1)
            s=s*a%p;
        a=a*a%p;
        b>>=1;
    }
    return s;
}
const int maxn=600000;
ll inv[maxn];
namespace ntt
{
    ll w1[maxn];
    ll w2[maxn];
    int rev[maxn];
    int n;
    void init(int m)
    {
        n=1;
        while(n<m)
            n<<=1;
        int i;
        for(i=2;i<=n;i++)
        {
            w1[i]=fp(g,(p-1)/i);
            w2[i]=fp(w1[i],p-2);
        }
        rev[0]=0;
        for(i=1;i<n;i++)
            rev[i]=(rev[i>>1]>>1)|(i&1?n>>1:0);
    }
    void ntt(ll *a,int t)
    {
        int i,j,k;
        ll u,v,w,wn;
        for(i=0;i<n;i++)
            if(rev[i]<i)
                swap(a[i],a[rev[i]]);
        for(i=2;i<=n;i<<=1)
        {
            wn=(t==1?w1[i]:w2[i]);
            for(j=0;j<n;j+=i)
            {
                w=1;
                for(k=j;k<j+i/2;k++)
                {
                    u=a[k];
                    v=a[k+i/2]*w%p;
                    a[k]=(u+v)%p;
                    a[k+i/2]=(u-v)%p;
                    w=w*wn%p;
                }
            }
        }
        if(t==-1)
        {
            u=fp(n,p-2);    
            for(i=0;i<n;i++)
                a[i]=a[i]*u%p;
        }
    }
    ll x[maxn];
    ll y[maxn];
    void copy_clear(ll *a,ll *b,int m)
    {
        int i;
        for(i=0;i<m;i++)
            a[i]=b[i];
        for(i=m;i<n;i++)
            a[i]=0;
    }
    void copy(ll *a,ll *b,int m)
    {
        int i;
        for(i=0;i<m;i++)
            a[i]=b[i];
    }
    void inverse(ll *a,ll *b,int m)
    {
        if(m==1)
        {
            b[0]=fp(a[0],p-2);
            return;
        }
        inverse(a,b,m>>1);
        init(m<<1);
        copy_clear(x,a,m);
        copy_clear(y,b,m>>1);
        ntt(x,1);
        ntt(y,1);
        int i;
        for(i=0;i<n;i++)
            x[i]=(2*y[i]-x[i]*y[i]%p*y[i])%p;
        ntt(x,-1);
        copy(b,x,m);
    }
    void integrate(ll *a,ll *b,int m)
    {
        int i;
        for(i=0;i<m-1;i++)
            b[i]=(i+1)*a[i+1]%p;
        b[m-1]=0;
    }
    void differential(ll *a,ll *b,int m)
    {
        int i;
        for(i=m-1;i>=1;i--)
            b[i]=a[i-1]*inv[i]%p;
        b[0]=0;
    }
    void ln(ll *a,ll *b,int m)
    {
        static ll c[maxn],d[maxn],e[maxn];
        integrate(a,c,m);
        inverse(a,d,m);
        init(m<<1);
        ntt(c,1);
        ntt(d,1);
        int i;
        for(i=0;i<n;i++)
            e[i]=c[i]*d[i];
        ntt(e,-1);
        differential(e,b,m);
    }
    void exp(ll *a,ll *b,int m)
    {
        if(m==1)
        {
            b[0]=1;
            return;
        }
        exp(a,b,m>>1);
        int i;
        for(i=m>>1;i<m;i++)
            b[i]=0;
        ln(b,y,m);
        init(m<<1);
        copy_clear(x,a,m);
        x[0]++;
        for(i=0;i<m;i++)
            x[i]=(x[i]-y[i])%p;
        copy_clear(y,b,m);
        ntt(x,1);
        ntt(y,1);
        for(i=0;i<n;i++)
            x[i]=x[i]*y[i]%p;
        ntt(x,-1);
        for(i=0;i<m;i++)
            b[i]=x[i];
    }
};
ll a[maxn];
ll b[maxn];
ll fac[maxn];
ll invfac[maxn];
int main()
{
    open("bzoj3456");
    int i,n;
    scanf("%d",&n);
    int m=1;
    while(m<=n+1)
        m<<=1;
    inv[0]=inv[1]=fac[0]=fac[1]=invfac[0]=invfac[1]=1;
    for(i=2;i<m;i++)
    {
        inv[i]=-(p/i)*inv[p%i]%p;
        fac[i]=fac[i-1]*i%p;
        invfac[i]=invfac[i-1]*inv[i]%p;
    }
    for(i=0;i<m;i++)
        a[i]=fp(2,(ll(i)*(i-1)/2)%(p-1))*invfac[i]%p;
    ntt::ln(a,b,m);
    ll ans=(b[n]*fac[n]%p+p)%p;
    printf("%lld\n",ans);
    return 0;
}

转载于:https://www.cnblogs.com/ywwyww/p/8511071.html

weixin_30745641
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
【XSY1332】【BZOJ3456】轩辕朗的城市规划 无向连通图计数 CDQ分治 FFT 多项式求逆 多项式ln
ez_yww的博客
09-05 621
题解　　设fif_i为nn个点组成的无向图个数，gig_i为nn个点组成的无向连通图个数　　经过简单的推导（枚举11所在的连通块大小），有： fi=2n(n−1)2 f_i=2^{\frac{n(n-1)}{2}} gi=fi−∑j=1i−1(n−1j−1)gjfi−j=fi−(i−1)!∑j=1i−1gj(j−1)!fi−j(i−j)! \begin{align} g_i&=f_i-\sum_{
【XSY1538】连在一起的幻想乡 数学 无向连通图计数
ez_yww的博客
08-15 1133
题目大意​　　给你n,pn,p，求nn个点组成的所有无向连通图图的边数的平方和​　　n≤2000,m≤109n\leq 2000,m\leq {10}^9题解​　　设m=n(n−1)2,h0n=nm=\frac{n(n-1)}{2},h0_n=n个点无向图的个数，h1n=nh1_n=n个点组成的所有无向图的边数之和，h2n=nh2_n=n个点组成的所有无向图的边数的平方和，f0n=nf0_n=n个点
【BZOJ3456】轩辕朗的城市规划 EGF+多项式求ln
diaoyoutun2652的博客
05-22 101
我们构造$f(i)$和$g(i)$。 其中$f(x)$表示由$x$个节点构成的无向简单连通图的个数。 $g(x)$表示有$x$个节点构成的无向简单图（不要求连通）的个数。 显然，由$x$个节点构成的无向简单图最多能有$\binom{x}{2}$条边，那么$g(x)=2^{\binom{x}{2}}$。 然后我们构造$f(x)$和$g(x)$的$EGF$： $F(x)=\sum_...
[n点无向图个数 分治FFT || 多项式求逆] BZOJ 3456 城市规划
雯舞
09-14 811
然后就分治FFT了 #include #include #include using namespace std; typedef long long ll; inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin)
题解 有标号无向连通图计数
im_zyINF的博客
06-01 822
计数dp
bzoj 3730. 震波（动态点分治 + vector树状数组）
01-06
（自己写的动态点分治巨垃圾，常数是别人的两倍） 用动态开点线段树死活过不去，学了一波大佬用 vector 开树状数组立马就卡过去了 考虑点分树的做法，在点分树上每个点以距离为下标建一棵线段树，每次询问查询子树的...
Node.js-ZOJCH是BZOJ题库的离线版
08-09
ZOJCH是BZOJ题库的离线版
CTSC 2011 无穷图的桥(BZOJ 2307) 题解.ppt
03-17
CTSC 2011 无穷图的桥(BZOJ 2307) 题解.ppt
wangjunrui666#bzoj-problem##5218. [Lydsy2017省队十连测]友好城市1
07-25
对于两座不同的城市a 和b，如果a 能通过这些单向道路直接或间接到达b，且b 也能如此到达a，那么它们就会被认为是一对友好城市。eland 的交通系统十分特殊，
wangjunrui666#bzoj-problem##3873. [Ahoi2014&Jsoi2014]拼图1
07-25
作为一个计算机科学家，JYY有一套黑白色的拼图，他希望通过合理的拼接，使得拼出的最终图案中，能包含面积最大的全白色子矩形。之后JYY发现，可以通过改变这S块拼图
【XSY1295】calc $n$个点$n$条边无向连通图计数 prufer序列
ez_yww的博客
09-19 1305
题目大意　　求nn个点nn条边的无向连通图的个数　　n≤5000n\leq 5000题解　　显然是一个环上有很多外向树。　　首先有一个东西：nn个点选kk个点作为树的根的生成森林个数为： (nk)×nn−k−1×k \binom{n}{k}\times n^{n-k-1}\times k 　　前面(nk)\binom{n}{k}是这些根的选编号的方案树，后面是prufer序列得到的：前面n−k
【XSY1295】calc n个点n条边无向连通图计数 prufer序列
weixin_30276935的博客
03-05 754
题目大意 　　求\(n\)个点\(n\)条边的无向连通图的个数 　　\(n\leq 5000\) 题解 　　显然是一个环上有很多外向树。 　　首先有一个东西：\(n\)个点选\(k\)个点作为树的根的生成森林个数为： \[ \binom{n}{k}\times n^{n-k-1}\times k \] 　　前面\(\binom{n}{k}\)是这些根的选编号的方案数，后面是prufer序列得到的：...
BZOJ1002 无向联通图的生成树计数
越努力越幸运—liupu
03-16 1029
分析：求无向连通图的生成树个数我们要用到Matrix-Tree定理，这里我只给出具体实现方法，至于怎么证明。。。。我现在也没耐心去看。。有兴趣的同学可以看看 周冬老师的《生成树的计数及其应用》。 好的，现在给出实现方法：点击打开链接，哈哈自己看吧，我是不太想写上来了。 按照算法将行列式构造出来之后我们可以发现这个行列式是有特点的，可以用行列式的展开来得到一个递推关系式（所以说线代还是要好好学滴
bzoj 3456 城市规划 - 图计数 - NTT - 多项式求逆 - NTT学习笔记 - 多项式求逆学习笔记
Mys_C_K的博客
04-12 376
设f[x]表示x个点的连通图，h[x]=2^{C(n,2)}为任意图的数量。 Hn=∑ni=1(n−1i−1)Fi×&amp;amp;nbsp;Hn−iHn=∑i=1n(n−1i−1)Fi×&amp;amp;nbsp;Hn−iH_n=\sum_{i=1}^n \binom{n-1}{i-1}F_i\times\ H_{n-i} Hn(n−1)!=∑ni=1Fi(i−1)!×&amp;amp;nbsp;Hn−i(n−i)!Hn(n−1)!=∑...
指数型生成函数（二）：贝尔数、无向连通图计数、弱连通DAG计数
hans774882968的博客
08-09 412
贝尔数：n个元素的集合的划分的个数 注意到S2的定义就是n个元素放入k个无区别的盒子的个数，“集合”是一种无区别的盒子，所以贝尔数等于一行S2的和。 定义dp[n]为所求，则考虑1号元素所在集合的元素个数： dp[n]=∑i=1ndp[n−i]∗Cn−1i−1=∑i=0n−1dp[n−1−i]∗Cn−1i,dp[0]=1 dp[n]=\sum_{i=1}^ndp[n-i]*C_{n-1}^{i-1}=\sum_{i=0}^{n-1}dp[n-1-i]*C_{n-1}^i,dp[0]=1 dp[n]=i=.
带标号的连通图计数
qichi_bj的专栏
12-10 2241
统计 n 个顶点的连通图有多少个，每个顶点有标号。记 f(n) 为带标记的连通图的个数，g(n) 为带标记的非连通图的个数，则 f(n) + g(n) = h(n) = 2^(n(n-1)/2) ；考虑 g(n) 的计数：按照结点 1 所在的联通分量的点数分类，设结点 1 所在联通分量包含 k 点（k  = 1 ... n-1）则有 C(n-1, k -1) 种不同情况，每种情况下，结点 1 所在
[带标号无向连通图计数 容斥原理 多项式求逆 多项式求ln 模板题] BZOJ 3456 城市规划
雯舞
02-21 2599
可以通过容斥求出答案的表达式fi=2C2i−∑j=1i−1Cj−1i−1∗fj∗2C2i−jf_i=2^{C_i^2}-\sum_{j=1}^{i-1} C_{i-1}^{j-1}*f_j*2^{C_{i-j}^2} 其中前一部分表示i个点任意连边 后半部分枚举1所在的连通块然后容斥掉∑j=1ifj(j−1)!∗2C2i−j(i−j)!=2C2i(i−1)!\sum_{j=1}^i {f_j \o
【Tinsen】A1493. 城市规划【简单无向连通图个数——NTT】
poursoul
05-15 2960
传送门：【Tinsen】A1493. 城市规划前言：今天很不在状态= =，NTT里面竟然同一个地方写错两次……而且脑补的时候明显比平时迟钝了很多……题目分析：首先我们可以轻易知道N个点能构成的无向图的个数：2C2n2^{C_n^2}，令其等于d(N)d(N)。然后假设我们固定一个点，枚举这个点所在连通块的大小，这样我们得到N个点的无向连通图公式： dp[n]=d(n)−∑n−2i=0Cin−1dp
BZOJ3456【分治FFT】
Lethelody的专栏
07-07 3488
然而会TLE. /* I will wait for you */ #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define make(a,b) make_p
BZOJ2179: FFT快速傅立叶 FFT实现高精度乘法
最新发布
05-14
这道题目是要求用 FFT 实现高精度乘法。我们可以先将两个高精度数转换为多项式，然后对这两个多项式进行相乘，最后将结果转换为高精度数输出即可。 具体来说，将两个高精度数 $a$ 和 $b$ 分别转换为多项式 $A(x)$ ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值