BZOJ 2749: [HAOI2012]外星人

最新推荐文章于 2018-09-25 16:50:57 发布

cgh_Andy

最新推荐文章于 2018-09-25 16:50:57 发布

阅读量429

点赞数 1

分类专栏： bzoj 数论

本文链接：https://blog.csdn.net/cgh_Andy/article/details/70878148

版权

bzoj 同时被 2 个专栏收录

151 篇文章 0 订阅

订阅专栏

数论

5 篇文章 0 订阅

订阅专栏

跪烂啊。。神题
直接贴网上dalao的吧。。
这里写图片描述

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=100000;
inline int read(){
    char ch=getchar(); int x=0,f=1;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}
    return x*f;
}
int p[N/10],pl,f[N+2];
bool v[N+2];
void pre(int n){
    f[1]=1;
    for(int i=2;i<=n;++i){
        if(!v[i]) p[++pl]=i,f[i]=f[i-1];
        for(int j=1;j<=pl && p[j]*i<=n;++j){
            int k=p[j]*i; v[k]=1;
            f[k]=f[p[j]]+f[i];
            if(i%p[j]==0)break;
        }
    }
}
int main()
{
    pre(N); int T=read();
    while(T--){
        int n=read(); LL ans=1;
        for(int i=1;i<=n;++i){
            LL p=read(),q=read();
            if(p==2)ans--;
            ans+=f[p]*q;
        }
        printf("%lld\n",ans);
    }
    return 0;
}