高仿圣罗兰包包价格
蓝气球机械表怎么样
高仿黑水鬼机械表怎么样
高仿帕玛强尼价格
高仿普拉达包包多少钱
皇家橡树机械表
江诗丹顿手表怎么样
高仿七个星期五机械表怎么样
高仿百年灵机械表价格
宝格丽机械表怎么样
普拉达女包
三叶草运动鞋价格
高仿新百伦篮球鞋
高仿劳力士机械表怎么样
真力时价格
高仿ASICS运动鞋
高仿宝玑机械表
高仿宝珀手表价格
宝格丽机械表多少钱
高仿宇舶男表价格
高仿七个星期五表价格
Dior包怎么样
宝珀表
陀飞轮男表
高仿浪琴机械表多少钱
高仿手表价格
loewe包价格
帕玛强尼男表价格
冠军鞋
高仿劳力士手表怎么样
高仿彪马运动鞋
高仿理查德米勒男士手表多少钱
高仿安德玛鞋价格
百年灵机械表怎么样
劳力士多少钱
积家多少钱
江诗丹顿怎么样
浪琴女表价格
高仿匡威篮球鞋怎么样
高仿泰格豪雅怎么样
高仿圣罗兰怎么样
理查德米勒男表价格
高仿欧米茄手表怎么样
高仿劳力士表
高仿卡地亚机械表怎么样
高仿江诗丹顿男表价格
七个星期五男士手表
高仿包
宝格丽机械表价格
高仿路易威登多少钱
欧米茄女表
圣罗兰
高仿男士手表价格
格拉苏蒂机械表怎么样
范思哲价格
高仿亚瑟士篮球鞋多少钱
高仿沛纳海表多少钱
高仿百年灵表怎么样
高仿江诗丹顿手表价格
理查德米勒表怎么样
芬迪包包价格
百达翡丽女士手表
高仿卡西欧男士手表
高仿亚瑟士
高仿阿迪运动鞋多少钱
浪琴男士手表怎么样
美度多少钱
高仿古奇女包多少钱
罗意威男包多少钱
高仿宝格丽女士手表多少钱
高仿欧米茄男表价格
高仿女包
高仿黑水鬼价格
高仿宝格丽怎么样
高仿格拉苏蒂手表怎么样
高仿理查德米勒手表怎么样
高仿浪琴
新百伦鞋怎么样
高仿瑞士
瑞士名表多少钱
高仿Armani包包多少钱
宝格丽女表多少钱
高仿欧米茄男表
万国机械表怎么样
爱马仕女包
高仿华伦天奴包包价格
高仿七个星期五表多少钱
高仿泰格豪雅男士手表怎么样
帕玛强尼机械表多少钱
真力时手表价格
高仿百达翡丽价格
高仿七个星期五男表
高仿卡西欧男士手表价格
宝格丽女表
绿水鬼怎么样
高仿七个星期五机械表价格
爱彼表
高仿蓝气球手表
真力时男士手表
高仿宇舶男士手表怎么样
宇舶手表价格
高仿爱彼手表价格
帕玛强尼机械表
高仿百达翡丽表
高仿黑水鬼机械表
阿迪篮球鞋怎么样
帕玛强尼表价格
高仿天梭男表价格
Dior包包怎么样
高仿VANS篮球鞋
高仿黑水鬼表价格
欧米茄表
香奈儿女士手表多少钱
卡西欧怎么样
江诗丹顿男士手表多少钱
高仿百年灵手表多少钱
nike怎么样
帕玛强尼机械表怎么样
高仿卡西欧男士手表怎么样
宇舶表怎么样
高仿圣罗兰包包怎么样
宇舶表多少钱
蓝气球男表多少钱
高仿芝柏男士手表
高仿卡地亚表多少钱
耐克多少钱
高仿香奈儿多少钱
高仿卡地亚机械表
高仿hermes女包多少钱
高仿七个星期五机械表
高仿帝舵价格
高仿沛纳海机械表价格
ugg
高仿圣罗兰包包怎么样
百达翡丽男士手表多少钱
高仿名表多少钱
高仿手表
新百伦篮球鞋价格
江诗丹顿男士手表怎么样
高仿爱马仕包包怎么样
高仿宝玑男士手表多少钱
高仿宝格丽
伯爵女表多少钱
高仿手表
高仿VANS多少钱
Dior包包价格
高仿爱彼男士手表价格
宝玑多少钱
高仿蓝气球表价格
高仿万国机械表价格
高仿古奇包包
沛纳海男士手表价格
高仿名表
高仿安德玛鞋
路易威登怎么样
彪马运动鞋价格
万宝龙男士手表价格
高仿Dior包包怎么样
伯爵女表
高仿耐克价格
高仿天梭表
高仿伯爵男表
高仿美度手表
高仿宝格丽机械表多少钱
高仿阿迪多少钱
高仿黑水鬼男表怎么样
沛纳海男士手表
prada男包价格
高仿prada包包怎么样
高仿七个星期五手表怎么样
高仿浪琴女士手表价格
高仿百达翡丽机械表
宝格丽包怎么样
COACH怎么样
高仿宝格丽男表价格
浪琴手表多少钱
高仿浪琴男士手表怎么样
PUMA运动鞋多少钱
高仿罗意威女包多少钱
伯爵表价格
斐乐鞋多少钱
高仿宝珀男表怎么样
爱彼男士手表
高仿爱马仕男包价格
芝柏手表怎么样
高仿陀飞轮机械表价格
高仿万国机械表怎么样
高仿乔丹运动鞋价格
绿水鬼手表价格
万国机械表价格
ugg女包价格
y3篮球鞋
高仿爱彼机械表
高仿劳力士表价格
七个星期五
浪琴男表价格
高仿宝格丽男表
沛纳海价格
高仿gucci包价格
七个星期五多少钱
百达翡丽女表怎么样
高仿罗意威男包价格
高仿格拉苏蒂男士手表怎么样
gucci包
高仿斐乐多少钱
欧米茄男士手表怎么样
高仿VANS运动鞋
阿玛尼女包怎么样
高仿泰格豪雅价格
高仿卡西欧男表多少钱
蓝气球女士手表怎么样
loewe女包
女包
路易威登男包
宇舶手表
百年灵男表
高仿天梭表怎么样
香奈儿女表多少钱
真力时男士手表价格
亚瑟士篮球鞋怎么样
高仿帕玛强尼男表
高仿江诗丹顿男表
爱彼机械表多少钱
百年灵手表价格
宝玑表价格
沛纳海表多少钱
高仿美度手表多少钱
宝格丽男士手表价格
黑水鬼
高仿耐克篮球鞋怎么样
高仿鞋
耐克运动鞋
高仿耐克篮球鞋
Fendi怎么样
格拉苏蒂机械表多少钱
高仿普拉达包价格
积家男士手表多少钱
mcm多少钱
罗意威包包怎么样
ugg包
高仿宝格丽机械表多少钱
COACH女包多少钱
Dior包包多少钱
高仿美度机械表
高仿积家女表多少钱
帕玛强尼
蓝气球男士手表
瑞士手表价格
黑水鬼多少钱
七个星期五男表怎么样
aj鞋
高仿万宝龙男士手表
迪奥女包多少钱
帝舵手表价格
蓝气球多少钱
高仿天梭机械表价格
欧米茄男士手表多少钱
高仿巴宝莉包包多少钱
高仿泰格豪雅多少钱
高仿ysl包包
沛纳海男士手表
高仿包
高仿loewe包
高仿普拉达男包价格
高仿伯爵男表
高仿皇家橡树机械表价格
泰格豪雅手表怎么样
七个星期五手表价格
古驰 朗格
英特尔与腾讯宣布战略合作,推动智慧教育行业发展,
齐鲁热线   2019-12-06 21:10   
luxurywatches:
VANS
河南通报21批次不合格食品 涉及微生物污染、食品添加剂、生物毒素污染等问题,  云南两期抽检通报8批次不合格食品,涉及微生物污染、兽药残留问题,  

  比照效果明隐的欧洲劣同青训体系,zhongguo足球的青训体系借正在初初阶段,但出有管若何的初初阶段,亦出有会对中小型职业俱乐部正在青训体系圆里收力组成阻碍——青训的重要性大家皆知,特别U15国少、U17国青、U19国青曾经沦降到亚洲三流水平,国奥队一个月后U23亚洲杯征程也是凶多凶少,zhongguo足球有足够的来由实正低头耕作青训,现正在问题的闭键是青训需供正确的理念战时间,即使足球文明基果强如本菲卡,也是正在苦守10年沉训以后才结出硕果。

  

  真力时“做为教死,我们秉持了志愿者‘奉献、友爱、相助、进步’的细神。看到特殊孩子们能取普通孩子一起踢球,那是值得我们欣喜的。”吴洪杰也讲出有上本果,便是觉得那群孩子的笑很有感染力,“当看到孩子们正在足球场上悲笑我们也会出有由自登时跟着笑起去”。

  而正在更多的时候,我们正在一些做品中看到的,恰恰是由于短少具有故事代价的场景,让那些故事项得空洞战空洞。好比做品家丁公的一个重要决意战让人正在乎的转变,出有展垫到位,实在出有实正在可托,当然也流于套路,出有挨动读者的力量。

