2019牛客暑期多校训练营（第一场）B.Integration（留数定理）

2019牛客暑期多校训练营（第一场）B.Integration（留数定理）

题目大意

总所周知
$${\int }_0^{\infty }\frac{1}{1+x^2}dx=\frac{\pi }{2}$$
,现求
$$\frac{1}{\pi }{\int }_0^{\infty }\frac{1}{{\prod }_{i=1}^n(a_i^2+x^2)}dx$$

解题思路

对原式进行裂项，根据留数定理，第j项就是
$$\frac{1}{\pi {\prod }_{i=1}^{n\&\&i!=j}(a_i^2-a_j^2)}{\int }_0^{\infty }\frac{1}{a_j^2+x^2}dx$$
而${\int }_0^{\infty }\frac{1}{a_j^2+x^2}dx=\frac{\pi }{2a_j}$

则可以得出答案

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int mod =1e9+7;
const int size=1e3+5;
int quick_pow(int a,int b)
{
	int ans=1;
	while(b)
	{
		if(b&1) ans=1LL*ans*a%mod;
		a=1LL*a*a%mod;
		b>>=1;
	}
	return ans%mod;
}
int a[size];
int k[size];
int main()
{
 	int n;
	while(~scanf("%d",&n))
	{	
		for(int i=1;i<=n;i++)
		{
			scanf("%d",&a[i]);
		}
		for(int i=1;i<=n;i++)
		{
			k[i]=1;
			for(int j=1;j<=n;j++)
			{
				if(j==i) continue;
				k[i]=1LL*k[i]*((1LL*a[j]*a[j]%mod-1LL*a[i]*a[i]%mod+mod)%mod)%mod;
			}
			k[i]=quick_pow(k[i],mod-2);
		}
		int ans=0;
		for(int i=1;i<=n;i++)
		{
			ans=(ans+1LL*k[i]*500000004%mod*quick_pow(a[i],mod-2)%mod)%mod;
		}
		printf("%d\n",ans);
	}
}