ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

[SDOI2008]仪仗队

2019-05-30 14:03:21  阅读:217  来源: 互联网

标签:ch const 仪仗队 int long SDOI2008 include 4e4


传送门:

分析完题目,不就是要找y=kx有多少条0.0,那就是要找 (n-1)内gcd(i,j)==1的个数吗,最后答案+2就行了,直接走莫比乌斯反演。。。(题解区并没有人用)

最后特判n==1输出0就行了

#include <cstdio>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <iostream>
using namespace std;
#define ll long long
#define re register
const int N=4e4+10;
const int M=4e4;
inline void read(int &a)
{
    a=0;
    int d=1;
    char ch;
    while(ch=getchar(),ch>'9'||ch<'0')
        if(ch=='-')
            d=-1;
    a=ch^48;
    while(ch=getchar(),ch>='0'&&ch<='9')
        a=(a<<3)+(a<<1)+(ch^48);
    a*=d;
}
int mu[N],pri[N],cnt;
bool vis[N];
inline void init()
{
    mu[1]=1;
    for(re int i=2;i<=M;i++)
    {
        if(!vis[i])
            pri[++cnt]=i,mu[i]=-1;
        for(re int j=1;j<=cnt&&i*pri[j]<=M;j++)
        {
            vis[i*pri[j]]=1;
            if(i%pri[j]==0)
                break;
            mu[i*pri[j]]=-mu[i];
        }
    }
    for(re int i=1;i<=M;i++)
        mu[i]+=mu[i-1];
}
inline ll solve(int n)
{
    ll ans=0;
    for(int l=1,r;l<=n;l=r+1)
    {
        r=n/(n/l);
        ans+=(mu[r]-mu[l-1])*(n/l)*(n/l);
    }
    return ans;
}
int main()
{
    init();
    int n;
    read(n);
    if(n==1)
        return puts("0"),0;
    n--;
    printf("%lld",solve(n)+2);
    return 0;
}

 

标签:ch,const,仪仗队,int,long,SDOI2008,include,4e4
来源: https://www.cnblogs.com/acm1ruoji/p/10949182.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有