ICPCアジア上海2019 Cトリプル(FFT +暴力)のための予備コンテスト
ポータル:https://nanti.jisuanke.com/t/41400
質問の意味:
あなたよう、どのように多くのトリプル(i、j、k)は、あなたを依頼する3つの配列、B、C、
\ [\ {配列} {開始左\ {1}を| A_ {I} {-B_ J} \権| \当量のC_ {K}、\テキスト{と}} \\ {\左| B_ {J} -C_ {K} \権| \当量のA_ {I}、\テキスト{と}} \ \ {\左| A_ {I } -C_ {K}右\ | \当量のB_ {J}} \端{アレイ} \]
ソリューション:
簡素化による不平等の上に、我々は得ることができます
私たちは、トリプルの数を要求\(a_iを、B_j、\ C_kを ) 三角形を形成することができます
したがって本発明の組成物は、三角形HDU4609に類似している(https://www.cnblogs.com/buerdepepeqi/p/11236100.html)
しかし、違いは、我々は3つの配列から選択する必要があるということです
だからここ複製に関連する問題の選択に、我々は、デエンファシスを考えます
+ bは、+でスティックの数の長さbの畳み込みを行うために再度それを取ると仮定されます
我々は、cがC + bの長さであり、三角形の最も長い辺、次いで、B、Cは、木材の三三角片を形成することができないが、数よりも大きいと仮定する
我々は、(A + B)は、この長さは、その後、フォームが現在の番組長構成することができない三角形の数(B)のC *の長さよりも大きい数を列挙する
したがって、これに応じて、我々はそれが三回の畳み込みを成し遂げることができ、cは最長時間はbが最も長い辺である場合には、最大長と法律の数の縁、そしてすべてのトリプルがあり、側面に列挙されていますマイナス違法トライアドグループの数は、トリプルの合法的な数であります
小規模な暴力はここでちょっとしたトリックだのですか??小規模な暴力場合は、Tません。
部屋の大虐殺について(ノスタルジア
のは、理由を分析してみましょうか?
我々は、一定のxとFFT複合体クラスと仮定します
T 100の場合には、我々は単に三回FFTを実行し、その三回正FFT、三回IDFTであり、それは一定の6
複雑さは最後です
\ [T * 2 ^ {ログ
\ X \\ = 100 * 2 ^ {18} * 18 * 6 = 2831155200 * 6 *ログ(2N)(2N)}] ただし、暴力の複雑さをn2
80小さなデータセット(1000)のセット20件の大規模なデータがあると仮定
複雑さは、その後で
\ [20 * 18 * 2 ^ {18} * 1000~1000 * 6 + 80 = 646231040 \]
したがって、大規模な暴力小さなFFTの複雑性が比較的良好でドロップ
複雑性分析のために(感謝NEの首長
コード:
/**
* ┏┓ ┏┓
* ┏┛┗━━━━━━━┛┗━━━┓
* ┃ ┃
* ┃ ━ ┃
* ┃ > < ┃
* ┃ ┃
* ┃... ⌒ ... ┃
* ┃ ┃
* ┗━┓ ┏━┛
* ┃ ┃ Code is far away from bug with the animal protecting
* ┃ ┃ 神兽保佑,代码无bug
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┗━━━┓
* ┃ ┣┓
* ┃ ┏┛
* ┗┓┓┏━┳┓┏┛
* ┃┫┫ ┃┫┫
* ┗┻┛ ┗┻┛
*/
// warm heart, wagging tail,and a smile just for you!
//
// _ooOoo_
// o8888888o
// 88" . "88
// (| -_- |)
// O\ = /O
// ____/`---'\____
// .' \| |// `.
// / \||| : |||// \
// / _||||| -:- |||||- \
// | | \ - /// | |
// | \_| ''\---/'' | |
// \ .-\__ `-` ___/-. /
// ___`. .' /--.--\ `. . __
// ."" '< `.___\_<|>_/___.' >'"".
// | | : `- \`.;`\ _ /`;.`/ - ` : | |
// \ \ `-. \_ __\ /__ _/ .-` / /
// ======`-.____`-.___\_____/___.-`____.-'======
// `=---='
// ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
// 佛祖保佑 永无BUG
#include <set>
#include <map>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <string>
#include <vector>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
typedef pair<int, int> pii;
typedef unsigned long long uLL;
#define ls rt<<1
#define rs rt<<1|1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define bug printf("*********\n")
#define FIN freopen("input.txt","r",stdin);
#define FON freopen("output.txt","w+",stdout);
#define IO ios::sync_with_stdio(false),cin.tie(0)
#define debug1(x) cout<<"["<<#x<<" "<<(x)<<"]\n"
#define debug2(x,y) cout<<"["<<#x<<" "<<(x)<<" "<<#y<<" "<<(y)<<"]\n"
#define debug3(x,y,z) cout<<"["<<#x<<" "<<(x)<<" "<<#y<<" "<<(y)<<" "<<#z<<" "<<z<<"]\n"
const int maxn = 5e5 + 5;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const double Pi = acos(-1);
LL gcd(LL a, LL b) {
return b ? gcd(b, a % b) : a;
}
LL lcm(LL a, LL b) {
return a / gcd(a, b) * b;
}
double dpow(double a, LL b) {
double ans = 1.0;
while(b) {
if(b % 2)ans = ans * a;
a = a * a;
b /= 2;
} return ans;
}
LL quick_pow(LL x, LL y) {
LL ans = 1;
while(y) {
if(y & 1) {
ans = ans * x % mod;
} x = x * x % mod;
y >>= 1;
} return ans;
}
LL res[maxn << 2], len;
struct Complex {
double r, i;
Complex(double r = 0, double i = 0) : r(r), i(i) {};
Complex operator+(const Complex &rhs) {
return Complex(r + rhs.r, i + rhs.i);
}
Complex operator-(const Complex &rhs) {
return Complex(r - rhs.r, i - rhs.i);
}
Complex operator*(const Complex &rhs) {
return Complex(r * rhs.r - i * rhs.i, i * rhs.r + r * rhs.i);
}
} va[maxn << 2], vb[maxn << 2];
void rader(Complex F[], int len) { //len = 2^M,reverse F[i] with F[j] j为i二进制反转
int j = len >> 1;
for (int i = 1; i < len - 1; ++i) {
if (i < j) swap(F[i], F[j]); // reverse
int k = len >> 1;
while (j >= k) {
j -= k;
k >>= 1;
}
if (j < k) j += k;
}
}
void FFT(Complex F[], const int &len, const int &t) {
rader(F, len);
for (int h = 2; h <= len; h <<= 1) {
Complex wn(cos(-t * 2 * Pi / h), sin(-t * 2 * Pi / h));
for (int j = 0; j < len; j += h) {
Complex E(1, 0); //旋转因子
for (int k = j; k < j + h / 2; ++k) {
Complex u = F[k];
Complex v = E * F[k + h / 2];
F[k] = u + v;
F[k + h / 2] = u - v;
E = E * wn;
}
}
}
if (t == -1) //IDFT
for (int i = 0; i < len; ++i)
F[i].r /= len;
}
void Conv(Complex a[], Complex b[], const int &len) { //求卷积
FFT(a, len, 1);
FFT(b, len, 1);
for (int i = 0; i < len; ++i) a[i] = a[i] * b[i];
FFT(a, len, -1);
}
void work() {
Conv(va, vb, len);
for (int i = 0; i < len; ++i)res[i] = va[i].r + 0.5;
}
int a[maxn], b[maxn], c[maxn];
int numa[maxn], numb[maxn], numc[maxn];
int suma[maxn], sumb[maxn], sumc[maxn];//长度为i的数量
LL resab[maxn], resbc[maxn], resac[maxn];//长度为i的 a,b的个数
int main() {
#ifndef ONLINE_JUDGE
FIN
#endif
int T;
int cas = 1;
scanf("%d", &T);
while(T--) {
int n;
scanf("%d", &n);
int maxlena = 0, maxlenb = 0, maxlenc = 0;
for(int i = 1; i <= n; i++) {
scanf("%d", &a[i]);
numa[a[i]]++;
maxlena = max(maxlena, a[i]);
}
for(int i = 1; i <= n; i++) {
scanf("%d", &b[i]);
numb[b[i]]++;
maxlenb = max(maxlenb, b[i]);
}
for(int i = 1; i <= n; i++) {
scanf("%d", &c[i]);
numc[c[i]]++;
maxlenc = max(maxlenc, c[i]);
}
LL ans = 1ll * n * n * n;
if(n <= 1000) {
for(int i = 1; i <= maxlena; i++) {
suma[i] = suma[i - 1] + numa[i];
}
for(int i = 1; i <= maxlenb; i++) {
sumb[i] = sumb[i - 1] + numb[i];
}
for(int i = 1; i <= maxlenc; i++) {
sumc[i] = sumc[i - 1] + numc[i];
}
for(int i = 1; i <= n; i++) {
for(int j = 1; j <= n; j++) {
resab[a[i] + b[j]]++;
resbc[b[i] + c[j]]++;
resac[a[i] + c[j]]++;
}
}
int mx = max(maxlena, max(maxlenb, maxlenc));
for(int i = 1; i <= mx; i++) {
if(i <= maxlenc) ans -= resab[i] * (sumc[maxlenc] - sumc[i]);
if(i <= maxlenb) ans -= resac[i] * (sumb[maxlenb] - sumb[i]);
if(i <= maxlena) ans -= resbc[i] * (suma[maxlena] - suma[i]);
}
memset(suma, 0, sizeof(int) * (maxlena + 2));
memset(sumb, 0, sizeof(int) * (maxlenb + 2));
memset(sumc, 0, sizeof(int) * (maxlenc + 2));
memset(numa, 0, sizeof(int) * (maxlena + 2));
memset(numb, 0, sizeof(int) * (maxlenb + 2));
memset(numc, 0, sizeof(int) * (maxlenc + 2));
memset(resac, 0, sizeof(LL) * (mx * 2 + 2));
memset(resab, 0, sizeof(LL) * (mx * 2 + 2));
memset(resbc, 0, sizeof(LL) * (mx * 2 + 2));
} else {
maxlena++, maxlenb++, maxlenc++;
len = 1;
int mxab = max(maxlena, maxlenb);
while(len < 2 * mxab) len <<= 1;
for(int i = 0; i < len; i++) {
if (i < mxab) {
va[i] = Complex(numa[i], 0);
vb[i] = Complex(numb[i], 0);
} else
va[i] = vb[i] = Complex(0, 0);
}
work();
for (int i = 0; i < len; i++) resab[i] = res[i];
len = 1;
int mxac = max(maxlena, maxlenc);
while(len < 2 * mxac) len <<= 1;
for(int i = 0; i < len; i++) {
if (i < mxac) {
va[i] = Complex(numa[i], 0);
vb[i] = Complex(numc[i], 0);
} else
va[i] = vb[i] = Complex(0, 0);
}
work();
for (int i = 0; i < len; i++) resac[i] = res[i];
len = 1;
int mxbc = max(maxlenb, maxlenc);
while(len < 2 * mxbc) len <<= 1;
for(int i = 0; i < len; i++) {
if (i < mxbc) {
va[i] = Complex(numb[i], 0);
vb[i] = Complex(numc[i], 0);
} else
va[i] = vb[i] = Complex(0, 0);
}
work();
for (int i = 0; i < len; i++) resbc[i] = res[i];
for(int i = 1; i <= maxlena; i++) {
suma[i] = suma[i - 1] + numa[i];
}
for(int i = 1; i <= maxlenb; i++) {
sumb[i] = sumb[i - 1] + numb[i];
}
for(int i = 1; i <= maxlenc; i++) {
sumc[i] = sumc[i - 1] + numc[i];
}
for (int i = 1; i <= 2 * mxab; ++i) {
if (i > maxlenc) break;
ans -= resab[i] * (sumc[maxlenc] - sumc[i]);
}
for (int i = 1; i <= 2 * mxbc; ++i) {
if (i > maxlena) break;
ans -= resbc[i] * (suma[maxlena] - suma[i]);
}
for (int i = 1; i <= 2 * mxac; ++i) {
if (i > maxlenb) break;
ans -= resac[i] * (sumb[maxlenb] - sumb[i]);
}
memset(suma, 0, sizeof(int) * (maxlena + 2));
memset(sumb, 0, sizeof(int) * (maxlenb + 2));
memset(sumc, 0, sizeof(int) * (maxlenc + 2));
memset(numa, 0, sizeof(int) * (maxlena + 2));
memset(numb, 0, sizeof(int) * (maxlenb + 2));
memset(numc, 0, sizeof(int) * (maxlenc + 2));
memset(resac, 0, sizeof(LL) * (mxac * 2 + 2));
memset(resab, 0, sizeof(LL) * (mxab * 2 + 2));
memset(resbc, 0, sizeof(LL) * (mxbc * 2 + 2));
}
printf("Case #%d: %lld\n", cas++, ans);
}
return 0;
}