codevs1735 方程的解数(meet in the middle)

题意

题目链接

Sol

把前一半放在左边,后一半放在右边

meet in the middle一波

统计答案的时候开始想的是hash,然而MLE了两个点

实际上只要排序之后双指针扫一遍就行了

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 7, MAX = 1e7 + 10;
int K[MAXN], P[MAXN], N, M, ans;
int a1[MAX], c1, a2[MAX], c2, cnt[MAX];
int fp(int a, int p) {
    int base = 1;
    while(p) {
        if(p & 1) base = base * a;
        a = a * a; p >>= 1;
    }
    return base;
}
void dfs(int x, int Lim, int opt, int sum) {
    if(x == Lim + 1) {
        if(!opt) a1[++c1] = sum;
        else a2[++c2] = -sum;
        return ;
    }
    for(int i = 1; i <= M; i++) dfs(x + 1, Lim, opt, sum + K[x] * fp(i, P[x]));
}
int main() {
    ios::sync_with_stdio(false);
    cin >> N >> M;
    for(int i = 1; i <= N; i++) cin >> K[i] >> P[i];
    if(N <= 2) {
        a1[++c1] = 0;
        dfs(1, N, 1, 0);
    } else {
        dfs(1, N / 2, 0, 0);
        dfs(N / 2 + 1, N, 1, 0);
    }
    sort(a1 + 1, a1 + c1 + 1);
    sort(a2 + 1, a2 + c2 + 1);
    int j = 1;
    for(int i = 1; i <= c2; i++) {
        if(i != 1 && (a2[i] == a2[i - 1])) {cnt[i] = cnt[i - 1]; continue;}
        while(a1[j] <= a2[i] && j <= c1) {
            if(a1[j] == a2[i]) cnt[i]++;
            j++;
        }
    }
    /*
    for(int i = 1; i <= c1; i++)
        for(int j = 1; j <= c2; j++)
            ans += (a1[i] == a2[j]);
    */
    for(int i = 1; i <= c2; i++) ans += cnt[i];
    cout << ans;
    return 0;
}

原文地址:https://www.cnblogs.com/zwfymqz/p/10246639.html

时间: 2024-10-10 06:47:12

codevs1735 方程的解数(meet in the middle)的相关文章

【poj1186】 方程的解数

http://poj.org/problem?id=1186 (题目链接) 题意:已知一个n元高次方程:  其中:x1, x2,…,xn是未知数,k1,k2,…,kn是系数,p1,p2,…pn是指数.且方程中的所有数均为整数. 假设未知数1 <= xi <= M, i=1,,,n,求这个方程的整数解的个数. solution  meet in the middle.移项,分两部分搜索,hash判断两次dfs的结果是否相同,统计结果. 代码: // poj1186 #include<alg

子集(状态压缩)(meet in the middle)

子集 [问题描述] R 君得到了?个集合,???共有 n 个正整数. R 君对这个集合很感兴趣. R 君通过努?钻研,发现了这个集合?共有 2n 个子集. 现在 R 君又对这个集合的?集很感兴趣. 定义?个集合的权值是这个集合内所有数字的和的话. 那么 R 君想问问你,这个集合的权值第 K小子集是多?. ps. 涉及到较少数字的 long long 输?输出,建议使用 cin/cout. [输入格式] 第??两个正整数 n,k. 接下来一行 n 个正整数,表?集合内元素. [输出格式] 输出?个

SPOJ4580 ABCDEF(meet in the middle)

题意 题目链接 Sol 发现abcdef是互不相关的 那么meet in the middle一下.先算出abc的,再算def的 注意d = 0的时候不合法(害我wa了两发..) #include<bits/stdc++.h> #define LL long long using namespace std; const int MAXN = 101, SS = 2e6 + 10; map<LL, LL> mp; int N; LL a[MAXN], ans; int a1[SS]

[CSP-S模拟测试]:毛一琛(meet in the middle)

题目描述 历史学考后,$MYC$和$ztr$对答案,发现选择题他们没有一道选的是一样的.最后他们都考了个$C$.现在问题来了,假设他们五五开,分数恰好一样(问答题分数也恰好一样,只考虑选择题).已知考题是$N$道选择题(第$i$题分数为$M(i)$).问$ztr$和$MYC$做对的题的并有多少种可能?众所周知,历史学考选择题有$25$题,但是$MYC$为了给你降低难度,$n$不超过$20$. 一句话题意:有多少个非空子集,能划分成和相等的两份. 原题见:$USACO\ 2012\ OPEN\ G

「10.13」毛一琛(meet in the middle)&#183;毛二琛(DP)&#183;毛三琛(二分+随机化???)

A. 毛一琛 考虑到直接枚举的话时间复杂度很高,我们运用$meet\ in\ the\ middle$的思想 一般这种思想看似主要用在搜索这类算法中 发现直接枚举时间复杂度过高考虑枚举一半另一半通过其他算法统计,保证两边互不影响 今天的题我们考虑枚举先枚举左半部分,然后每个物品有三种取值情况 选入A集合,选入B集合,不选,系数不同 考虑完左半部分再去考虑右半部分,那么我们可以用哈系表先从将左半部分的答案统计出来 然后右半部分查询他的相反数注意去重 也可以用将两边状态都用结构体存下来 注意去重 思

折半搜索(meet in the middle)

折半搜索(meet in the middle) ? 我们经常会遇见一些暴力枚举的题目,但是由于时间复杂度太过庞大不得不放弃. ? 由于子树分支是指数性增长,所以我们考虑将其折半优化; 前言 ? 这个知识点曾经在模拟赛中出现过,所以这里稍微提一下; ? 讲的很浅显,但是不要D讲者; 入门 ? dfs搜索树是指数性增长,如果将指数减少一半,就将会有量的飞跃,所以在遇见暴力枚举太大时,我们可以考虑这种算法; ? 总体思想即,dfs搜素通常从一个点出发,遍历所有深度,那么我们考虑将深度减半,从两个点出

【gym102222K】Vertex Covers(高维前缀和,meet in the middle)

题意:给定一张n点m边的图,点带点权,定义点覆盖的权值为点权之积,问所有点覆盖的权值之和膜q n<=36, 1<=a[i]<=1e9,1e8<=q<=1e9 思路:n<=36,考虑middle in the middle分成两个点数接近的点集L和R 对于L,枚举其子集S,判断S能否覆盖所有L内部的边,预处理出所有合法的S的超集的贡献 对于R,枚举其子集T,判断T能否覆盖所有R内部的边,如果可以则可以推出L,R之间在确定R中选T的前提下左边至少需要选点集T’,答案即为T的

方程的解数(codevs_1735)——hash

这题挺简单的吧,其实只需要一步--就可以把暴力搜索的m^6降到m^3--那就是--请看代码 #include<iostream> #include<cstdio> #include<cstring> using namespace std; inline int read(){ int t=1,num=0;char c=getchar(); while(c>'9'||c<'0'){if(c=='-')t=-1;c=getchar();} while(c>

bzoj 1770 [Usaco2009 Nov]lights 燈 meet in the middle

题面 题目传送门 解法 题解里都是高斯消元然后dfs,蒟蒻表示不会 直接分两半dfs即可 时间复杂度:\(O(2^{\frac{n}{2}})\) 代码 #include <bits/stdc++.h> #define LL long long #define N 110 using namespace std; template <typename node> void chkmax(node &x, node y) {x = max(x, y);} template &