本文编辑:北京热线010

雅克德罗
  Montblanc MCHR又双叒获奖了!CEO李超先生荣获2019中国人力资源科技影响力TOP人物大奖 ,
  呼ね    
 
 
 
  理查( )BellRoss luxurywatches  miumiu    
 
BVLGARI
   
泰格豪雅这8省注意!2020年二建不限专业就能报考!,
阿玛尼亨达外汇:观望非农就业数据美元短线偏软,
帕玛强尼2020年广东省选调生成绩单在哪看:http://ggfw.gdhrss.gov.cn/yxks,
ysl这8省注意!2020年二建不限专业就能报考!,
古奇全力以赴抓好重点领域安全生产工作 ,
CHANEL郭姓股民向东旭光电发起索赔 张志旺律师接单,
理查长江三峡通航管理局物业服务项目开标时间延期通知其它,
施华洛九都镇开展《宪法》知识宣传活动,
普拉达半开放式设计+主动降噪:华为FreeBuds 3创造TWS耳机新体验,
blancpain高通骁龙765/骁龙765G详解:首次集成5G、最高速3.7Gbps,
雅典湖南“一带一路”绿色博览会暨绿色产业发展论坛开幕,
帝舵长江三峡通航管理局物业服务项目开标时间延期通知其它,
万国两企业红线重叠致项目无法推进,
JaquetDroz伊顿EGE3040FFG新优惠价格 价格咨询请联系:毛女士13671505969微信同号,
 
