[BZOJ3513][MUTC2013]idiots(FFT)

题目:

我是超链接

题解:

一个合格的三角形——–两小边之和>第三边
总情况数我们知道 C3n C n 3 ,这里要用容斥,总的-不合法的
我们设a[i]表示i长度的木棍有多少,两个木棍和为i的方案数为F[i]
那么 F[i]=a(ij)a(j) F [ i ] = ∑ a ( i − j ) ∗ a ( j ) ,对于这个首先当i为偶数的时候,可能选择了两个一样的木棍,要先减去选择一样的数量;然后要/2,因为ab选了ba也选了。
然后设g[i]表示长度大于等于i的木棍数量,不合法的数量就是 g[i]F[i] ∑ g [ i ] ∗ F [ i ]

为什么我们一定要用容斥呢?g[i]表示长度小于i的木棍数量不行吗?其实是不行的,因为长度小于i的时候,就有可能选上和这两个木棍一样的木棍,但是长度大于等于的话就不可能和这两个木棍一样了

一定不要忘记卷积的数组要开两倍以上。。还有a数组r数组要全面清零

代码:

#include <cmath> 
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
using namespace std;
const int N=300005;
const double pi=acos(-1.0);
struct complex
{
    double x,y;
    complex(double X=0,double Y=0){x=X;y=Y;}
}a[N];
int n,r[N],g[N];LL F[N];
complex operator +(complex a,complex b){return complex(a.x+b.x,a.y+b.y);}
complex operator -(complex a,complex b){return complex(a.x-b.x,a.y-b.y);}
complex operator *(complex a,complex b){return complex(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
void FFT(complex *a,int id)
{
    for (int i=0;i<n;i++)
      if (i<r[i]) swap(a[i],a[r[i]]);
    for (int k=1;k<n;k<<=1)
    {
        complex wn=complex(cos(pi/k),id*sin(pi/k));
        for (int i=0;i<n;i+=(k<<1))
        {
            complex w=complex(1,0);
            for (int j=0;j<k;j++,w=w*wn)
            {
                complex x=a[i+j],y=w*a[i+j+k];
                a[i+j]=x+y; a[i+j+k]=x-y;
            }
        }
    }
}
int main()
{
    int T;scanf("%d",&T);
    while (T--)
    {
        memset(g,0,sizeof(g));memset(r,0,sizeof(r));
        for (int i=0;i<N;i++) a[i]=complex(0,0);
        int fn,maxx=0,c;scanf("%d",&fn);
        for (int i=1;i<=fn;i++) scanf("%d",&c),g[c]++,a[c].x=g[c],maxx=max(maxx,c);
        double fz,fm; fm=(double)fn*(fn-1)*(double)(fn-2)/6.0; fz=fm; 
        fn=maxx*2;int L=0;
        for (n=1;n<=fn;n<<=1) L++;
        for (int i=0;i<n;i++) r[i]=(r[i>>1]>>1)|((i&1)<<L-1);
        FFT(a,1); 
        for (int i=0;i<=n;i++) a[i]=a[i]*a[i];
        FFT(a,-1);
        for (int i=0;i<=n;i++) 
        {
            F[i]=(LL)(a[i].x/n+0.5); 
            if (i%2==0) F[i]-=g[i>>1];
            F[i]>>=1; 
        }
        for (int i=maxx;i>=1;i--) g[i]+=g[i+1];//>=i的个数 
        for (int i=1;i<=maxx;i++) fz-=F[i]*g[i];
        printf("%.7lf\n",fz/fm);
    }
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值