Topcoder 酱油记

　　SRM 632 Div.2 (58/688):

　　这场是除了我第一场之外最简单的一场了，前两道太简单，一下就秒了。最后一题复杂度很微妙，感觉深搜会超时，结果貌似没超...看了一下别人的代码，觉得那些用map的人的代码很好，也学着写了一下。现在感觉map真是太有用了。

// BEGIN CUT HERE

// END CUT HERE
#line 5 "GoodSubset.cpp"

#include <string>
#include <vector>
#include <map>
#include <iostream>
using namespace std;
typedef long long LL;
#define foreach(it,v) for(__typeof((v).rbegin()) it=(v).rbegin();it!=(v).rend();it++)
const int size = 50001;

map<int, int> ans;
class GoodSubset {
    public:
    int numberOfSubsets(int goodValue, vector <int> d) {
        for (int i = 0; i < d.size(); i++)
        {
            foreach(it, ans)
            {
                LL y = it->first;
                if (y * d[i] >goodValue) continue;
                ans[y * d[i]] = (ans[y * d[i]] + it->second) % 1000000007ll;
            }
            ans[d[i]]++;
        }
        return ans[goodValue];
    }
};

　　SRM 633 Div.2 :

　　还是没升到Div.1...

　　1000是道简单的搜索题，被我想复杂了。注意到a*b=gcd(a,b)*lcm(a,b)，所以知道了a之后，我们就可以求出b了。所以对于每个连通分量，我们只需要枚举其中一个数的大小就行了。复杂度O(n*10000)。（这题我用最小最大值优化后挂了...所以对于竞赛来说，能不优化就不优化吧...）

　　做算法题不能只是注重算法，还是要分析所给数据范围。

// BEGIN CUT HERE

// END CUT HERE
#line 5 "GCDLCMEasy.cpp"

#include <string>
#include <vector>
#include <iostream>
using namespace std;
typedef long long LL;
typedef vector<int> vi;
#define pb push_back
#define foreach(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++)
#define rep(i,n) for(int i=0;i<(n);i++)
#define FOR(i,l,r) for (int i=(l);i<=(r);i++)
const int maxn = 501, maxm = 500001;

int x[maxn];
int g[maxn], l[maxn];
vi mat[maxn];

bool vis[maxn];
int dl[maxn], po;
void bl(int v, vector <int> A, vector <int> B)
{
    int l = 1, r = 1;
    dl[r] = v; vis[v] = 1;
    for (; l <= r; v = dl[++l])
        for (auto i : mat[v])
            if (!vis[A[i] ^ B[i] ^ v])
                dl[++r] = A[i] ^ B[i] ^ v, vis[dl[r]] = 1;
    po = r;
}
class GCDLCMEasy {
    public:
    int gcd(int a, int b)
    {
        int temp;
        while (b != 0) temp = b, b = a % b, a = temp;
        return a;
    }
    string possible(int n, vector <int> A, vector <int> B, vector <int> G, vector <int> L) {
        for (int i = 0; i < A.size(); i++)
            mat[A[i]].pb(i), mat[B[i]].pb(i);
        bool ans = 1;
        for (int u = 0, v; u < n; u++) if (!vis[u])
        {
            bl(u, A, B);
            for (x[u] = 1; x[u] <= 10000; x[u]++)
            {
                for (int z = 2; z <= po; z++)
                    x[dl[z]] = 0;
                bool exist = 1;
                for (int z = 1; z <= po; z++)
                {
                    v = dl[z];
                    for (auto i : mat[v])
                    {
                        int ww = A[i] ^ B[i] ^ v;
                        if (!x[ww]) x[ww] = G[i] * L[i] / x[v];
                        if (gcd(x[v], x[ww]) != G[i]
                            || x[v] * x[ww] / gcd(x[v], x[ww]) != L[i]) { exist = 0; break; }
                    }
                    if (!exist) break;
                }
                if (exist) break;
            }
            ans &= (x[u] != 10001);
        }
        return ans ? "Solution exists" : "Solution does not exist";
    }
};