ASICS
   
卡西欧共享自习室行业发展亟须规范,
Glashutte科技早报 小米电视5Pro于12日开售;华米X概念手表明年量产,
rolex致个人所得税扣缴单位的一封信,
CHANEL原油交易提醒:OPEC+額外減產50萬桶!100%履行減產?壓力山大,美油徘徊兩個月高點,
雅克德罗华为nova6发布会直播地址 5g新品手机发布会观看入口,
古驰河北唐山发生4.5级地震 北京天津等地有震感,
法穆兰天津滨海新区建设三级志愿服务管理体系,
阿迪达斯吉喆在北京效力时的众多队友发文缅怀,
彪马2020款丰田酷路泽4000丐版最低配最低价59,
卡地亚工作计划书格式书写,
亚瑟士用 python 分析基金!让赚钱赢在起跑线!,
万国【工程进度】邵阳这个公园即将开园!美景提前曝光!,
iwceclipse 克隆 https 地址的 Git 仓库报错:cannot open git-upload-pack,
mcm8岁女童吃“伟哥”续命背后:数百万“蓝嘴唇”病人艰难求生,
阿迪达斯
  
阿特金森:普林斯很优秀 恩瓦巴会得到更多的机会,
2019甘肃省生态环境科学设计研究院编外招聘16人公告(第二期),
宝珀
 
Breguet2020泰安事业单位面试热点:为“网生代”寻找新思路——未成年人保护,
dwMCHR又双叒获奖了!CEO李超先生荣获2019中国人力资源科技影响力TOP人物大奖 ,
BVLGARI湛江等单位消防车项目公开招标公告(GPCGD19B800HG057F),
loewe美国海军珍珠港造船厂发生枪击事件致3人受伤 枪手疑自杀,
百达翡丽黄河颂:百折不挠,奔涌向前12邱承彬诗邱木兮诗邱颂,
Breguet苹果开始购买无碳铝 计划在部分产品中使用,
爱马仕好消息,“二师兄”又降价了!7000亿猪肉股却慌了,
圣罗兰
  
耐克北京奥林匹克公园2020年春节部分道路景观亮化项目更正公告,
A.Lange&S?hne体坛:除卡拉斯科外,奥斯卡&伊哈洛也想重返欧洲,
阿玛尼墨菲:阿里展现出了最好的一面,他对曼联的进球真是太精彩了,
沛纳海《实况足球2020》更新包3.0已推出 含球场/球员面貌更新,
mk苹果开始购买无碳铝 计划在部分产品中使用,
JaquetDroz永州:全面推进中心城区周边乡镇和城乡结合部环境卫生整治攻坚工作,
古琦莫丞雨:12.5黄金晚间初请会继续破新高吗?明日非农如何操作?,
法穆兰
 
