叉姐的魔法训练 （一） 初级魔法练习

https://vjudge.net/contest/32389 这是链接地址，顺着作者就能找到其他的训练了
POJ - 2443 Set Operation
这道题的题意是：现在有N个集合，每个集合中的数字可能有重复的，现在对于每个询问i，j,判断这两个数是否同时出现在一个集合中。

题意很简洁，但是数据范围就让我们无法暴力过题。我们可以使用Bitset来做这个事情。我们给每个数字都开一个bitset，其中如果这个数出现在第k个集合中，那么就这个数的bitset对应位置就是1.维护完之后，我们就可以用求交集的办法来接这个问题了

#include <cstdio>
#include <map>
#include <bitset>
#include <cstring>
#define fi first
#define se second
using namespace std;

int main()
{
    map<int,bitset<1001> > in;
    int n;
    while(scanf("%d",&n) != EOF)
    {
        in.clear();
        for(int i = 1; i<=n; i++)
        {
            int c;

            scanf("%d",&c);
            while(c--)
            {
                int x;
                scanf("%d",&x);
                in[x][i] = 1;
            }
        }
        int q;
        scanf("%d",&q);
        while(q--)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            //cout<<in[x]<<endl;
            bitset<1001> temp = in[x] & in[y];
            if(temp.any()) printf("Yes\n");
            else printf("No\n");
        }

    }


    return 0;
}

POJ - 3244 Difference between Triplets

这道题是一道推公式的题，提议就是就N个三元组之间的距离D的和，两个集合的D表示为

$$D(a,b) = \max ( I_a - I_b ,J_a - j_b ,K_a - K_b) - \min (I_a - I_b,J_a - -J_ b,K_a - K_b )$$
我们先设a = ia - ib,b = ja - jb,c = ka - kb，然后把他们画在数轴上，我们会发现其实求D也就办成了求这三个点之间最远的距离，也就是 $D =\frac {|a - b| + |b- c| +| c - a|}{2}$.然后我们把ia，ib带回去，然后就又能得出一个式子
$$D = \frac{|I_a - I_b - J_a + J_b|+|J_a - J_b - K_a + K_b| +|K_a - k_b - I_a +I _b|}{2}$$
这样的话，我们再设a = i - j;b = j - k;c = k - i;
为了去掉绝对值，我们把a，b，c排序，就可以算出他们的贡献。对于任意一个a[i]，在比他小的i个数中加了i次，比他大的n - i -1次中被减了n- i -1次，合并一下就是2i - n +1次。注意要开的大一点不然会爆

#include <cstdio>
#include <map>
#include <bitset>
#include <cstring>
#include <algorithm>
#define ll __int64
using namespace std;
const int maxn = 2e5 + 5;
int a[maxn],b[maxn],c[maxn];
int main()
{
    int n;
    while(scanf("%d",&n) != EOF)
    {
        if(n == 0) break;
        for(int i = 1;i<=n;i++)
        {
            int x,y,z;
            scanf("%d%d%d",&x,&y,&z);
            a[i] = x - y;
            b[i] = y - z;
            c[i] = z - x;
        }
        sort(a+1,a+n+1);
        sort(b+1,b+n+1);
        sort(c+1,c+n+1);
        ll ans = 0;
        for(int i =1;i<=n;i++)
        {
            ans += (ll)(i * 2 - n + 1)*a[i];
            ans += (ll)(i * 2 - n + 1)*b[i];
            ans += (ll)(i * 2 - n + 1)*c[i];
        }
        printf("%I64d\n",ans/2);
    }
    return 0;
}

POJ - 3685 Matrix
这道题的题意就是给出一个公式和参数的范围，求第M小的数。看到求第M小的数肯定是要2分了。我们首先能想到二分最后的答案，那么怎么写判断呢。首先这个公式是一个二元的函数，如果我们固定一个变量j，那么关于i的函数值是递增的，这个可以借助matlab画图看一下。我们枚举j，然后二分当前j中有多少个i对应的函数值，小于我们二分的函数值。最后把每个j的结果加起来就能拿去判断了。这个思路真的巧妙啊，学习一个~

#include <cstdio>
#include <iostream>
#include <climits>
#define ll long long
using namespace std;
const ll INF = 0x3f3f3f3f3f3f3f;    //数字可能很大
ll n,k;
ll calc(ll i,ll j)
{
    return i*i + j*j + 100000*(i - j) +i *j;
}
ll check(ll m)
{
    ll cnt = 0;
    ll ans = 0;
    for(ll j = 1;j<=n;j++)
    {
        ll l = 1,r =n;
        while(l <= r)
        {
            ll mid = (l + r) >> 1;
            if(calc(mid,j) <= m)
            {
                ans = mid;
                l = mid+1;
            }
            else
            {
                r = mid-1;
            }
        }
        cnt += ans;
    }
    return cnt;
}
int main()
{
    int ca;
    scanf("%d",&ca);
    while(ca--)
    {
        scanf("%I64d%I64d",&n,&k);
        ll l = -INF,r = INF;
        ll ans = 0;
        while(l <= r)
        {
            ll mid = (l + r) >> 1;
            if(check(mid) >= k)
            {
                ans = mid;
                r = mid-1;
            }
            else
            {
                l = mid +1;
            }

        }
        printf("%I64d\n",ans);
    }


    return 0;
}