小猪佩奇玩游戏，LGR-060-div1-1，简单的数论

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接： https://blog.csdn.net/Deep_Kevin/article/details/102586346

正题

      Portal

      考虑每一个数被自己删去的概率，因为如果被别人删去了，那么这个贡献就不用算了。

       那么就是对于每一个i，若它是x个数的幂，那么概率就是，答案就是

       这样就做完了，我们用枚举i来遍历他的次幂，然后更新答案，因为除了1次幂，这样的数并不多。

#include<bits/stdc++.h>
using namespace std;

int T,n;
map<int,int> mp;
map<int,int>::iterator it;
struct node{
	int n,id;
	double ans;
	bool operator<(const node q)const{
		return n<q.n;
	}
}s[110];
double ans[110];

int main(){
	scanf("%d",&T);
	for(int i=1;i<=T;i++) scanf("%d",&s[i].n),s[i].id=i;
	sort(s+1,s+1+T);
	for(int i=1;i<=T;i++) s[i].ans=s[i].n-s[i-1].n;
	for(int i=2;i*i<=s[T].n;i++){
		for(long long j=i*i;j<=s[T].n;j*=i){
			node*x=lower_bound(s+1,s+T+1,(node){j,0,0});
			if(!mp[j]) mp[j]=0;
			(*x).ans-=1.0/(mp[j]+1);
			mp[j]++;
			(*x).ans+=1.0/(mp[j]+1);
		}
	}
	for(int i=1;i<=T;i++) s[i].ans+=s[i-1].ans,ans[s[i].id]=s[i].ans;
	for(int i=1;i<=T;i++) printf("%.8lf\n",ans[i]);
}