【XSY3042】石像 拓扑排序 状压DP 洲阁筛

最新推荐文章于 2019-09-05 16:26:37 发布
最新推荐文章于 2019-09-05 16:26:37 发布
阅读量103 收藏 1
点赞数
文章标签： 数据结构与算法
原文链接：http://www.cnblogs.com/ywwyww/p/9229200.html
版权

题目大意

　　有 \(n\) 个整数 \(a_1,a_2,\ldots,a_n\)，每个数的范围是 \([1,m]\)。还有 \(k\) 个限制，每个限制 \(x_i,y_i\) 表示 \(a_{x_i}\leq a_{y_i}\)。

　　问有多少种不同的情况，以及所有情况中 \({\sigma_0(\gcd(a_1,a_2,\ldots,a_n))}^3\) 的和。

　　\(n\leq 20,m\leq {10}^{10}\)

题解

　　记 \(f(x)\) 为当 \(m=x\) 时第一问的答案。

　　记 \(g(x)\) 为当 \(a_1,a_2,\ldots,a_n\) 有 \(x\) 种不同的值时第一问的答案。

　　\(g(x)\) 要怎么求呢？

　　tarjan 求强连通分量，然后缩点，状压DP。

　　直接枚举子集是 \(O(n3^n)\) 的。

　　因为缩点后是一个 DAG，所以按拓扑序来转移就好了。

　　设 \(h_{i,j,S}\) 为 当前考虑的值为 \(i\)，正在考虑 \(a_j\) 这个数，已经确定数的集合为 \(S\) 的方案数。

　　每次决定 \(a_j\) 这个数的值是不是 \(i\) 就可以了。

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

　　\(f(x)\) 要怎么求呢？

　　注意到 \(a_1,a_2,\ldots,a_n\) 最多只有 \(n\) 种不同的值，所以可以用 \(g(1),g(2),\ldots,g(n)\) 计算出 \(f(x)\)：
\[ f(x)=\sum_{i=1}^n\binom{x}{i}g(i) \]
　　求组合数时要记下 \(2\) 的次数和非 \(2\) 的倍数的乘积两部分。

　　或者你求出前 \(O(n)\) 项然后拉格朗日插值也是可以的。

　　这样就求出第一问的答案了。

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

　　那第二问要怎么求呢？
\[ \begin{align} &\sum_{a_1=1}^m\sum_{a_2=1}^m\cdots\sum_{a_n=1}^m\sigma_0({\gcd(a_1,a_2,\ldots,a_n)}^3)\\ =&\sum_{i=1}^m\sigma_0(i^3)^3\sum_{a_1=1}^m\sum_{a_2=1}^m\cdots\sum_{a_n=1}^m[\gcd(a_1,a_2,\ldots,a_n)=i]\\ =&\sum_{i=1}^m(\sum_{j\mid i}\sigma_0(j^3)^3\mu(\frac{i}{j}))\sum_{a_1=1}^m\sum_{a_2=1}^m\cdots\sum_{a_n=1}^m[i\mid\gcd(a_1,a_2,\ldots,a_n)]\\ =&\sum_{i=1}^m(\sum_{j\mid i}\sigma_0(j^3)^3\mu(\frac{i}{j}))f(\lfloor\frac{m}{i}\rfloor) \end{align} \]
　　那么前面那部分要怎么求呢？

　　容易发现 \(h(n)=\sum_{d\mid n}\sigma_0(d^3)^3\mu(\frac{n}{d})\) 是一个积性函数，满足
\[ \begin{align} h(1)&=1\\ h(p)&=1\times (-1)+4^3\times 1=63\\ h(p^c)&={(3(c-1)+1)}^3(-1)+{(3c+1)}^3(1)\\ &=(27c^3+27c^2+9c+1)-(27c^2-54c^2+36c-8)\\ &=81c^2-27c+9 \end{align} \]
　　所以可以大力洲阁筛/min_25筛求。

　　可能需要卡一下常。

　　时间复杂度：\(O(n^22^n+n\sqrt m+\frac{m^\frac{3}{4}}{\log m})\)

代码

#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,b=0;
    while(((c=getchar())<'0'||c>'9')&&c!='-');
    if(c=='-')
    {
        c=getchar();
        b=1;
    }
    do
    {
        s=s*10+c-'0';
    }
    while((c=getchar())>='0'&&c<='9');
    return b?-s:s;
}
void put(int x)
{
    if(!x)
    {
        putchar('0');
        return;
    }
    static int c[20];
    int t=0;
    while(x)
    {
        c[++t]=x%10;
        x/=10;
    }
    while(t)
        putchar(c[t--]+'0');
}
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;
}