爱马仕[影像标杆 索尼A7R3浙江售价仅14200元,
blancpain
  
连云港客厅风景自粘壁画(尺寸定制)(销售公司+欢迎您 ,
莫丞雨:12.5黄金晚间初请会继续破新高吗?明日非农如何操作?,
中卫沥青麻筋)厂家批发/有限公司,
莫丞雨:12.5黄金晚间初请会继续破新高吗?明日非农如何操作?,
法穆兰
   
亨达外汇:观望非农就业数据美元短线偏软,
三亚市人力资源和社会保障局-三亚市人力资源和社会保障局服务大厅和办公楼室内装修改造修缮工程-竞争性磋商公告,
NBA新秀榜:帕斯卡尔升至榜首 八村塁第4,
湛江等单位消防车项目公开招标公告(GPCGD19B800HG057F),
  bally   
  dw   
Glashutte
美国最高金融小组希望监控数字资产,
外贸沟通谈判技巧,
---labbasenet.cn---
首页 > 其他好文 > 详细

EOJ Monthly 2019.11 E. 数学题(莫比乌斯反演+杜教筛+拉格朗日插值)

时间:2019-12-06 23:52:43      阅读:15      评论:0      收藏:0      [点我收藏+]
百度新闻是包含海量高仿手表资讯的网上购物新闻服务平台,真实反映每时每刻的>高仿包包新闻热点。您可以搜索新闻事件、>高仿手表热点话题、电子商务、产品资讯等,快速了解它们的最新进展byalex。

标签:inline   void   ons   const   复杂度   time   second   iostream   return   

传送门

题意:
统计\(k\)元组个数\((a_1,a_2,\cdots,a_n),1\leq a_i\leq n\)使得\(gcd(a_1,a_2,\cdots,a_k,n)=1\)
定义\(f(n,k)\)为满足要求的\(k\)元组个数,现在要求出\(\sum_{i=1}^n f(i,k),1\leq n\leq 10^9,1\leq k\leq 1000\)

思路:
首先来化简一下式子,题目要求的就是:
\[ \begin{aligned} &\sum_{i=1}^n\sum_{j=1}^n\cdots \sum_{k=1}^n gcd(i,j,\cdots, k,n)=1\=&\sum_{i=1}^n\sum_{j=1}^n\cdots \sum_{k=1}^n\sum_{d|i,j,\cdots,k,n}\mu(d)\=&\sum_{d|n}\mu(d)\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{n}{d}}\cdots \sum_{k=1}^\frac{n}{d}1\=&\sum_{d|n}\mu(d) (\frac{n}{d})^k \end{aligned} \]
套路到此结束~然后观察到这个式子其实是一个狄利克雷卷积的形式,\(f(i)=\mu(i),g(i)=i^k\),上式则为:\(f*g_{(n)}\)
那么题目要求的就是这个卷积的前缀和,注意两个积性函数的卷积也是积性函数,因为\(\mu*I=\varepsilon\),所以我们再构造一个积性函数\(h=I\),直接上杜教筛就行了。最后的式子是:
\[ h(1)\cdot S(n) = \sum_{i=1}^ng(i)-\sum_{d=2}^n h(d)S(\lfloor\frac{n}{d}\rfloor) \]
后半部分直接整除分块,至于\(\sum_{i=1}^ng(i)\),拉格朗日插值能在\(O(k)\)的时间复杂度解决。
代码如下(比赛的时候写得稍微有点乱):

/*
 * Author:  heyuhhh
 * Created Time:  2019/11/29 21:03:32
 */
#include <iostream>
#include <algorithm>
#include <vector>
#include <cmath>
#include <set>
#include <map>
#include <iomanip>
#define MP make_pair
#define fi first
#define se second
#define sz(x) (int)(x).size()
#define all(x) (x).begin(), (x).end()
#define INF 0x3f3f3f3f
#define Local
#ifdef Local
  #define dbg(args...) do { cout << #args << " -> "; err(args); } while (0)
  void err() { std::cout << '\n'; }
  template<typename T, typename...Args>
  void err(T a, Args...args) { std::cout << a << ' '; err(args...); }
