ZOJ 4029 Now Loading!!!

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=4029

这道题比赛的时候本想着能够跳过取整操作来快速得到答案

但事实上是跳不过的。。。

先看分母，因为是log操作，所以其实分母上的答案最大只有30+左右

那么就可以排序二分来处理分母了

至于怎么操作 可以自己想想。挺简单的

#include<bits/stdc++.h>
using namespace std;
long long n,m;
long long a[555555];
struct node {
	long long data;
	long long id;
}b[555555];
long long qz[555555];
long long c[555555];
long long mod=1e9;
bool cp(node x,node y)
{
	return x.data<y.data;
}
int main(){
	int T;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%lld%lld",&n,&m);
		for(int i=0;i<n;i++)
		scanf("%lld",&a[i]);
		for(int i=0;i<m;i++)
		scanf("%lld",&b[i].data),b[i].id=i+1;
		sort(b,b+m,cp);
		qz[m-1]=b[m-1].id;
		for(int i=m-2;i>=0;i--)
		qz[i]=qz[i+1]+b[i].id;
		for(int i=0;i<m;i++)
		c[i]=b[i].data;
		long long ans=0;
		for(int i=0;i<n;i++)
		{
			long long zs=-1;
			for(long long j=1;j<=30;j++)
			{
				double mb=pow(a[i],1.0/(j*1.0));
				long long mb1=mb;
				if(mb1!=mb)mb1++;
				long long pos=lower_bound(c,c+m,mb1)-c;
				//cout<<pos<<" "<<mb<<endl;
				if(pos==zs&&zs==0)break;
				if(pos>=m)continue;
				if(zs==-1)
				{
				     zs=pos;
					 long long sl=qz[pos];
					 long long an=a[i]/j;
					 ans=(ans+sl*an%mod)%mod;	
				}
				else{
					long long sl=qz[pos]-qz[zs];
					 long long an=a[i]/j;
					 ans=(ans+sl*an%mod)%mod;
					 zs=pos;
				} 
			}
		}
		cout<<ans%mod<<endl;
	}
	return 0;
}