int fp(int a,ll b)
{
    int s=1;
    for(;b;b>>=1,a*=a)
        if(b&1)
            s*=a;
    return s;
}

namespace gao1
{
    const int N=30;
    int n,m;
    vector<int> g[N],g2[N];
    int lx[N*N],ly[N*N];
    int dfn[N],low[N],b[N],st[N],top,c[N],ti,ccnt;
    int f[1<<20],a[N];
    int s[N];
    int d[N][N];
    int inv[N],pw[N];
    void dfs(int x)
    {
        dfn[x]=low[x]=++ti;
        st[++top]=x;
        b[x]=1;
        for(auto v:g[x])
            if(b[v]!=2)
            {
                if(!b[v])
                    dfs(v);
                low[x]=min(low[x],low[v]);
            }
        if(low[x]>=dfn[x])
        {
            ccnt++;
            int v;
            do
            {
                v=st[top--];
                b[v]=2;
                c[v]=ccnt;
            }
            while(v!=x);
        }
    }
    void add(int &a,int b)
    {
        a+=b;
    }
    void gao()
    {
        scanf("%d",&m);
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&lx[i],&ly[i]);
            g[lx[i]].push_back(ly[i]);
        }
        for(int i=1;i<=n;i++)
            if(!b[i])
                dfs(i);
        for(int i=1;i<=m;i++)
            if(c[lx[i]]!=c[ly[i]])
            {
                g2[c[lx[i]]].push_back(c[ly[i]]);
                a[c[lx[i]]]|=1<<(c[ly[i]]-1);
            }
        n=ccnt;
        f[0]=1;
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
                for(int k=0;k<1<<n;k++)
                    if(!((k>>(j-1))&1)&&(k&a[j])==a[j])
                        add(f[k|(1<<(j-1))],f[k]);
            s[i]=f[(1<<n)-1];
        }
        for(int i=0;i<=n;i++)
        {
            d[i][0]=1;
            for(int j=1;j<=i;j++)
                d[i][j]=d[i-1][j-1]+d[i-1][j];
        }
        for(int i=n;i>=1;i--)
        {
            for(int j=i-1;j>=1;j--)
                s[i]+=((i-j)&1?-1:1)*d[i][j]*s[j];
        }
        for(int i=1;i<=n;i++)
        {
            int x=i;
            while(!(x&1))
            {
                pw[i]++;
                x>>=1;
            }
            inv[i]=fp(x,(1ll<<31)-1);
        }
    }
    int query(ll x)
    {
        int res=0;
        int prod=1;
        int num=0;
        for(int i=1;i<=n&&i<=x;i++)
        {
            ll y=x-i+1;
            while(num<=20&&!(y&1))
            {
                y>>=1;
                num++;
            }
            num-=pw[i];
            prod*=inv[i];
            prod*=y;
            res+=(prod<<num)*s[i];
        }
        return res;
    }
}

namespace gao2
{
    const int M=100010;
    const int m=100000;
    ll n;
    int pri[M],b[M],cnt;
    int f1[M],f2[M],s1[M],s2[M];
    void init()
    {
        for(int i=2;i<=m;i++)
        {
            if(!b[i])
                pri[++cnt]=i;
            for(int j=1;j<=cnt&&i*pri[j]<=m;j++)
            {
                b[i*pri[j]]=1;
                if(i%pri[j]==0)
                    break;
            }
        }
    }
    void gao()
    {
        init();
        int mx=n/(m+1);
        for(int i=2;i<=m;i++)
            f1[i]=i-1;
        for(int i=1;i<=mx;i++)
            f2[i]=n/i-1;
        for(int i=1;i<=cnt;i++)
        {
            int x=pri[i];
            int y=f1[x-1];
            int mx1=min(n/x/(m+1),n/x/x);
            int mx2=min((ll)mx,n/x/x);
            ll p2=(ll)pri[i]*pri[i];
            ll o=n/x;
            for(int j=1;j<=mx1;j++)
                f2[j]-=f2[j*x]-y;
            for(int j=mx1+1;j<=mx2;j++)
                f2[j]-=f1[o/j]-y;
            for(int j=m;j>=p2;j--)
                f1[j]-=f1[j/x]-y;
        }
        for(int i=2;i<=m;i++)
        {
            f1[i]*=63;
            s1[i]=f1[i];
        }
        for(int i=1;i<=mx;i++)
        {
            f2[i]*=63;
            s2[i]=f2[i];
        }
        for(int i=cnt;i>=1;i--)
        {
            int x=pri[i];
            int mx1=min((ll)mx,n/x/x);
            for(int j=1;j<=mx1;j++)
            {
                s2[j]-=63;
                ll x2=n/j/x;
                for(int k=1;x2;x2/=x,k++)
                    s2[j]+=((x2<=m?s1[x2]-f1[min((int)x2,x)]:s2[n/x2]-f1[x])+1)*(81*k*k-27*k+9);
            }
            ll mn=(ll)pri[i]*pri[i];
            for(int j=m;j>=mn;j--)
            {
                s1[j]-=63;
                int x2=j/x;
                for(int k=1;x2;x2/=x,k++)
                    s1[j]+=(s1[x2]-f1[min(x2,x)]+1)*(81*k*k-27*k+9);
            }
        }
    }
    int query(ll x)
    {
        if(x<=0)
            return 0;
        return (x<=m?s1[x]:s2[n/x])+1;
    }
}