#else
  #define dbg(...)
#endif
void pt() {std::cout << '\n'; }
template<typename T, typename...Args>
void pt(T a, Args...args) {std::cout << a << ' '; pt(args...); }
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
//head
const int N = 1e4 + 5, MOD = 998244353;

int n, k;
ll qpow(ll a, ll b) {
    ll ans = 1;
    while(b) {
        if(b & 1) ans = ans * a % MOD;
        a = a * a % MOD;
        b >>= 1;   
    }
    return ans;   
}
struct Lagrange {
    static const int SIZE = 1005;
    ll f[SIZE], fac[SIZE], inv[SIZE], pre[SIZE], suf[SIZE];
    int n;
    inline void add(ll &x, int y) {
        x += y;
        if(x >= MOD) x -= MOD;
    }
    void init(int _n) {
        n = _n;
        fac[0] = 1;
        for (int i = 1; i < SIZE; ++i) fac[i] = fac[i - 1] * i % MOD;
        inv[SIZE - 1] = qpow(fac[SIZE - 1], MOD - 2);
        for (int i = SIZE - 1; i >= 1; --i) inv[i - 1] = inv[i] * i % MOD;
        //设置f初值,可以根据需要修改
        f[0] = 0;
        for (int i = 1; i <= n; ++i)
            f[i] = (f[i - 1] + qpow(i, k)) % MOD;
    }
    ll calc(ll x) {
        if (x <= n) return f[x];
        pre[0] = x % MOD;
        for (int i = 1; i <= n; ++i) pre[i] = pre[i - 1] * ((x - i) % MOD) % MOD;
        suf[n] = (x - n) % MOD;
        for (int i = n - 1; i >= 0; --i) suf[i] = suf[i + 1] * ((x - i) % MOD) % MOD;
        ll res = 0;
        for (int i = 0; i <= n; ++i) {
            ll tmp = f[i] * inv[n - i] % MOD * inv[i] % MOD;
            if (i) tmp = tmp * pre[i - 1] % MOD;
            if (i < n) tmp = tmp * suf[i + 1] % MOD;
            if ((n - i) & 1) tmp = MOD - tmp;
            add(res, tmp);
        }
        return res;
    }
}lagrange;
int mu[N], p[N];
bool chk[N];
int pre[N];
void init() {
    mu[1] = 1;
    int cnt = 0;
    for(int i = 2; i <= N - 1; i++) {
        if(!chk[i]) p[++cnt] = i, mu[i] = -1;
        for(int j = 1; j <= cnt && i * p[j] <= N - 1; j++) {
            chk[i * p[j]] = 1;
            if(i % p[j] == 0) {mu[i * p[j]] = 0; break;}
            mu[i * p[j]] = -mu[i]; 
        }
    }
    for(int i = 1; i <= N - 1; i++) {
        int res = 0;
        for(int j = 1; 1ll * j * j <= i; j++) {
            if(i % j == 0) {
                int d1 = j, d2 = i / j;
                res = (res + 1ll * mu[d1] * qpow(d2, k) % MOD) % MOD;
                if(d1 != d2) res = (res + 1ll * mu[d2] * qpow(d1, k) % MOD) % MOD;
                if(res < 0) res += MOD;
            }
        }   
        pre[i] = (pre[i - 1] + res) % MOD;
    }
}
map <int, ll> mp;
ll djs(int n) {
    if(n < N) return pre[n];
    if(mp.find(n) != mp.end()) return mp[n];
    ll ans = lagrange.calc(n);
    for(int i = 2, j; i <= n; i = j + 1) {
        j = n / (n / i);
        ans -= 1ll * (j - i + 1) * djs(n / i) % MOD;
        if(ans < 0) ans += MOD;   
    }
    return mp[n] = ans;
}

void run(){
    lagrange.init(k + 1);
    init();
    int ans = djs(n);
    cout << ans << '\n';
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0); cout.tie(0);
    cout << fixed << setprecision(20);
    while(cin >> n >> k) run();
    return 0;
}

EOJ Monthly 2019.11 E. 数学题(莫比乌斯反演+杜教筛+拉格朗日插值)

标签:inline   void   ons   const   复杂度   time   second   iostream   return   

原文地址:https://www.cnblogs.com/heyuhhh/p/11973681.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
? 2014 mamicode.com 版权所有 京ICP备13008772号-2
迷上了代码!
labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网
labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网
labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网
labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网 labbasenet.cn新闻网