int main()
{
    open("statue");
    scanf("%d%lld",&gao1::n,&gao2::n);
    gao1::gao();
    gao2::gao();
    int ans=0;
    ll n=gao2::n;
    for(ll i=1,j;i<=n;i=j+1)
    {
        j=n/(n/i);
        ans+=(gao2::query(j)-gao2::query(i-1))*gao1::query(n/i);
    }
    printf("%u %u\n",gao1::query(n),ans);
    return 0;
}

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

weixin_30451709
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
OFDM.rar_XSY_ofdm qpsk awgn_ofdm+qpsk_两径_两径信道
09-21
OFDM-QPSK系统在不同信道（AWGN、AWGN+两径衰落信道、AWGN+三径衰落信道）下的接收端误码性能
洲阁筛/Min25筛
Freopen的博客
12-16 1216
1：完整的Min25筛学习方案 2：虽然不是但我也挂
洲阁筛学习总结
热门推荐
SemiWaker的博客
06-27 1万+
洲阁筛学习总结By Semiwaker 积性函数求和
洲阁筛详解
a2850363633的博客
03-01 540
你还真信了 丢链接 这筛对积性函数的要求不同于杜教筛，只消函数在自变量为质数或质数整数幂时是一个低阶多项式即可。以下n<=1e11。 首先有一个性质：1~n的每个数，大于$\sqrt{n}$的质因子只有一个。根据是否有大于$\sqrt{n}$的质因子，再根据他是积性函数，得 $\sum_{i=1}^nf(i)=\sum_{i=1}^{n}[i \ \ has \ \ no ...
【Learning】积性函数前缀和——洲阁筛（min_25写法）
weixin_30731287的博客
06-14 171
问题描述 　　 　　洲阁筛解决的问题主要是\(n\)范围较大的积性函数前缀和。 　　 ​　　已知一积性函数\(f(i)\)，求\(\sum_{i=1}^nf(i)\)。 　　 　　\(n\leq10^{12}\). 　　 　　 　　 求解方法 　　 　　如果\(f(i)\)在质数处的取值比较简单，那么可以运用洲阁筛来求解。 　　 ​　　我们需要两个辅助数组。 　　 \(g_{i,j}\) 　　 　...
洲阁筛 学习笔记
weixin_30684743的博客
06-12 189
洲阁筛 对具有如下性质的积性函数 \(F\) \(p\in prime\) 时，\(F(p)=G(p)\) $p\in prime $ 时，\(F(p^c)=T(p^c),c>1\)(如无特殊说明，和式里面的 \(p\) 同样指代质数) 其中，\(G,T\) 为关于 \(p,c\) 的低阶多项式或幂函数 洲阁筛的作用是在 \(O(\frac{n^{\frac{3}{4}}}{\log n...
三菱FX PLC 定位控制三菱伺服简要总结XSY.pdf
06-04
三菱FX PLC 定位控制三菱伺服简要总结XSY.pdf 图文 含程序截图 可编程控制器： FX3GA-60MT-CM（三菱） 伺服马达： HG-KR13J（三菱） 伺服驱动器： MR-J4-10A（三菱） 伺服参数设置软件：MR Configurator2（三菱） ...
Xsy-Light-KO.zip_单片机开发_C++_
08-12
单片机用十字路口灯控制，绿灯40s，后3s闪亮，红灯30s，可修改
安川机器人FS100相关总结20170727XSY.pdf
08-16
安川机器人FS100相关总结20170727XSY.pdf
FANUC机器人与西门子TIA Portal V15之间PROFINET通讯XSY.pdf
06-04
一、 设置机器人侧 ② 进入PROFINET通讯页面，选择2频道（1频道为机器人做主站时用，2频道为机器人做从站用），将【1频道】禁用，启用【2频道】（显示“有效”表示当前为禁用状态，显示“无效”即表示当前为启用状态...
【XSY2691】中关村 卢卡斯定理 数位DP
ez_yww的博客
01-07 363
题目描述 　　在一个kk维空间中，每个整点被黑白染色。对于一个坐标为(x1,x2,…,xk)(x_1,x_2,\ldots,x_k)的点，他的颜色我们通过如下方式计算： 如果存在一维坐标是00，则颜色是黑色。 如果这个点是(1,1,…,1)(1,1,\ldots,1)（每一维都是11），这个点的颜色是白色 如果这个点的kk个前驱（任取一维坐标减11）中的白点有奇数个，那么这个点的颜色就是白色
杜教筛 狄利克雷卷积入门
Stydwn的博客
09-05 478
补下暑假的坑 杜教筛看了两天勉勉强强才看懂了一天 杜教牛逼 杜教筛意义 其实，对于一般的数论题，线性筛已经非常的优秀了。 但是就是有那些duliu出题人，硬是要把数据出到1e10之类的，就看你会不会杜教筛，min_25筛，洲阁筛等各种神奇的筛法。(PS:后面这两个筛法我是真的不会了QAQ) 了解杜教筛之前必须得了解一下狄利克雷卷积 听名字，似乎是一个很高深的东西，其实就是一个公式 证明莫...
min_25筛入门
某司机的黑车
07-11 1148
前言 我好菜啊现在才会 好像踩爆了洲鸽（划掉）阁筛？ （貌似min_25筛就是洲阁筛的另一种写法） 如果有错欢迎dalao指出因为我特么在学的时候就被坑过N次了 墙裂推荐在线 LaTeX 公式编辑器 min_25筛 用来求积性函数的前缀和 ∑i=1nF(i)\sum_{i=1}^{n}F(i)∑i=1n​F(i)，并且F(i)满足一些性质： ①当i为质数时F(ic)要存在一个较快的计算方法 ②当i...
min_25筛
myjs999的博客
12-24 480
看一遍绝对懂的博客。在我学完这个之后。
洲阁筛法学习小计
samjia2000的博客
01-14 3481
张俊的课件如是说：　　 一个对单个n有效的方法： 　　 　　f[n]表示n以内的素数个数 　　 　　c[i]表示第i个素数 　　 　　g[n][m]表示n以内不被c[1..m]整除的数的个数 　　　　　　(1) 　　　　　　f[n]=f[sqrt(n)]+g[n][f[sqrt(n)]]-1 　　　　　　//大于sqrt(n)的合数必然会被小于等于sqrt(n)的素数去掉，而1要
洲阁筛
xumingyang0的博客
01-10 691
来自
【spoj】DIVCNTK
weixin_30776863的博客
06-17 84
Portal -->Spoj DIVCNTK Solution 　　这题的话其实是。。洲阁筛模板题？差不多吧 　　题意就是给你一个函数\(S_k(x)\) \[ S_k(n)=\sum\limits_{i=1}^{n} \sigma_0(i^k) \] 　　其中\(\sigma_0(x)\)表示的是\(x\)的约数个数，现在已知\(n\)和\(k\)求\(S_k(n)\) mod \(2^{...
Min_25筛
Demon_Rieman的博客
07-10 891
前言： 因为我不会州阁筛，所以就只有学习一个他的一个相对简单的但是功能大致相同的筛法：min_25筛。 首先，其时间复杂度是n34lnnn34ln⁡n\frac{n^{\frac{3}{4}} }{\ln n}的。其功能是求出某个积性函数f(x)f(x)f(x)的前缀和， 并且需要满足如下几个条件： {1.f(p)是一个多项式(Let p be a prime)2.f(pc)容易被算出{1...
http://localhost/api/xsy-auth-security/loginUser/login 403 (Forbidden)
最新发布
07-28
- *2* [c/c++使用libcurl库做http客户端及封装(HTTP_GET和HTTP_POST)](https://blog.csdn.net/xsy29000/article/details/103181267)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值